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1 Introduction 



Quantum Chromodynamics (QCD) is today considered the fundamental theory of 
strong interactions. The adimensional couphng of the theory Qq, once it is renor- 
mahzed, will depend on the energy scale of the considered physical process, and it 
gives a measure of the strength of the interaction at that energy scale. The out- 
standing property of QCD, asymptotic freedom, tells us that the coupling decreases 
with increasing energies. This allows the usage of perturbation theory to make phe- 
nomenological predictions for processes with large momentum transfered. However 
the increase of the coupling with decreasing energies does not allow to use pertur- 
bative methods to compute physical quantities in the low energy region such as 
the mass spectrum or hadronic matrix elements. A possible strategy in this case 
is to use a non-perturbative regularization of the theory, introducing a discretized 
space-time (lattice) [1,2]. This strategy has two main advantages: first it provides 
an ultraviolet regularization and secondly it reduces the degrees of freedom of the 
theory to a numerable infinity. Considering the theory in a finite volume lattice, it is 
possible to perform numerical simulations of QCD through Monte Carlo methods. 
The continuum QCD action has to be discretized in a sensible way, and a simple 
and attractive lattice action is the one proposed by Wilson long time ago [1,2]. 

The relation between the coupling constant 5^0 and the lattice spacing a in physical 
units is given, in lattice QCD, by the renormalization group equations, and we are 
naturally interested in the continuum limit of the theory. The simulations of lattice 
QCD are always performed for values of go corresponding, in the renormalized theory, 
to a finite and non zero value of a. This introduces in the results of the simulations, 
using the Wilson fermion action, errors of order a and numerically these 0(a) errors 
could be of the order of 20 — 30%. The 0(a) discretization errors could be eliminated, 
in principle, following the Symanzik's improvement program [3-5], where the 0(a) 
cutoff effects in on-shell quantities are canceled by adding local 0(a) counterterms to 
the lattice action and to the composite fields of interest [6-10]. A technical difficulty 
is that the improvement coefficients multiplying these counterterms are not known 
a priori, and they should be all computed using Monte Carlo simulations. 

A new intriguing possibility is the so called automatic 0(a) improvement [11], where 
none of the improvement coefficients are needed in order to have O(a^) cutoff effects 
in physical observables. The basic idea is that the Wilson theory for fermions with 
a suitable infrared cutoff is in the massless limit free from 0(a) errors. We will see 
that to extend automatic 0(a) improvement to a theory in infinite volume with a 
non zero mass term we have to add the so called twisted mass, keeping the standard 
quark mass to be zero. The twisted mass term, that in a way will act also as a sharp 
infrared cutoff, can be obtained in continuum QCD via a non-anomalous chiral 
rotation from the standard mass term. To be specific, if we consider QCD with a 
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field X describing a flavour doublet, the twisted mass term looks like 



(1.1) 



where is the Pauli matrix in flavour space and /iq is what is called twisted mass. 

The twisted mass term in a lattice QCD action appears to my knowledge for the 
flrst time in ref. [12], where it is given an ansatz for the phase structure of Wilson 
fermions in the parameter space rriQ — gQ, where ttiq is the bare quark mass. Based on 
the analysis of the lattice Gross-Neveu model, and on the strong coupling expansion 
of Wilson lattice QCD, the author suggested that there are regions in the parameter 
space of Qq and mo where the true vacuum has a non zero expectation value of 
iXl5X signalling the spontaneous breaking of parity symmetry. It was then natural 
to propose, in order to pick up the real vacuum from numerical simulations, to 
add an external fleld iHtp'-f^ip to the original lagrangian and to perform the limit 
i7 — > 0"*". The H is what now we would call twisted mass. The twisted mass in this 
case is then just an external fleld used to probe the structure of the vacuum of the 
theory and it has to be removed at the end of the computation. We will come back 
in deatil to the chiral phase structure of the Wilson theory in sect. 5. Here I just 
would like to mention that in the same paper a new method to improve the scaling 
behaviour of the chiral condensate was proposed based on the observation that the 
scaling violations of the condensate are odd under a change of sign of the coefficient 
of the Wilson term, and they can be easily averaged out. We will see in app. E that 
this is a possible starting point to understand automatic 0(a) improvement. 

The twisted mass term breaks parity and flavour symmetry. It is a natural question 
whether this mass term changes also the continuum action or just the discretization 
errors of the theory. In this report we will show that with a Wilson fermion lattice 
action the twisted mass term generates parity and flavour violating cutoff effects (in 
most cases of 0(a^)) which go away performing the continuum limit. 

The fact that the twisted mass actually induces only flavour and parity breaking 
cutoff effects, and it is actually equivalent to QCD, can be understood considering 
an old remark made by Gasser and Leutwyler. In fact it was noticed many years 
ago [13] that the usage of what we would now call a twisted mass term is irrelevant 
in continuum QCD. The fermionic part of the 2 flavours QCD Lagrangian is usually 
given in the form 



where m is the mass matrix. The quark masses in the standard model originate from 
the asymmetries of the electroweak vacuum. Since the electroweak interactions do 
not preserve parity there is no reason a priori for the quark mass term of QCD to be 
parity invariant, and can be generically written as x(m + ^75M)x- We assume here 
that /X is a traceless matrix to avoid an unnecessary discussion of the QCD vacuum 
angle. With a suitable non-anomalous chiral transformation of the quark flelds the 



>Cqcd = X It^D^ + m X 



(1.2) 
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general mass term x(m + i75/x)x can always be brought to the standard form, where 
the mass matrix is diagonal with real positive eigenvalues and m^. The remaining 
part of the Lagrangian constrained by the requirement of renormalizability is left 
invariant by this chiral transformation. In brief, with a change of variables in the 
functional integral, that leaves the measure invariant, one can show that even if 
a general parity and flavour violating mass term is allowed, the request of having 
a renormalizable theory preserving gauge and Lorentz invariance, generates these 
"accidental symmetries". This simple example shows that the specific form of the 
mass term in the 2 flavours continuum QCD Lagrangian is actually irrelevant. The 
reason for this is the fact that the massless theory is invariant under the chiral 
non-anomalous transformation that changes the form of the mass term. 

This is just an example of a more general phenomenon. Renormalizable theories 
that describe electro-weak and strong interactions, can be considered as low energy 
effective theories of more general not necessarily renormalizable high energy the- 
ories. The condition of having low energy renormalizable field theories can be so 
stringent that the corresponding Lagrangian may turn out to obey extra accidental 
symmetries, that were not symmetries of the higher energy theories [14]. 

This observation becomes important on the lattice. If we discretize the continuum 
QCD action with Wilson fermions [1], the Wilson term explicitly breaks chiral sym- 
metry and the lattice action is not invariant anymore under the field rotations men- 
tioned before. But we still have the freedom to choose the Wilson term and the mass 
term to point in different relative "directions" in the Dirac and flavour space [15]. 
This freedom is the key to constrain the form of the cutoff effects induced by the 
Wilson term. 

The observation that physical observables computed with the Wilson lattice action 
are automatically 0(a) improved in the "infrared safe" (i.e. with no spontaneous 
symmetry breaking) chiral limit is relevant also for the renormalization properties 
of local operators that depend on the breaking of chiral symmetry induced by the 
Wilson term. 

To summarize: Wilson twisted mass QCD is a lattice regularization that allows 
automatic 0(a) improvement only tuning one parameter. The bare untwisted quark 
mass mo has to be tuned to the so called critical mass in order to maximally disalign 
the Wilson term and the mass term. In this approach the renormalization of local 
operators relevant for phenomenological applications is significantly simplified with 
respect to the standard Wilson regularization. The price to pay is the existence of 
O(a^) cutoff effects that break parity and flavour symmetry. All these statements 
will be demonstrated explicitly in this report. 

This is not the only review on twisted mass QCD. A set of lectures has been pre- 
sented by S.Sint at the School "Perspectives in Lattice Gauge Theories" [16]. In 
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these lectures a nice introduction on the basic setup, exceptional configurations, 
automatic 0(a) improvement together with few applications of Wtm is given. Par- 
ticular emphasis is also put in finite volume renormalization schemes with chirally 
twisted bounday conditions. In our report we enlarge the topics covered by Sint, and 
we elaborate on the ones already there. On the other side we touch only marginally 
finite volume renormalization schemes with chirally twisted bounday conditions. For 
this reason we believe that the present review is in many respects complementary 
to the one of Sint, and together they can be used as a complete introduction to all 
the topics connected with twisted mass QCD. 

The paper is organized as follows. In section 2 I use classical considerations in 
continuum QCD to show the equivalence of twisted mass (tm) QCD and QCD. I 
also describe the rigorous theoretical properties of Wilson twisted mass QCD (Wtm 
QCD), for degenerate and non-degenerate quarks. In section 3 I discuss the 0(a) 
discretization errors of Wilson-like lattice actions. In particular I show several proof 
of automatic 0(a) improvement, with a particular emphasis on the choice of the 
critical mass. Numerical results confirming this property will conclude the section. 
In section 4 I derive again some of the results obtained in the previous sections 
using a different fermion basis. Hopefully this could help the reader in a better 
understanding of the subject of this review. In section 5 I analyse O(a^) parity 
and isospin violating discretization errors and O(a^) cutoff effects responsible for 
the non trivial chiral phase structure of the lattice theory. In section 6 I discuss 
selected numerical results obtained with Wtm QCD and show different methods to 
ease the renormalization of local operators. In section 7 I make a short digression on 
algorithms to simulate light Wilson-like quarks. This is an important prerequisite 
for part of the numerical results presented in this review. Conventions, notations 
and more technical discussions are deferred to the Appendix. 
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2 Basic properties 



In this section we introduce twisted mass QCD in the continuum using classical 
arguments for a doublet of degenerate quarks. This academic exercise allows the 
reader to get acquainted with twisted mass QCD and to learn how to relate cor- 
relation functions from QCD to twisted mass QCD. To extend this concepts at a 
quantum level we discretize the theory on a lattice using Wilson fermions. The re- 
sulting theory Wilson twisted mass (Wtm) QCD is ultra-local, unitary, reflection 
positive and renormalizable to all orders in perturbation theory. 

2.1 Twisted mass QCD in the continuum 

We first consider the continuum limit of the twisted mass (tm) QCD action for 
N{ = 2 degenerate quarks. We will always work in euclidean space and all the 
definitions and conventions are collected in app. A. We call the set of fermion fields 
{Xi x} the twisted basis. In this basis the tm QCD action reads 

Sf[x,X, G]= j ^^^X {it^D^ + mq + i/iq75r^) X, (2-1) 

where = c?^ + denotes the covariant derivative in a given gauge field G^, 
and is the third Pauli matrix acting in flavour space. To better understand the 
structure of the action we write explicitly all the indices. The fermion fields are 

XA,aA^) (2.2) 

where A is the colour index, a the Dirac index and i the flavour index. The action 
then reads 

Sf[X,X,G] = J (i'^XXA,a,ii^'> {{ll^)ap{Df,)AB5ij + {m^)5AB5ap5ij 

+ {^^^q){'lb)o,|3{r^)ij^AB) XB,f3,j{x). (2.3) 

In the following we will use the compact notation in order to keep the reading not 
too heavy. The mass term of the tm action can be written as 

+ ifi^jnT^ = Me*"^^"' (2.4) 

where 

M = ^ml + /i2 (2.5) 

is the so called polar mass. At the moment the tmQCD action is just a rewriting of 
standard QCD in a different basis. In fact performing the following axial transfor- 
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mat ion 



the form of the action is left invariant, but the mass term transforms into 



(2.6) 



(2.7) 



In particular, the standard QCD action for Nf = 2 degenerate quarks 



/ 



(2.8) 



is obtained if u; = a , i.e. if the twist angle uj satisfies the relation 



tancu = fiq/niq. 



(2.9) 



We call the physical basis {ip, ip} the basis where the continuum QCD action takes 
the standard form. The two basis are related by the rotation (2.6) with uj satisfying 



In the following we will mainly use the twisted basis since this is the basis used 
in the numerical simulations. Although the physical interpretation of the fermionic 
correlation functions is most transparent in the physical basis the renormalization 
of gauge-invariant correlation functions, including those with insertions of local 
operators, looks simpler in the twisted basis. This will become clear later when we 
will discretize the action (2.1) using Wilson fermions [1,2]. In sect. 5 I derive again 
some results in the physical basis. 

Twisted mass QCD (2.1) and standard QCD (2.8) actions are exactly related by 
the transformation (2.6) and therefore share all the symmetries. The symmetry 
transformations in the twisted basis are simply the transcription of the standard 
symmetry transformations using eq. (2.6). These symmetry transformations will be 
from now on called twisted symmetries and they are collected in app. B. The twisted 
vector symmetry SUy(2)i^ defined in eq. (B.l) is a symmetry of the tmQCD action 
(2.1), while the mass term M breaks the twisted axial symmetry SUa{2)^^ defined in 
eq. (B.2). This is equivalent to the transformation properties of the QCD action (2.8) 
under the standard SUy{2) and SUx{2) symmetries transformations. 

In order to give the correct physical interpretation of the interpolating fields in the 
twisted basis it is important to have available the relations between the currents in 
the two basis. The axial and vector currents are 



eq. (2.9). 




^ Xl^l 2 



(2.10) 



while the pseudoscalar and scalar densities are defined by 



P' 



<a 





(2.11) 
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The axial transformation (2.6) of the quark and anti-quark fields induces a trans- 
formation of the composite fields. We indicate with a calligraphic symbol the corre- 
sponding currents in the physical basis. For example, the rotated axial and vector 
currents read, 



cos(cj)A;'^ + e^"'' sm{uj)V^ (a = 1, 2) 
Ai (a = 3), 



ya^^r_^ = r°^W + ^'"'^i^(^K (« = 1'2), 
' '2 [V^ (« = 3), ^ • ^ 

and similarly, the rotated pseudo-scalar and scalar densities are given by 

P^-^75-^= 1'^" = (2 14) 

^^'2^ [cos(u;)p3 + ,sin(cj)i^o (a = 3), ^' ' 

=^^ = cos(cj)S° + 2i sm{iu)P^. (2.15) 

The same procedure to change basis can be used for any local operator. In app. F 
we give the example for a proton interpolating field. 

The form of the Ward identities in the twisted basis is slightly different from the 
standard form. The local SUy{2) x SU\{2) chiral transformations of the fermionic 
fields are defined as follows 

6x{x) = ^[a^(x)y + «^(x)y75]x(a;) (2.16) 

Sx{x) =tx{x)[-a1.{x)^ + a'Xix)^-f5\, (2.17) 

and the symmetry at the classical level {SS = where S indicates the variations on 
the fermion fields of eqs. 2.16 - 2.17) gives the so-called partially conserved axial 
current (PCAC) and partially conserved vector current (PCVC) relations 

a^A;^ = 2mqP'^ + 2/iq53'^^o^ ^2.18) 

d,V^ = -2fi^e'''''P'. (2.19) 

It is easy to verify that the rotated currents and densities satisfy the Ward identities 
in their standard form, 

d^A; = 2MV\ 9^V; = 0, (2.20) 

if uo is related to the mass parameters as in eq. (2.9). The PCAC and PCVC masses 
that appear in eqs. (2.18,2.19) are the 2 components of the physical mass which is 
given by the polar mass M (see eqs. 2.5 and 2.20). A particular case is when one of 
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the 2 masses vanishes: then the physical mass is given by the non-vanishing mass. 
We will see in the following that a very interesting case is when we work at full twist, 
also called maximal twist. Full twist corresponds at the classical level to mq = or 
equivalently to = 7r/2. In this case the role of the physical mass is fully played by 
the twisted mass yUq. 

2.2 Beyond the classical theory 

These classical considerations, based on the possibility of performing the axial trans- 
formation (2.6) on the fermion fields, show that there is a one-to-one correspondence 
between tmQCD and standard QCD. 

Denoting the sum of the gauge and fermion QCD actions by 5* = 5'g + Sp, the 
physical content of the quantum theory can be extracted from its n-point correlation 
functions 

. . . , Xn)) = J D[^, ^D[U] e-^Oix,, . . . , (2.21) 

with 

Z = l D[x,x]D[U]^'^, (2.22) 

where 0{xi, . . . product of local gauge invariant composite fields which are 

localised at the space time points xi, . . . ,Xn- For the discussion of this subsection we 
concentrate on the dependence of the correlation functions on the fermionic fields 
and on the mass. Formally in the functional integral (2.21) we can make the axial 
change of variables (2.6). This change of variables is non anomalous, i.e. it does 
not change the integration measure. As we have already discussed in the previous 
subsection only the mass term of the action is affected. Then the relation between 
correlation functions in the two basis reads 

W,i^])iM,o) = {0[x,x])i^m,,,,) (2.23) 

where the index of the correlation function in the l.h.s indicates that it has been 
computed in standard QCD with quark mass M and the index of the correlation 
function in the r.h.s indicates that it has been computed in tmQCD with quark 
masses and yUq. The relation between the arguments of the correlation function 
is given exactly by eq. (2.6). An explicit example of eq. (2.23) reads 

{Alix)V\y))^M,o) = cosu{Alix)P\y))(^m.^,,^) + smu{V^ix)P\y))^^^,,^y (2.24) 

Correlation functions in QCD can be written as a linear combination of correlation 
functions computed in tmQCD at a given twist uj. 
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This equivalence remains valid at finite lattice spacing for Wilson fermions up to 
discretization errors, if the theory is correctly renormalized in a mass independent 
scheme [15]. 

To carry over these formal considerations to a rigorous level, we have to regularize 
the theory with a regulator which preserves the chiral symmetry of the massless 
theory, i.e. which preserves the axial symmetry containing the transformation (2.6). 
A chiral invariant regularization is provided by Ginsparg- Wilson (GW) quarks [17] 
on the lattice. 

With GW fermions we can repeat the same steps performed formally in the contin- 
uum theory. In particular with GW fermions eq. (2.23) is valid in the bare theory, 
and analogously the bare Ward identity masses coincide with the bare mass param- 
eters of the action. Using then Ginsparg- Wilson fermions as a theoretical tool it is 
possible to prove [15] that tniQCD and QCD are equivalent, i.e. given a lattice reg- 
ularization that preserves chiral symmetry, the one-to-one correspondence between 
tmQCD and standard QCD is preserved at finite lattice spacing: consequently they 
have the same continuum limit. There is only one condition that has to be fulfilled 
for this statement to be true: all multiplicative renormalization constants have to be 
independent of the angle u, not only up to cutoff effects [15]. A simple example of 
a suitable renormalization scheme of this kind is a mass-independent scheme which 
is obtained by renormalizing the theory in the chiral limit. Under this condition 
then the relations between renormalized correlation functions take the analogous 
form given by classical considerations. In other words in the continuum or with GW 
fermions, there is no reason to introduce a twisted mass term: it can be rotated 
away by a change of variable in the functional integral, because in both cases the 
massless theories are chirally symmetric. 

Based on universality one then expects that tmQCD in other, not necessarily chirally 
invariant, regularizations can again be renormalized such that equivalence between 
renormalized correlation functions computed in tmQCD and standard QCD is sat- 
isfied up to cutoff effects [15]. In particular if the regulated theory breaks chiral 
symmetry even in the massless limit the twisted mass term cannot be rotated away 
and the twist angle u parametrizes a family of regularizations, which differ at finite 
lattice spacing, but have the same continuum limit. It may hence be advantageous 
to choose the twist angle in a suitable way, e.g. to reduce discretization errors. 

The relations between local fields in the classical continuum theory, e.g. eqs. (2.12,2.13), 
thanks to the equivalence between tmQCD and QCD at the quantum level, can be 
extended to the renormalized theory. Thus the classical theory may be used as a 
guide to establish the relations between renormalized correlation functions. This is 
a very important result and we will come back on it in sect. 2.8, in order to better 
understand which are the relations between tmQCD and QCD correlation functions. 
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2.3 Wilson twisted mass QCD 



We are interested to exploit the freedom given to us by the choice of the twist angle 
in order to improve the properties of chirally non-invariant Wilson regularization of 
lattice QCD. 

The setup is a hypercubic infinite lattice, with lattice spacing a. Standard reference 
books for lattice field theories are [18-22]. The gauge group is SU (iVc) and the gauge 
field on the lattice is an SU (N^) matrix U{x; fi) that depends on the lattice point 
and on the four directions. Fermion fields reside on the lattice sites and as we have 
explained in the previous section, they carry colour, Dirac and flavour indices. The 
lattice action for A^f = 2 degenerate flavours is of the form 

S[x,X, U] = Sg[U] + Sf[x,X, U], (2.25) 
where So denotes the Wilson plaquette action 

SG[U] = f^ E tr{l-pixi(x;/x,z/)}, (2.26) 

with P = 6/gQ, go being the bare gauge coupling and P^^^(x; /i, z/) the parallel 
transporter around a plaquette P 

P^^^(x; /i, 1^) = U{x; fL)U{x + afi; u)U{x + aP; fiy^U{x; u)-\ (2.27) 

The sum runs over all the oriented plaquettes P on the lattice. With Wilson quarks 
the tmQCD Dirac operator is given by 



X 



Xix), (2.28) 



where 



1 



= + V;) - arV;V J, (2.29) 



V^, V* are the standard gauge covariant forward and backward derivatives defined 
in app. A, mo and fiq are respectively the bare untwisted and twisted quark masses. 
The parameter r is the so called Wilson parameter and it is always set to 1 unless 
specified. We will call the action (2.28) the Wilson twisted mass (Wtm) QCD action. 
Sometimes Wilson actions are also written as 
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a; 



1 + z2Ka/iq75r^ x{x) 
+ U{x - fl; +-f^)x{x- ft) 



At=0 



where we define the rescaled dimensionless fermion field 



X 



,3/2 ^' 



X II ^ dXi, 



and the hopping parameter 



(2.30) 



(2.31) 



(2.32) 



8 + 2amo 

The hopping parameter is an alternative way to label the bare untwisted quark mass, 
as f3 defined before is an alternative way to label the bare gauge coupling. 



The first term of eq. (2.29) is a standard symmetric discretization of a lattice deriva- 
tive and the second term 

axV;V,x (2.33) 

is the so-called Wilson term. It can be easily checked that this term is not invariant 
under all the axial transformations, but it is needed in order to remove from the 
spectrum of the theory other spurious particles, called doublers, introduced by the 
lattice discretization. In the case of naive fermions, without the Wilson term, the 
existence of the doublers is related to the so called spectrum doubling symmetry [23] , 
that can be seen as an exchange symmetry among the corners of the Brillouin zone 
in the reciprocal momentum space. It can be shown that the doubling phenomenon 
is a more general feature of ultra-locaf^ actions that goes under the name of Nielsen- 
Ninomiya theorem [24,25]. We will not go further into this topic, but since we are 
here interested in an action with a next-neighbour interaction (i.e. ultra-local) we 
need to insert a Wilson term. 



The Wilson term (2.33) breaks explicitly the axial symmetry (2.6). This implies 
that the twisted mass term cannot be rotated away by a chiral transformation and 
the exact equivalence between the Wilson action with vanishing and non-vanishing 
twisted mass is lost: Wilson and Wilson twisted mass are different lattice regular- 
ization. As we have discussed in sect. 2.2 the exact equivalence is recovered only in 
the continuum limit. 



With ultra-local actions we think of actions where the interaction range is spread over 
a finite number of points of the lattice. 
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2-4 Symmetries 



The introduction of a twisted mass and Wilson terms requires an analysis of sym- 
metries. In sect. 2.1 we have already seen that at the classical level the symmetries 
of tmQCD are the transcriptions of the standard QCD symmetry transformations 
via the axial transformation (2.6). They are collected in app. B. In the regulated 
theory, the action given in (2.28) breaks some of the symmetries of the classical 
action (2.1). The Wilson term (2.33) breaks twisted parity V^j (B.3), twisted time 
reversal 7^ (B.4), and twisted vector symmetry SUy{2)^ (B.l). Wtm shares with 
standard Wilson fermions the following symmetries: gauge invariance, lattice rota- 
tions, translations and charge conjugation C (see app. B for the definition). The 
ordinary parity symmetry transformation 

{f/(a;o,x; 0) ^ f/(xo, — x; 0), U{xq,x] k) ^ U^^{xo, —x — ak; k), /c = 1, 2, 3 
x(xo,x) 7oX(a;o, -x) 
x(a;o,x) x{xo, -x)7o 

(2.34) 

is only a symmetry if combined either with a discrete flavour rotation 

{f/(a;o, x; 0) ^ f/(xo, — x; 0), f/(xo, x; fc) ^ [/~^(a;o, — x — afc; fc). A; = 1, 2, 3 
x(xo,x) hoTi^2Xi{xo, -x) 
x(xo,x) -ixi^o, -x)ri,27o 

(2.35) 

or with a sign change of the twisted mass term 

V = Vx[fi^^ -/iq] . (2.36) 
The same holds for ordinary time-reversal 

{f/(xo, x; 0) ^ [/^^(— Xo — a, x; 0), f/(a;o,x; /c) ^ f/(— a;o,x; /c), /c = 1, 2, 3 
x(xo,x) i7o75x(-a;o,x) 
x(a;o,x) -x(-a;o,x) 27570 

(2.37) 

which is only a symmetry if combined either with a discrete flavour rotation 



^1,2 



[/(xo, x; 0) ^ f/ "'^(— Xo — a, x; 0), f/(xo, x; fc) ^ f/(— xq, x; A;), A; = 1, 2, 3 
x(xo,x) i7o75n,2X«(-a^o,x) 
^x(xo,x) -ix(-Xo,x)ri,2757o 

(2.38) 

or with a sign change of the twisted mass term 

r = r X K ^ -^q] . (2.39) 
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Incidentally this implies that CVT is a symmetry of the lattice action (2.28). Con- 
cerning the continuous symmetries the ordinary isovector SUy(2) symmetry is bro- 
ken explicitly by the fiq term down to the 1/^(1)3 subgroup with diagonal generator 
T3 

which is a symmetry of the lattice action (2.28) together with the U{1) transforma- 
tions associated with fermion number conservation. 



2.4-1 Chiral symmetry at lu = ^ 

In this section we analyze the case of full twist, i.e. u = 7r/2, because for axial 
and vector transformations this is a special case. We will discuss in sect. 2.8 and 
extensively in sect. 4.3 how to tune, in a non-perturbative way, the twist angle u, 
i.e. the bare untwisted quark mass ttiq. For what follows it is enough to assume that 
mo, i.e. LJ, has been tuned appropriately. 

The twisted "charged" axial transformations for a = 1, 2 and = 7r/2 read 



1,2 




i-f5T^) exp ( i^75T^'M 7f (1 + ^751-^) Xix 



Xix)^ (1 + ^15T') exp(i^75ri'2)-L (1 



hence 



1,2 




exp ±i^-T ' I xix) 



1,2 

2,1^ 



x{x) exp(=Fz^r 



The twisted "charged" vector transformations for a 



X{x) — ^ exp ( ±2-|-75r2'^ 



1 , 2 read 
X{x) 



1,2 



Xix) 



X{x) exp(±i^75r 



2,1^ 



(2.41) 



(2.42) 



For = f and for a = 1, 2 the form of the vector and axial transformations is re- 
versed compared with the ordinary transformations. In the continuum this is really 
just a different transcription of the chiral transformations, coming from a differ- 
ent choice of the fermionic basis. On the contrary in the regulated theory different 
terms of the action break different sectors of the axial and vector transformations. 
In particular we observe that the Wilson term (2.33) breaks the "charged" twisted 
vector symmetry (2.42) and it is invariant under the "charged" twisted axial sym- 
metry transformation (2.41). The twisted mass term has effectively an orthogonal 
behaviour because it is invariant under the "charged" twisted vector symmetry trans- 
formation (2.42) and it breaks the "charged" twisted axial symmetry (2.41). 
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To be more specific if we set tlie twisted mass to zero tlie Wilson theory is invariant 
under tlie group 



(2.43) 



while the twisted mass term is invariant under the group 



SUy{2)= f/v(l)z 



f/v(l). 



(2.44) 



This means that for the "charged" axial transformations the lattice action (2.28)0 



1,2 



has a continuum-like behaviour: the twisted "charged" axial symmetry t/A(l) 

is only softly broken by the mass term. This exact symmetry of the massless theory 
protects the charged pion from chiral breaking cutoff effects (see eq. 4.74). This 
result is obviously independent on the choice of the basis and we will discuss it 
again in sect. 5. A consequence of this consideration is also that the charged vector 
current in the twisted basis is protected from renormalization and it is at the same 
time the current that defines the pseudoscalar decay constant for the charged pion 
(see eqs. 2.61,7.9,7.10). The neat result being that at full twist the pseudoscalar 
decay constant for the charged pion does not need to be renormalized (see sect. 7). 

This analysis also shows in which sense the Wilson and the mass term at full twist 
are maximally disaligned: they are maximally disaligned concerning chiral symme- 
try. While the Wilson term breaks twisted "charged" flavour symmetry (2.42) the 
mass term breaks as in continuum QCD the full axial group (B.2), or to phrase it 
differently the Wilson term preserves a subgroup of twisted axial symmetry (2.41) 
while the mass term does not. What is relevant is that the Wilson and the mass 
term are "orthogonal" concerning the "charged" subgroup of chiral symmetry and 
this is achieved at full twist u = 7r/2. 



2.5 Tree-level 



To get some more insights on Wtm and to understand better the importance of 
disaligning mass term and Wilson term we compute here the tree-level Wtm prop- 
agator, that can be written as an integral over the first Brillouin zone of 

g . (2.45) 

where we have defined 

1 7^ 2 (XT) 

Pf, = - sin(ap^), M{p) = mo + -apl, P^, = - sin(— ^). (2.46) 

Qi Z: Qj 



I acknowledge a very useful discussion with G.C. Rossi on this point. 



15 



The poles of the propagator give us the spectrum of the theory. If we now make an 
expansion for small lattice spacing a neglecting all the terms of O(a^) we obtain 

+ mQ + arriQrp'^ + fi'^, where = p^p^. (2.47) 

The leading 0(a) discretization effects of the dispersion relation are given by the 
term arriQrp^. We can already make some interesting remarks. This term vanishes if 
we set mo = 0. This means that the chiral limit of the plain Wilson theory (/iq = 0) 
does not have 0(a) effects. Of course we are only considering the tree- level and if 
we would switch on the interaction the dynamics of the theory could change this 
result. We will see that if we consider the theory in a finite volume with suitable 
boundary conditions, where the mass dependence of the theory is smooth and there 
are no phase transitions, this result is still true. 

From eq. (2.47) we also see that even if we set mo = we can add a mass to the 
theory without introducing 0(a) effects. To understand how this can happen it is 
good to understand the origin of the 0(a) term amorp^. It comes from the cross term 
between the Wilson term and the mass term in M. (p) . This cross term is absent with 
a twisted mass because twisted mass and Wilson term point in different "directions" 
in the chiral-flavour space. 

Even if mo does not vanish but mo = 0(a) this will not change our conclusion since 
the term with mo will only change the O(a^) terms. 

All these considerations are only at tree- level, but we will see in sect. 4 that they 
remain true if we switch on the interaction between quarks and gluons. In particular 
setting the Wtm action at full twist allows to have physical observables that are 
automatically 0(a) improved. 

2.6 Transfer matrix 

Necessary and sufficient conditions under which the physical content of the theory 
in Minkowski space can be reconstructed from Euclidean Green's functions are the 
so called Osterwalder-Schrader conditions [26,27]. One of the required properties 
that does not obviously hold in a lattice theory is physical positivity. This condition 
states that given a gauge invariant polynomial of positive time (xq > 0) fundamental 
fields O one should have 

(6(0^)0) > (2.48) 

where G denotes euclidean time reflection and is the Hermitian conjugate of O. 
An explicit expression for the transfer matrix that in particular is strictly positive, 
i.e. all its eigenvalues are bigger than zero, was given for Wilson fermion and gauge 
actions in ref. [28]. This allows to prove the positivity condition (2.48) for the Wilson 
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action at finite lattice spacing. 



In the app. C we briefly repeat the steps of the proof with the extension to Wtm, 
because they become important for non-degenerate quarks. Here we simply list the 
main result: adding a twisted mass term for degenerate quarks does not change [29] 
the constraint on k (|k| < |) valid [28] for Wilson fermions. 



2. 7 Renormalization 



In perturbation theory, it has been shown that Wilson lattice QCD is renormal- 
izable [30-34], and we shall assume that this remains true beyond perturbation 
theory. Since the twisted mass term can be viewed as a super-renormalizable inter- 
action term which does not modify the power counting, this result can be extended 
also to Wtm. 

To understand the structure of the counterterms we use the symmetries of the Wtm 
lattice action (2.28), treating the masses as spurion flelds which transform under 
these symmetries. The counterterms to the action with dimension less or equal four 
are 

tr{F^^F^^}, XX, moxx, ^l^qXl^r^X, (2-49) 
where F^jy is the gluon fleld strength tensor. The flrst counterterm gives a multi- 
plicative renormalization of the bare gauge coupling. The others enter in the renor- 
malization of the bare quark masses. The continuum renormalized quark action can 
then be written as 



So = Sg[A] + I d^xxix) 



X(x), (2.50) 



where Sq now is the continuum gluon action. The renormalized parameters are given 
by 

9l = 9oZg{gl,afi), (2.51) 

mR = mqZ^{gl, a/i), (2.52) 

fiK = fiqZ^{g^,a^), (2.53) 

where fi denotes the renormalization scale dependence of the renormalization con- 
stants Z, and 

rrig = tuq — rricT- (2.54) 

It is well known that due to the loss of chiral symmetry at flnite lattice spacing, the 
bare untwisted quark mass renormalizes also additively, with a linearly divergent 
(with the lattice spacing) counterterm m^. The critical line rricr is the value of mo 
where the untwisted quark mass vanishes. We will extensively discuss in sect. 4.3 
how to deflne non-perturbatively the critical mass. For the moment we just assume 
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that such a value exists. Because of the V symmetry defined in eq. (2.39) the fiavour- 
parity violating operator xi^'^^X comes with a coefficient odd in /iq, and opposite to 
what happens to the untwisted quark mass, the twisted mass is renormalized only 
multiplicatively. To state it differently, for zero twisted mass, parity is a symmetry 
of the bare action. The residual [/v(l)3 symmetry (2.40) forbids bilinears containing 
flavour matrices r^'^, and the parity flavour symmetry P^'^ requires that parity and 
flavour are violated together, so it forbids flavour singlet parity violating terms xi^X 
and efj,upaFfj,^Fpa, as well as the flavour violating, parity even, operator x^sX- It is 
easy to check that all the dimension four operators which violate parity or isospin 
or both are ruled out by Vp'^. The flnal continuum theory (2.50) has now only an 
apparent flavour and parity breaking. As we have already explained in sect. 2.2, 
provided that all the renormalization constants are deflned in a mass independent 
scheme, this theory is equivalent to standard = 2 degenerate flavours QCD with 
a mass M = Jm\ + and no parity-flavour breaking. 



2.8 Correlation functions 



In sects. 2.1 and 2.2 we have seen that both in the classical theory, and in the 
quantum theory regularized in a chiral invariant way, there is an exact equivalence 
between QCD and tmQCD. This equivalence reflects itself in a correspondence be- 
tween correlation functions computed in the two theories. In particular eq. (2.23) 
shows the relation between correlators, valid in the bare lattice theory regularized 
with GW fermions. Based on universality arguments we expect this equivalence to 
be true also with Wtm fermions between renormalized correlation functions. 

The equations which relate correlation functions in the two theories depend on how 
the twist angle is deflned. The twist angle can be deflned in the renormalized theory, 
analogously to the continuum theory, by 

tana; = ^ = ^ ^'^ . (2.55) 

mR Zm mo - rricr 

To tune the twist angle it is then necessary to compute the ratio of renormalization 
constants ^ and the critical mass nicr- We will discuss extensively in sect. 4.3 how 
practically to determine the critical mass. One possible way to determine the critical 
mass is using the PCAC relation 

Za (goAg(x)P"(0)) 
mK=— mpcAc with mpcAc = 2^p.(^)p.(o)^ « = 1>2. (2.56) 

Measuring directly the PCAC mass, the twist angle is obtained by 

iimuj = —^ , (2.57) 

ZyimpcAc 
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where we have used the fact, imphed by the PCVC (7.2), that Zp = l/Z^. 

To tune the twist angle to = 7r/2 it is not necessary to determine any renormal- 
ization constant, but only the critical mass. 

Now that we have renormalized the theory in a mass independent scheme and we 
have determined the twist angle, using universality arguments, we can conclude that 
the relation between correlation functions in QCD and tniQCD can be inferred by 
the transformation of the integration variables (2.6), i.e. 

((9[^,^])(M,,o) = (0[x,x])(,nH„m) (2-58) 

and at finite lattice spacing is valid up to cutoff effects. The index of the correlation 
function in the l.h.s indicates that it has been computed in standard QCD with 
renormalized quark mass Mr and /xr = 0, and the index of the correlation function 
in the r.h.s indicates that it has been computed in tmQCD with renormalized quark 
masses mR and /ir satisfying 

Ml = ml + ^il. (2.59) 

Hence a given standard correlation function in QCD can be written as a linear 
combination of correlation functions computed in tmQCD at a given twist uj. 

To summarize the procedure: 

• start with the QCD correlation function you are interested in; 

• perform the axial rotation that in the continuum brings the action from the phys- 
ical basis to the twisted basis (see eq. (2.6)) on the fields appearing in the corre- 
lation function; 

• compute the resulting correlation function with the Wtm lattice action in the 
twisted basis, with a choice of quark masses; 

• perform the continuum limit. 

The final result will be exactly the desired QCD correlation function in the contin- 
uum with quark mass Mr given by eq. (2.59). 

We give here an explicit example, which is relevant for the extraction of the pseu- 
doscalar decay constants and will be analyzed again in sect. 7.1. 

The pion decay constant /t^ can be determined in the standard Wilson case from 
the correlation function 

((^r)J(x)P^(i/))(m„o). (2.60) 
With Wilson fermions the axial current is not protected by chiral symmetry, hence 
it needs to be renormalized by the scale independent renormalization constant Za-, 
which has to be determined to extract the decay constant from (2.60). If we want 
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to compute the same correlation function in tmQCD we perform first the axial 
rotation (2.6) on the fermion fields in (2.60), obtaining 

((^R)o(a^)^R(^)){MR,,o) = cos(w)((AR)J(a;)Pj^(y))(„,j,,^j,)+sin(aj)((Va)o(a;)iR(?/))(mR,^.R)- 

(2.61) 

This relation is very useful because, as we will see in sect. 7.1, there is a definition of 
the vector current which is protected from any renormalization. As a consequence at 
uj = 71 /2, the decay constant can be computed with Wtm without the computation 
of any renormalization constant. 



2.9 Exceptional configurations 



One of the historical reasons why a twisted mass term was introduced is the so called 
problem of exceptional configurations that we are briefly going to explain. 

We have seen that lattice QCD with Wilson quarks [1] breaks explicitly chiral sym- 
metry. To deal with this breaking we have to add suitable counterterms [35] in order 
to restore chiral symmetry in the continuum limit. We have seen, e.g., in sect. 2.7, 
that the ordinary untwisted quark mass is renormalized also additively. As a conse- 
quence the value of mo which corresponds to physical light quark masses is typically 
negative. This could have further practical consequences. In fact it implies that the 
Wilson-Dirac operator is not protected against zero modes. The Wilson operator 
Dw (2.29) fulfills the property Ts/^wTs = ^w- We can then define the Hermitian 
Wilson operator 

Qw = 75(^w + "^o) Qw = Qw- (2-62) 
Qw can have in general, for a given gauge configuration, a very small eigenvalue, 
even at values of rriQ which correspond to a not so small quark mass. These modes 
are expected to disappear in the continuum limit, but they can still be dangerous 
in numerical simulations 

To understand this we write a pseudoscalar density propagator 

Cp{x) = -(V;(x)75y^(x)^(0)75y^(0)) (2.63) 

in terms of eigenfunctions and eigenvalues of Qw- Performing the functional integral 
(see eq. 2.21) over the fermion fields we obtain 

Cp{x) = ^Z-' J Dp] det(Q^)tr [Q^\0,x)Q^\x,0)] e'^^^, (2.64) 

This example shows the well-known fact that a functional integral over Grassmann 
variables cannot diverge. In fact if we write the the r.h.s of eq. (2.64) in terms of 
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eigenfunctions and eigenvalues Aj of Qw we obtaiiO] 

^ 0, (0) f 0*, (0) . (2.65) 

j,k 

After integration over the quark fields, a small eigenvalue of the Wilson operator 
appears both in the fermionic determinant and in the quark propagators entering 
the correlation functions. If a small eigenvalue occurs in the course of the integration 
over the gauge fields, the contributions from the denominators are always exactly 
compensated by the same eigenvalues in the expression for the determinant. 

A problem arises however in the so-called quenched model, which consists in neglect- 
ing the fermionic determinant. The contribution of a small eigenvalue to a fermionic 
correlator is then not balanced by the determinant, leading to large fluctuations in 
some of the observables which completely compromise the ensemble average. Gauge 
field configurations where this happens are called "exceptional". The approach to 
the chiral limit in the quenched model with ordinary Wilson quarks is then limited 
by exceptional configurations. 

In refs. [36,37], to solve the problem of exceptional configurations, was suggested to 
perform a chiral rotation of the mass term. But in ref. [36] was suggested to send the 
twisted mass to zero at the end of the computation, i.e. treating the twisted mass 
as an external source not as the real mass term, and in both references the axial 
rotation was anomalous, i.e. flavour singlet. Hence it would have not been possible 
to use it as an alternativ discretization of lattice QCD. 

Adding the twisted mass term (2.28) to the Wilson action solves in a straightforward 
way the problem [15]. For the 2 flavours Wtm operator D = Dw + + ^/^qTsT^ we 
have 

Q = -i^D = Qw + ij^qT^ Q^Q = QwQw + fJ'l = Qw + l4- (2-66) 

This means that the Wtm operator does not have exceptionally small eigenvalues 
for arbitrary gauge fields: they can only appear in the massless (/iq = 0) theory. 

These considerations are true if we assume that the distribution of the eigenval- 
ues of Qw are either mass independent (like in the quenched model) or the mass 
dependence is analytic near the chiral point. We will see in sect. 6 that in infinite 
volume this is not true because of the non trivial chiral phase structure for Wilson 
fermions. On the contrary, approaching the chiral limit at fixed finite volume could 
in principle be advantageous using Wtm. 



Strictly speaking we are working in a finite volume where the spectrum of Qw is discrete. 
This is practically the case in all the numerical simulations. 



Cp{x) oc J D[U] 
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3 Non-degenerate quarks 



In this section we show how twisted mass QCD can be generahzed to a doublet of 
non-degenerate quarks. We extend the topics covered in the previous section to the 
non-degenerate case, emphasizing the main differences with the degenerate case. 

3.1 Continuum actions 

The continuum action we have discussed in sect. 2.1 describes two degenerate hght 
quarks. To add a further doublet of non-degenerate quarks in order to describe the 
heavier (c,s) quarks, two proposals have been made [38,39]. Both proposals can of 
course be used to describe a possible non- degeneracy also in the light sector. The 
first proposal [38] is based on a flavour off-diagonal splitting 

'S'f[X,X,G'] = j d^XX (Tm^M + "^q + «/^q75r^ + CqT^) X, (3.1) 

where we take /iq > and eq > 0. The off-diagonal splitting is particularly interesting 
because, as we will see, it retains all the nice properties of tmQCD at full twist and 
it keeps the quark determinant real and positive if y"T|~+7^ > Cq (see below and 
sect. 3.2). 

In order to change from the twisted basis to the physical basis the following fleld 
transformations are needed. First we need an isovector rotation of uj2 = n /2 

x' = exp(«a;2rV2)xL2=V2 = + (3.2) 

X = xexp(-z^2rV2)U2=V2 = ^^{^ - ir^)- (3-3) 

This vector transformation leaves invariant the kinetic term and transforms the mass 
terms as 

nig + ifiq-f5r^ + egT^ ^ nig - ifiq-f5r'^ + e^r^. (3.4) 
Now we perform an axial rotation as before in the direction of the twisted mass term 

^p = exp{-iuJi'y5T^/2)x, ijj = x expi-iui'y^r^ /2). (3.5) 

This transformation leaves the form of the action invariant, but transforms only the 
mass term with /iq, and, if we want to have the standard action, the rotation angle 
Ui has to satisfy eq. (2.9). The action now looks like 

Sf[iP,^] = J d^x^^-f^D^ + M + e^r')^, (3.6) 
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where M = ^Jfng + /i^ is again the polar mass. We remark here that these are 
the transformations needed given the action in eq. (3.1). Choosing a different basis 
for the action one is interested in (typically depending on the details of the lattice 
simulations), will induce different field transformations in order to connect the initial 
basis with the physical one. 



Analogously to the degenerate case we can derive the partial conservation laws 

.lab 



(9 V"- 



-2/iq e' 



3ab pb 



where 



= X75X, 



(3.8) 
(3.9) 



If we have in mind to describe with this action the heavy doublet (c,s) we will 
naturally associate the quark mass in the following way: 



rrir 



M + e, 



M 



(3.10) 



The fermion determinant will then be positive if M > eq. This constraint will be 
reconsidered when introducing the Wilson term in the lattice actions, and taking 
into account how the bare masses are renormalized. 

Another way to extend the tmQCD action to four flavours has been proposed in 
ref. [39]. In this proposal the continuum action reads 

Sf[x,X,G]= J d^xxilf,D^ + m + in-f5)x, (3.11) 

[u, d, s, c) and the mass matrices have 

\ / A/r.rns/,), n n n \ 



where now x collects four quark fields 
the form 



rn 

















rud 














nis 





V 








m 










o\ 
























I 











( Mu cos ijJi 

Md cos uJi 



\ 



( Mu sin uji 



V 








-Md sin LOi 








Ms cos oJh 

Mccosujh J 

\ 



Mc sin w/i J 



(3.12) 







Ms sin OJh 




(3.13) 



Hence the theory has six independent parameters, namely the four polar quark 
masses Mj {i = u,d,s,c), with Mf = mf + /if, and the two twist angles uJi,Uh- 
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In other words, the four standard mass parameters and the four twisted mass 
parameters fii are constrained by 



tana;/ = — = , 

rriu rrid 



tancu/i = — = . 



(3.14) 



This framework extends the degenerate two-flavour tniQCD to non-degenerate quarks 
with the property that the quark mass terms remain flavour diagonal. At vanishing 
twist angles uJi and uJh one recovers the standard QCD action of four quark flavours, 
while for ujh = ^ and = the two-flavour version of tmQCD in eq. (2.1) is 
reproduced, with two additional untwisted quark flavours s and c. 

Given the form of the mass matrices (3.12,3.13) the rotation that has to be performed 
to go back to the physical basis is 



X = exp{-iuja5Tj^/2 - iuJhl5'rh/'2)iJ, 
with 



/l 







-1 













o^ 




0/ 



X = ^exp(-iu;;75rf/2 - iuJhl^Tl/2), 



(3.15) 


















1 











and uoi and uo^ satisfying eqs. (3.14). The four flavour QCD action now takes the 
standard form 

^f[^,?,G']= Y. Jd^xi^,{x){'y^D^ + M,)^lj,{x). (3.16) 



i=u,d,s,c 



The rotation in eq. (3.15) will also give the relations among local fields in the 2 
basis. 



For a generic four flavour QCD theory with non-degenerate quarks, chiral and flavour 

symmetries are broken explicitly and this is expressed by the PC AC and PCVC 
relations (for i ^ j) 

d,,Af,^ij = {mi + mj)Pij + z(/ii + ^lj)Sij, (3.17) 

^M^M-ij = i^i - + - l^j)Pij, (3.18) 
where the bilinear fields are defined by 

Sij XiXj^ P'ij XjTsXj; AfJ.,ij Xi1^J,15Xj ^ Xil^iXj- (3.19) 
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3.2 Lattice actions 



We can now write also the fermionic action proposed in ref. [38] for two non- 
degenerate quarks in the twisted basis 



Sf[x,X,U] = a^J2x{'^ 



(3.20) 



The introduction of the off-diagonal splitting leaves intact, for example, the sym- 
metries V and Vp, while Vp is a symmetry only if combined with a sign change of 
e„ 



(3.21) 



In sect. 3.1 we have argued that the determinant is always real and positive provided 
the constraint /iq > eq is fulfillecfl. This is an important practical issue in order to 
perform dynamical simulations with the currently available algorithms, since the 
fermionic determinant is usually included with the gauge action to form an effective 
Boltzmann weight (see sect. 8 for a discussion on recent algorithmic developments). 

To understand how the condition on the bare masses translates to the quantum 
theory, we anticipate here that the renormalization factors of /iq and eq are related 
to the renormalization factors of the pseudoscalar and scalar currents, i.e. 



(3.22) 



It is then easy to show that the constraint /iq > eq augmented with the defini- 
tions (3.22) gives 

^ (/^c)r - (/is)R ^2 23) 

(/ic)R + (/is)R' 

where, having in mind phenomenological applications, we have defined 

(Atc)R = /iR + cr, (Ats)R = AiR - eR. (3.24) 

To give an example, fixing the values of the renormalized strange and charm quark 
masses, gives the following constraints 



(/ic)R ~ l.SGeV (//s)r ^ O.lGeV 



> 0.875. 



(3.25) 



The fermionic action proposed in ref. [39] for four non-degenerate flavours reads 



X] 



X{x) 



(3.26) 



For simplicity we consider here the full twist case. 
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with = (m, d, s,c), m and /j, defined in eqs. (3.12,3.13). To discuss the properties 
of the fermionic determinant with the action (3.26) we assume that the hght doublet 
is mass degenerate, 



Then the integration over the hght fermion fields yields to 



det 

Nf=2 



(Dw + mi)l + ifia^T^ 



det 



{Dw + rriiy {Dw + mi) + 



(3.27) 



(3.28) 



where the indices indicate in which fiavour space the determinant is taken [15]. 
Hence the determinant of the light doublet is positive at non-zero /i/, irrespective of 
the background gauge field. Integrating over strange and charm quarks one obtains 



det 

Nf=l 



{Dw + m^YiDw + - /i^/Xc + ifi^biDw + m^) + iHs'^^{Dw + m^ 



(3.29) 



which is real and positive only for degenerate strange and charm quarks. The other 
possibility to ensure the reality of the determinant is to employ untwisted strange 
and charm quarks, /i^ = /ic = 0, as the fermion determinants for individual Wilson 
quark fiavours are real. Even if in this case the positivity of the determinant is not 
guaranteed, recent numerical results [40] indicate that this is indeed the case. 

Practically, with the action (3.26), there are two options: do not include the strange 
and the charm quarks in the determinant, or include them but without twist. The 
first option would correspond to a partially quenched simulation with Nf = 2 light 
dynamical quarks, i.e. the strange and the charm would remain quenched. The 
second option would correspond to a Nj = 3,4 dynamical simulation with two 
light twisted quarks and two heavier non-degenerate untwisted quarks. 



3.3 Transfer matrix 



In sect. 2.6 we have seen that provided |k| < | Wtm fulfills physical positivity. In 
app. C we briefiy repeat the steps of the proof with the extension to Wtm, because 
they are important for non-degenerate quarks. 

In fact in app. C we show that if we add a non-degenerate doublet as in eq. (3.20) 
the constraint on k has to be changed into 

I'^l < a A ' > 0- (3.30) 

6 + 2aeq 

This difference can be understood by the different Hermiticity property of the twisted 
mass term ifiq'j^r^ and the splitting term e^r^. 
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First simulations with non-degenerate twisted quarks [41] indicate that suitable 
values of eq are rather small and since, in the continuum limit and for a value of 
the untwisted bare quark mass close to the critical mass, n is near 1/8, for practical 
purposes the constraint (3.30) should not cause any limitation. 

To conclude Wtm, like the pure Wilson theory, is reflection positive for all the 
relevant values of the bare parameters. 



3.4 Renormalization 



Here we shortly discuss the structure of the counterterms for the non-degenerate 
action (3.20) and we refer to the original paper [39] for the counterterm structure 
of the action (3.26). The only additional counterterm to those in eq. (2.49) allowed 
by the lattice symmetries is 

eqXT^X- (3.31) 
Because of the Vp symmetry defined in eq. (3.21), the flavour violating operator x^^X 
comes with a coefficient odd in eq, and the twisted mass splitting is renormalized 
only multiplicatively. All the dimension four operators which violate parity or isospin 
or both are absent because Vp^ are still symmetries when eq = 0. The dimension 
three operators that are absent in eqs. (2.49,3.31) are all ruled out by Vp symmetry, 
with the exception of X75T^X that is ruled out by charge conjugation C which is still 
a symmetry of the action (3.20). To conclude, the continuum renormalized quark 
action for non degenerate quarks is 

5*0 = Sg[A] + J rf^xx(x) 7^D^ + m^ + «/iR75^^ + crT"^ xi^), (3.32) 

where in addition to the degenerate case we just have to add 

eu = eqZ,{gQ,an). (3.33) 
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4 0(a) improvement 



The continuum limit of lattice QCD is of fundamental importance to relate numerical 
simulations with experimental results. Practically to perform the continuum limit it 
is crucial to simulate at several values of the lattice spacing a, and it is also possible 
(and sometimes mandatory) to improve the rate of the discretization errors from a 
to 0?, using a suitable lattice QCD action. A possibility is to apply Symanzik's 
improvement program [3-5], where the 0(a) cutoff effects in on-shell quantities 
are cancelled by adding local 0(a) counterterms to the lattice action and to the 
composite fields of interest [6-10]. A technical difficulty is that the improvement 
coefficients multiplying these counterterms are not known a priori. An alternative is 
to use Wtm at full twist. By this we mean Wtm with bare parameters tuned in order 
to have in the continuum limit a vanishing untwisted quark mass. We will show that 
in this case physical observables will be automatically 0(a) improved without the 
knowledge of any improvement coefficient. The only parameter tuning required is 
that of the bare untwisted quark mass to its critical value. 



4--1 Symanzik expansion 



The form of the unimproved lattice action is 

S[x,X.U]= Sg[U] + Sf[x,X.U] . 



(4.1) 



In this section we will analyse the 0(a) effects, so we leave unspecified the form of 
the gauge lattice action Sg, since it would only change the theory at O(a^). Sp is 
the Wtm quark action defined in eq. (2.28), that we rewrite here for convenience 



Sy[x,X,U] = a^^x( 



x\ 



(4.2) 



Following the program of Symanzik, the long distance properties of Wtm close to 
the continuum limit may be described in terms of a local effective theory with action 



5'cfT = So + aSi + 0^82 + ■ ■ ■ 



(4.3) 



The key constraint is that each term of the Lagrangian has to be invariant under 
the symmetries of the regularized theory, i.e. the lattice theory. The leading term. 
So, is the action of the target continuum theory 



So = Sg[A] + J d^xxix) 



(4.4) 



discussed already in sect. 2.7. The remaining operators 5*^ have to be interpreted as 
operator insertions in the continuum theory. The continuum theory can be defined 
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employing a lattice with spacing much smaller than a, or using a chiral invariant 
regularization that fulfills the Ginsparg- Wilson relation. The terms in the effective 
action are of the form 

S, = J d'yCkiy) (4.5) 

where the Lagrangians Ck{y) are linear combinations of local composite fields of 
dimension 4 + k. 

Cutoff effects come also from the local composite fields one is interested in. A generic 
renormalized local gauge invariant field (j)fi{x), constructed from quark and gluon 
fields on the lattice, is represented in the effective theory by an effective field of the 
form 

<Pcsix) = (j)o{x) + a0i(x) + a^(t)2{x) + . . . (4.6) 

where the fields (^k should have the appropriate dimension and should transform 
under symmetries as the lattice field. 

All on-shell quantities in QCD can be extracted from correlation functions of lo- 
cal composite fields. These correlation functions are needed at non-zero physical 
distance. We take a generic connected lattice correlation function made of a multi- 
plicatively renormalized multilocal field 

G(xi, ...,Xn) = (0r(xi) ■ ■ ■ 0ij(x„)) = ($) (4.7) 

and we always consider xi 7^ X2 7^ ■ ■ ■ 7^ a;„. In the effective theory up to order a it 
will be described by 

($) = ($0)0 -aj A($oA(y))o + a(<l'i)o + 0(a2) (4.8) 

where the expectation values on the right hand side are to be taken in the continuum 
theory with action Sq and 

('^'0)0 = (0o(a;i) 0o(a;n))o (4.9) 

n 

($1)0 = E(<^o(a:i) ■ ■■(t>iM ■ ■■Uxn))o (4.10) 

k=l 

The second term in the r.h.s. of eq. (4.8) develops potentially divergent contact 
terms whenever y = Xk- An important remark [9] is that these contact terms do not 
spoil the form of the expansion. Any contact term coming from (j)Q{x)Ci{y) when 
y X will be given by an operator that has the same dimension and symmetry 
properties as 0i(x). Since we leave the expression of 0i unspecified, the way used to 
subtract the divergence from the contact term will not change the form of eq. (4.8). 
The explicit a dependence in eq. (4.8) is not the full a dependence of the lattice 
correlation function: (pi and Ci{y) are linear combinations of fields, the coefficients 
of which depend on a logarithmically, as shown in perturbation theory [5] . Additional 
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0(a) effects can arise if one integrates the lattice correlation functions over short 
distances, and these cutoff effects will not be described by the effective theory. We 
remind here that this is not a crucial restriction since hadron masses and matrix 
elements are computed from correlation functions at non-zero distance. This remark 
is also important because it allows further simplifications in determining the set of 
operators Oi contributing to Ci. In app. D, I briefly summarize how to construct Ci 
and I give as an example the operators contributing to (pi for the currents A^, 
and P". 

The result of this analysis gives as the leading term of the effective Lagrangian 

Cl = j2c^Oi (4.11) 

1=1 

where 

Oi = tw^uF^uX, (4.12) 

O2 = mqtr{F^,F^,}, (4.13) 

O3 = mlxx, (4.14) 

C4 = mq/iqixTs^^X, (4.15) 

O5 = f^lxx- (4.16) 



Now that we know the form of the leading corrections to the effective action we can 
add to the Wtm Lagrangian the suitable counterterms in order to remove the 0(a) 
terms from the lattice action. This will be already enough to improve all the spectral 
quantities like the hadron masses. The counterterms to add to the Wtm action are 

a'EEc,a, , (4.17) 

X i=l 

where the fields Oi will be some lattice representation of the continuum Oi. In gen- 
eral the form of the lattice fields Oi is not fixed because this amounts to change the 
O(a^) terms of the theory. These discretization ambiguities allow to represent the 
gauge strength field and the local scalar density in the way they already appear in the 
Wtm action. The 0(a) counterterms O2 to O5 amount then to a reparametrisation 
of the twisted and untwisted quark masses together with the reparametrisation of 
the bare coupling qq. These reparametrisations are important if one chooses a mass 
independent renormalization scheme. We have seen in sect. 2 that it is important to 
make such a choice, in order to have equivalence between twisted mass and standard 
QCD correlation functions. We assume now that a mass-independent renormaliza- 
tion scheme has been chosen. The 0(a) improved action, also called clover action, 
is given by 

Snr^Ax, X, U] = Six, X, U] + 6S[x, X, U] , (4.18) 
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8S[x,X,U] = a^^Cs^x{x)ja-f,^F^^{x)x{x) , (4.19) 

where F^i. is a lattice representation of the gluon field strength tensor and Csw is 
the so-called Sheikholeslami-Wohlert parameter [7]. This parameter depends on the 
bare gauge coupling and all the details of the lattice action (but not on the quark 
masses) and has to be tuned in order to achieve on-shell 0(a) improvement. To 
define, consistently with 0(a) improvement and a mass independent renormalization 
scheme, the renormalized coupling and masses we have to define the theory around 
the chiral point. We have seen in sect. 2.7, in the renormalization procedure, that in 
the plane of the bare parameters the massless point is given by (mo,/iq) = (mcr,0) 
where the critical line rric is the value of tuq where the physical quark mass vanish. 
Then is natural to introduce the subtracted bare quark mass 

mq = mo — nicr (4.20) 

and to define the renormalized 0(a) improved masses and coupling constant as 

gl=~glZ,{~glafi), (4.21) 

mR = rhqZ^{gl, a/i), (4.22) 

= jx^Zf,{gl,an), (4.23) 

where /i denotes the renormalization scale dependence of the renormalization con- 
stants Z, and the improved bare parameters are given by 

~gl = gl{l + h^am^), (4.24) 

fhq = mg{l + byaaniq) + b,aafi^, (4.25) 

yiiq = /iq(l + fe^aTTiq) . (4.26) 

It easy to recognize the 4 terms which correspond in the effective theory to Oi with 
z = 2, 3, 4, 5. To summarize, to achieve full 0(a) improvement it is necessary to tune 
not only Cgw, but also the b and b parameters defined in eqs. (4.24-4.26), and the 
improvement coefficients related to the operators one is interested in. To be specific 
the renormalized 0(a) improved axial current will look like 

(Ar); = Za(1 + ^Aamq) 

where we have eliminated the operator (Os)^ (D-8) using the equations of motion. In 
principle all these improvement coefficients would need to be computed to have an 
0(a) improved evaluation of an hadronic matrix element including an axial current. 
At this point the introduction of a twisted mass has added the 6^ and bm parameters 
together with some b parameters for the improved operators. Still the number of 
improvement coefficients needed to be computed even without twisted mass is not 
negligible, especially if one considers non-degenerate quarks [42]. There are two 



a; + acA^^P'^ 



afigbAe'^'V. 



(4.27) 



31 



technical remarks to make. The set of 0(a) counterterms that have been introduced 
are shghtly redundant [29]. This generic feature of tmQCD can be traced back 
to the equivalence of correlation functions of tmQCD and standard QCD in the 
continuum limit. 0(a) improved Wtm is a one-parameter family of 0(a) improved 
theories: one of the improvement coefficient can be chosen arbitrarily. The choice is 
usually to set 6m = —1/2 because with this choice almost all the other improvement 
coefficients vanish at tree-level of perturbation theory. The second remark is to 
note that the case = is particularly interesting. In the most general case 
this corresponds to mq being an 0(a) (see eq. 4.22), and in the spirit of 0(a) 
improvement, where O(a^) effects are neglected, from eqs. (4.24-4.26) we infer that 
one remains with only one parameter and moreover mass dependent 0(a) effects 
in the bare coupling are absent. Given also the remark of the redundancy of the 
improvement coefficients Wtm at full twist is 0(a) improved just tuning the clover 
term and all the improvement coefficients needed to improve the local operators. 
The b parameters are associated with opposite parity operators and it is in principle 
possible to get rid of them in a quantum mechanical analysis (we will come back to 
this topic in sect. 6). So implementing the standard 0(a) Symanzik's improvement 
program with Wtm at full twist (m^ = 0) is cheaper in the number of improvement 
coefficients to compute, with respect to standard Wilson fermions. 

We will see in sec. 4.2 that the situation of full twist when = is of capital 
importance, and it has tremendous consequence in the cutoff effects of correlation 
functions. It is already clear then that it becomes extremely important to discuss 
the way of practically implementing on the lattice the condition = 0. 

Before discussing the consequences of working at full twist, I briefly summarize 
some numerical results obtained with Wtm adopting the Symanzik's improvement 
program just discussed. 



4.1.1 Numerical tests 

The Symanzik program requires the knowledge of a set of improvement coeffi- 
cients, that we have just discussed. A first possibility is to compute them in per- 
turbation theory. A one loop computation would leave physical observables with 
0{agQ) discretization errors and if the gauge coupling (the lattice spacing) is small 
enough it is reasonable to hope that the quantity of interest will numerically scale 
with O(a^) corrections. A further check could be to change by factors of 0(1) 
the one loop improvement coefficients and check if the scaling violations change 
dramatically or not. A better approach is to compute the improvement coeffi- 
cients non-perturbatively. Some of the improvement coefficients have been computed 
non-perturbatively, within the ordinary Wilson framework, both in the quenched 
model [43,44], in the Nf = 2 [45] and in the Nf = 3 theory [46]. 
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The improvement coefficients computed in a mass independent renormalization 
scheme do not depend on the form of the mass term and they can be used also 
with Wtm fermion. Only the additional improvement coefficients, specific to Wtm, 
need to be determined in addition. In [29] it has been shown how to implement the 
standard Symanzik program for Wtm, and many improvement coefficients have been 
computed at one loop in perturbation theory, especially the improvement coefficients 
related to Wtm. 

One is certainly interested in checking numerically if the 0(a) improvement has been 
successfully implemented and if the remaining O(a^) effects are small. In particular 
quantities which have a finite continuum limit can be computed at several values 
of the lattice spacing to check the amount of scaling violations. In particular these 
scaling studies have to be performed on a line of constant physics. With this we mean 
that changing the lattice spacing the bare parameters have to be changed keeping 
a number of (physical) quantities, corresponding to the number of bare parameters, 
fixed. 

In the quenched model it is standard to tune the quark mass keeping fixed an 
hadronic mass and to tune the gauge coupling keeping fixed the so called Sommer 
parameter [47]. tq is an intermediate distance (usually fixed to be 0.5 fm) where 
the force between two static quark is evaluated. While this quantity can be measured 
on the lattice very precisely, it has a rather uncertain phenomenological value. In the 
quenched model this is not a big problem, since the systematic error of neglecting 
the fermionic determinant is anyhow unknown, and the precise determination of the 
value of the lattice spacing allows very careful analysis of scaling violations. In the 
results presented here the values of ro/a, where needed, are taken from [48]. 

A scaling test with 0(a) improved Wtm has been performed [49] in the quenched 
model and in a finite volume (L^ x T ^ 0.75^ x 1.5) fm'^ with Schrodinger functional 



The outcome of the study is that even if some improvement coefficients are known 
only at one loop in perturbation theory, for lattice spacings a > 0.093 fm the scaling 



° These boundary conditions allow to perform in intermediate volumes scaling studies of 
the lattice theory close to the continuum limit. In particular this framework allows very 
precise determinations of several quantities which have a well defined continuum limit and 
can be used to study scaling violations. Moreover these quantities become phenomenolog- 
ically relevant in the infinite volume limit. 

^ It is interesting to note that in this study, even if the untwisted quark mass is small, we 
are not at full twist. 




0.131. 



(4.28) 
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Fig. 1. fps vs. at at a fixed value of the pseudoscalar mass Mps ~ 1.2Mj^± for non-pertur- 
batively (NP) improved Wilson fermions (□) [54] and non-perturbatively (NP) improved 
Wilson twisted mass (■) [53]. The continuum extrapolation by ref. [54] is also shown. 



behaviour of some renormalized and improved quantities (which in large volume yield 
the mass and the decay constant of pseudoscalar and vector mesons) is consistent 
with 0(a) improvement, with residual cutoff effects at a = 0.093 fm ranging from 
0.5% to 9%. 



The same setup and observables have then been employed for a study in large 
volumes [53]: L = 1.5 to 2.2 fm and T = (2 — 3)L. This study was restricted to 
two lattice resolutions, a = 0.093 and 0.068 fm, and, for each of them, three sets 
of quark mass parameters, which correspond to = 7r/2 + 0(a) and pseudoscalar 
meson masses Mps in the range 1.85 > (Mps/Mk±)^ > 0.85. 



The scaling behaviour of the pseudoscalar decay constant /ps in large volume is 
presented in fig. 1 for the 0(a) improved Wilson formulations with \uj\ ~ 7r/2 [53] 
(■) and a; = [54] (□). In the latter case four lattice resolutions were considered 
to allow for continuum extrapolation. The results for /ps Tq that are obtained from 
the two lattice formulations should agree in the continuum limit: this seems to be 
the case within the statistical errors shown in the figures. It is very interesting to 
note that Wtm at full twist, for this physical quantity, shows a much weaker a^ 
dependence compared with 0(a) improved Wilson fermions. 
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4-2 Automatic 0(a) improvement 



We have shown that 0(a) improved Wtm requires the computation of less improve- 
ment coefficients than Wilson fermions. In particular we have seen that setting Wtm 
at full twist only Csw has to be computed in order to remove the 0(a) from the action 
and from the reparametrisation of the quark masses and gauge coupling. While to 
improve operator matrix elements the standard improvement coefficients are needed. 
All the improvement coefficients depend on the details of the lattice action, in par- 
ticular they depend for example on the choice of the gauge action or on the way it 
is discretized the lattice derivative. 

In a remarkable paper of Frezzotti and Rossi [11] a step forward was made. It was 
proved that correlation functions made of parity even multiplicatively renormalizable 
fields are free from 0(a) effects, and so no improvement coefficients are needed, if 
in the continuum limit rrtR = (see eq. 4.4 and eq. 4.22), i.e. if the physical quark 
mass is given solely by the twisted mass /xr. 

We will call this property automatic 0(a) improvement. The ffist remark is that 
the theory itself is not improved, but the physical correlation functions are. To say 
it in another way, from all the possible sets of correlation functions that define the 
regularized theory, the 0(a) effects are all in those which vanish in the continuum 
limit. In particular the parity odd correlators, which vanish in the continuum limit, 
will have a ffist contribution of 0(a) while the parity even correlators will have as 
a ffist correction to the continuum value an O(a^) error. 

4.2.1 Proof 

Automatic 0(a) improvement can be demonstrated in many ways, and many proofs 
appeared in the literature [11,55-60]. The ffist proof was given in ref. [11], which 
is sketched in app. E. This proof is based on a set of spurionic symmetries of the 
lattice action. 

Automatic 0(a) improvement can be proved in a different way just considering the 
symmetries of the continuum action. In the following I try to summarize the proof 
which emphasize the role of the symmetries of the continuum action [57,58], the 
automatic 0(a) improvement of the massless Wilson operator in a finite volume [59], 
and the usage of symmetries which are not spontaneously broken in infinite volume 
continuum QCD [60]. 

Before going in details in the proof I would like to give an historical remark, in 
Minkowski space, that is not directly connected to twisted mass, but it helps to 
clarify the nature of the cutoff effects in the Wilson theory. The leading non- 
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renormalizable corrections to QED would be those interactions of dimension 5, which 
are suppressed by only one factor ofl/E, where E is some high energy scale. Accord- 
ing to Lorentz, gauge and CV invariance there is only one term allowed: the Pauli 
term proportional to 1/ Eipa^yF^^il). The contribution of this term to the magnetic 
moment of the electron or muon would give a constraint on the value of E. It is well 
known, see for instance sec. 12.3 of [14], that this constraint can be strengthened 
using the following argument. The QED Lagrangian is symmetric under the follow- 
ing chiral transformation '-y^ip and m —m, were m is the lepton mass. Then 
assuming that also the inclusion of a Pauli term should respect this symmetry we see 
immediately that the Pauli term in the Lagrangian would have to appear with an 
extra multiplicative factor m/E, i.e. m/ E'^ipa^yF^^ip. We will see now that this is 
essentially the mechanism responsible for automatic 0(a) improvement, because the 
Pauli term describes the leading discretization errors with an energy scale i? ~ 1/a. 

We first consider the massless Wtm action (eq. 2.28 with mo = mcr, yUq = 0) in 
a finite volume with suitable boundary conditions for all the fields. The choice to 
work in a finite volume is done to keep the mass dependence smooth and avoid any 
complications with possible phase transitions. 

The Symanzik effective action reads 

^eff = So + aSi + ... (4.29) 

and we are interested in a massless continuum target theory. 

5o = | d^xx{x)[i^D,]x{x) (4.30) 

The correction terms in the effective action are given by 

Sr = j d'yCM CM = Y.<^^o,iy)■ (4-31) 

i 

In the massless case the only operator contributing is 

Oi = ixa^^uF^uX, (4.32) 

which is the usual clover term. We consider now a general multiplicatively renor- 
malizable multilocal field that in the effective theory is represented by the effective 
field 

$cfT = $0 + a$i + . . . (4.33) 
A lattice correlation function of the field $ to order a is given by 

($) = ($o)o -a J dS{%C,{y))o + a{^,)o + ... (4.34) 
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where the expectation values on the r.h.s are to be taken in the continuum the- 
ory with action 5*0. The key point is that the continuum action (4.30) is chirally 
symmetric, e.g. the following discrete chiral symmetry 




«75i"^'^x(a;o,x 



X(xo,x)z75r ' 

is a symmetry of the continuum action, while all the operators in eq. (4.32), of the 
Symanzik expansion of the lattice action, are odd under the discrete chiral symmetry 
7^.5'^ of the continuum actionQ If the operator $ is a lattice representation of the 
continuum chirally even field $0; then the second term in the r.h.s. of eq. (4.34) 
vanishes. To show that also the $1 term vanishes we have to show that an operator 
of one dimension higher than the original one but with the same lattice symmetries 
has opposite chirality. To do this we introduce a symmetry that essentially counts 
the dimensions of the operators [11] 

W{—x — afi; fi), 
g3i7r/2^(_a;) (4.36) 

The gauge lattice action is invariant under I) while in the fermion lattice action 
the terms that break chiral symmetry are odd. But in particular the lattice action 

1 2 

is invariant under TZ^' x D. So the operators in $1 will necessarily have opposite 
chirality to $0- Given the fact that the continuum action is chirally symmetric 
also $1 vanishes. The conclusion is then: the chiral limit of the Wilson theory is 
automatically 0(a) improved, if we stay in a finite volume where no symmetry can 
be spontaneously broken [59]. We remark that the case of the Schrodinger functional 
is different since there the standard boundary conditions [51, 52] break the chiral 
symmetry so automatic 0(a) improvement does not appljo. 

We add now a standard mass term mRx(x)x(a^) to the action (4.30). The dimension 
5 operators contributing to Ci are now given by the operators in eqs. (4.12-4.14). 




^ Strictly speaking in the massless case the non trivial flavour structure of TZ^''^ it is not 
needed to prove automatic 0(a) improvement. It will become necessary in the massive 
case. 

* Chirally twisted boundary conditions [59] have been introduced in order to obtain a 
bulk automatic 0(a) improvement, with remaining 0(a) cutoff effects stemming solely 
from the boundaries. 
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All these operators are odd under the spurionic symmetry 

n'f = nl'^ X (mR -mR) (4.37) 

while the continuum action is even under TZl'^. This does not mean that the inser- 
tions of these operators in the Symanzik expansion vanish, but since we are in finite 
volume where the mass dependence is smooth this means that the cutoff effects are 
all of the kind O(amq). 

It seems like a similar argument could be given in infinite volume. This indeed is not 
the case because of the spontaneous breaking of chiral symmetry. Now the quark 
mass dependence around the chiral limit does not need to be smooth and even 
if all the leading 0(a) effects are odd in the quark mass because of the spurionic 
symmetry 7ll'\ a possible non-analyticity in the quark mass can generate pure 0(a) 
cutoff effects. [1^ To say it differently, the insertion of chirally odd operators does not 
vanish in the chiral limit because of the spontaneous symmetry breaking of chiral 
symmetry. 

To obtain automatic 0(a) improvement in infinite volume including a mass term, 
we have to consider as a target continuum theory for the fermion fields 



X(x), (4.38) 



where the physical mass is given by the twisted mass term. While in the massless 
case it is not possible to make a distinction between vector and axial symmetries, in 
the massive case one usually associates the symmetry broken by the mass term with 
the axial symmetry. In the following we will refer generically to chiral symmetry, the 
form of which will depend on the form of the mass term. We remind the reader that 
in the twisted basis the symmetry left unbroken by the mass term is the "twisted" 
vector symmetry (2.44), and the symmetry broken by the mass term is the "twisted" 
axial symmetry (2.44) (cf. sec. 2.4.1). 

We can repeat the same steps done in the first proof adding in Ci the term 

a = f^lxx- (4.39) 

The reason to have the physical mass term fully given by the twisted mass is that 
the continuum action (4.38) is still invariant under the discrete symmetry TZl''^, and 



^ Actually also the operators in eqs. (D.6) and (D.7), which are eliminated using the 

— 12 

equations of motion, are odd under the spurionic symmetry TZr^' . This means that in 
principle the equation of motion are not really needed to eliminate them in the context of 
automatic 0(a) improvement. 

A simple example is given by a sign(?7iq) where the sign function is still odd under 
m-q — > — m,q but it is non analytic in niq = 0. 
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if we now count the dimensions of the operators including the mass 

P = P X K ^ -/iq] (4.40) 

then TZl''^ xV is also still a symmetry of the lattice action [60]. We can then conclude 
that the second term in the r.h.s. of eq. (4.34) vanishes, and $i, being of one dimen- 
sion higher, is odd under a 7^5'^ transformation: for the same reason the third term 
in the r.h.s of eq. (4.34) vanishes. Possible contact terms coming from the second 
term amount to a redefinition of $1 as we have discussed in sec. 4.1, and so do not 
harm the proof. 

The reason why now spontaneous breaking of chiral symmetry does not spoil the 
proof, is that the twisted mass term breaks chiral symmetry in an orthogonal di- 
rection compared to the breaking of the Wilson term. Spontaneous chiral symmetry 
breaking is in the chiral flavour direction of the mass term, while the Wilson term, 
and accordingly the relevant dimension five operators Oi and are "orthogonal" , 
so are not affected by spontaneous chiral symmetry breaking. 

1 2 

It is interesting to note that TZ^' corresponds to a discrete twisted "charged" vector 
transformation defined in eq. (2.42). The actual symmetry used to prove automatic 
0(a) improvement is what would correspond in the physical basis to the standard 
charged flavour symmetry which is well known not to be spontaneously broken in 
continuum QCD [61]. 

The same proof can be repeated identically using, instead of chiral symmetry, the 
twisted parity symmetry [57, 58] 

{f/(a;o, x; 0) ^ f/(xo, — x; 0), t/(xo, x; A;) f/~^(xo, — x — a/c; /c), /c = 1, 2, 3 
x(xo,x) 7o(«75T^)x(a;o, -x) 
x(a;o,x) -> x{xo, -x) (275x^)70, 

(4.41) 

which is also not spontaneously broken in continuum QCD [62]. The only difference 
would be that for matrix elements with non-vanishing momenta, the twisted par- 
ity proof requires, an average of matrix elements computed with momentum p and 
— p, while this is not required using chiral symmetry. It is then clear that in order 
to achieve automatic 0(a) improvement, the continuum target theory must have a 
vanishing untwisted quark mass m^, otherwise the standard mass term mjiXX will 
break the residual "twisted" vector symmetry 7^5'^ (or the twisted parity symmetry 
P| ) of the continuum action. The most natural way to achieve this on the lattice is 
by setting the untwisted bare quark mass to its critical value mo = rUcr- The proof 
also shows that a possible uncertainty of 0(a) in the critical mass does not invali- 
date automatic 0(a) improvement since these uncertainties are odd under "twisted" 
vector symmetry (or twisted parity). 
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The careful reader may wonder what happens to automatic 0(a) improvement if 
the critical mass is fixed such that the untwisted quark mass is of 0(a). We recall 
that the theory itself it is not improved, but only the physical correlators are. In 
principle to obtain automatic 0(a) improvement we have seen it is necessary to have 
mR = 0, only in the continuum limit, which means at finite lattice spacing at most 
TTiq = 0(a) (see eq. (4.22,4.25)). This uncertainty can be described by a dimension 
5 operator 

Oo = A'XX, (4.42) 

where is some energy scale squared which depends on the way the critical mass 
is determined, e.g. it could be of the order of the QCD scale Aqqj-,, or it could 
be something proportional to Aqco/^q- In particular to have a massless continuum 
theory in general we have rriq = 0(a). We can say that the operator Oq parameterizes 
0(a) uncertainties in the critical mass. This uncertainty does not harm automatic 
0(a) improvement because is described by an operator which is odd under the 
transformation TZ^''^ which is symmetry of the continuum action (4.38). 

We can conclude that correlators which are even under a twisted parity or twisted 
vector transformation, are automatically 0(a) improved without the knowledge of 
any improvement coefficient, and by just tuning the critical mass such that the 
untwisted quark mass is at most rriq = 0(a). 

From the proof apparently there are no constraints on the values of the quark masses 
where automatic 0(a) improvement is at work. The presence of the Wilson term 
in the lattice action enforces us to perform a continuum limit first at a fixed value of 
the renormalized quark mass, and then to study the quark mass dependence, but it 
is important to understand how low we can go with the quark mass at fixed lattice 
spacing. 



To have a first guess we can take the polar mass 
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(4.43) 



where the rji term parameterizes the mass independent 0(a) uncertainties in the 
value of the untwisted quark mass rriq (see eqs. 4.22 and 4.25). Expanding in powers 
of a we have 



AtR 



1 + 



(4.44) 



We observe immediately that if numerically /xr < aA^, even if parametrically 0(a) 
terms are absent in (4.44), there is a term of O(a^) with a coefficient that tends 
to diverge as soon as /iR is made smaller and smaller. From this example we can 
conclude that to have an effective automatic 0(a) improvement, without big O(a^) 
effects, with a generic choice of the critical mass, such that the uncertainties in the 
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untwisted quark mass are of order aA^, we need to have the constraint 

/iR > als?. (4.45) 

From a practical point a view this constraint can be very strong. If we take the 
reasonable value A = 300 MeV and a lattice spacing a = 0.1 fm then the minimal 
quark mass that can be simulated without being affected by large O(a^) effects is 
/iR = 45 MeV corresponding, in the pseudoscalar sector, roughly to the mass of a 
kaon made up by two degenerate quarks. It is then clear that in order to go closer 
to the physical point corresponding to the up and down quark masses the constraint 
has to be weakened. From the example of the pole mass (4.44) we immediately 
understand that the crucial issue is the determination of the critical line and the 
understanding of its 0(a) uncertainties. 

4-3 The critical mass 

In this section, to let the interested reader have a general background, we give the 
basic information on the two main theoretical frameworks used in the discussion on 
the determination of the critical mass: Symanzik expansion and Wilson chiral per- 
turbation theory (W^PT). We then list the main results and discuss them, omitting 
some of the technical details of their derivations. 

The issue of the choice of the critical line was raised by the work of Aoki and Bar [55] 
and by the numerical results obtained in [63] . This problem has been further analyzed 
in several aspects in [56,57,64]. 

We have seen that to obtain automatic 0(a) improvement the untwisted quark mass 
has to be set to its critical value, i.e. to a value such that in the continuum limit 
TTiR = 0. To understand how to impose this, it is enough to understand which are the 
symmetries that are recovered in the continuum if rriR = and to impose suitable 
identities on the lattice. 

The symmetries are the twisted parity defined in eq. (4.41) and twisted vector sym- 
metry in the isospin direction 1 and 2 (2.42). One way to impose the restoration of 
twisted vector symmetry is using the PCAC relation, i.e. determining the critical 
mass setting the PCAC quark mass to zero. To restore twisted parity it is enough 
to use a twisted parity violating matrix element like a correlator between charged 
axial and pseudoscalar currents, and setting it to zero. We remark at this point 
that actually imposing the restoration of twisted parity automatically restores the 
twisted vector symmetry and vice versa. This can be understood observing that the 
discrete version of twisted vector symmetry (4.35) times twisted parity (4.41) is a 
symmetry of the lattice action, i.e. P^'^ (2.38). 
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All these options have been investigated both analytically and numerically. There 
are two possible approaches to show that the restoration of twisted parity or twisted 
vector symmetry is enough to ensure automatic 0(a) improvement down to quark 
masses satisfying a weaker constraint than eq. (4.45). One is the analysis, using the 
Symanzik expansion, of the 0(a) effects in the correlators used to define the critical 
mass. The second one is the use of a suitable modified chiral expansion in order to 
include discretization errors in the effective Lagrangian describing pion interactions. 



4-3.1 Symanzik expansion 

One possible way to restore twisted vector symmetry is to tune the bare untwisted 
quark mass rriQ to a critical value rricr such that the PCAC mass 

vanishes for large euclidean times. 

It is possible to show that the cutoff effects of the PCAC mass in infinite volume 
can be schematically written as 

r/iaAgcD + ?72a/i| + ?73aA/iR. (4.47) 

We are here implicitly assuming that the only physical scales of the theory are 
Aqcd and /ir, i.e. we have already a rough estimate of the critical mass such 
that rrtR =0(aAQCD)- Practically at this point one has several possibilities. For 
example [65, 66] for each value of it is possible to determine the critical mass 
mo = mcr(/iq) (see left panel of fig. 2) tuning the PCAC mass to zero, and then to 
extrapolate the obtained set of mcr(yUq) to /iq = (see right panel of fig. 2). This ex- 
trapolated value of TTT-cr (or equivalently Kc) can then be used to perform simulations 
for all the values of /iq. In fact from eq. (4.47) we observe that the critical mass has 
been tuned such that for all the values of fiq the PCAC mass has at most 0(a/iq). 
The slope of the curve in the right plot of fig. 2 is proportional, as it has been dis- 
cussed in [57,64,67], to 0(a) cutoff effects related to the discretization errors of the 
PCAC mass. In other words this slope is proportional to the 773 term in eq. (4.47). 
We remind that it is not surprising that the PCAC mass is not automatically 0(a) 
improved since it is an odd quantity under the twisted parity transformations (4.41) 
and discrete symmetry (4.35). 

Another possibility would be after the determination of the function m^{fiq) to use 
a different value of rricT such that mpcAc vanishes for each value of fi, 
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In practice one could smoothly interpolate the curve nicrifJ-q) at the desired value of the 
twisted mass 
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Fig. 2. Determination of the critical mass nif. (k^^ = 2amQ + 8) for a given value of /iq at 
a = 0.093 fm (left plot), and extrapolation to /iq = (right plot). The red point in the 
right plot is the critical mass determined using the "Wilson pion" definition (see sect. 4.3.2 
and 4.4 for details). 

This procedure has been used in [68] but using a slightly different correlator. One 
possible way to impose twisted parity restoration, is to tune mo to a critical value 
rricr such that the twisted parity violating correlator 

a'j:{Al{^,Xo)P\0)) (4.48) 

vanishes for large euclidean times Xq. The Symanzik expansion of the correlator (4.48) 
is, in form, identical to the one of the PCAC mass (4.47). So the same considerations 
done for the PCAC mass apply here. In particular if one performs an extrapolation 
of mcr(/iq) to /iq = 0, the cutoff effects of the critical mass would be at most of 
0(a/iq), while if is tuned to have twisted parity restoration for each value of /Xq 
then the critical mass is tuned such that mpcAc = for each lattice spacing and 
each twisted mass value. 

Another possible way to fix the critical mass, especially practical for expensive dy- 
namical simulations, is to compute the critical mass, using the PCAC relation at 
the smallest value of /iq = /imim and then use this critical mass for all the simulation 
points at heavier twisted masses. This method has been used in a recent work [69] 
where for the first time large scale dynamical simulations have been performed with 
Wtm. This is justified if at different lattice spacings the values of /xr where itIct 
is computed are properly matched. Using again eq. (4.47) one sees that for all the 
values of /iq > /imin the PCAC mass has cutoff effects at most of 0(a/iinin)- 

We immediately observe that all these methods tune the critical mass such that the 
PCAC mass is identically zero or at most an O(ayUq). 

To obtain the same result it is possible to use a clover term in the action with a 
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non-perturbatively tuned value for Cgw This will remove by definition the O(aA) 
in the PCAC mass, leaving again the PCAC mass with 0(a/iq). The critical mass 
can be determined also in the standard Wilson framework, with /iq = 0. In this 
case a non-perturbatively tuned value for Cs^, will again eliminate the O(aA) errors, 
leaving only O(amq) cutoff effects in the PCAC mass. Performing then a chiral 
extrapolation at fixed lattice spacing would implicitly determine the critical mass 
rricr up to O(a^). The only disadvantage of using this method with /iq = being 
that for quenched computations a long extrapolation to the chiral point is needed 
because of the occurrence of exceptional configurations. 

The same considerations apply if the critical mass is determined in large volume 
simulations using Schrodinger Functional boundary conditions. If the standard Schro- 
dinger Functional is used in a small volume (L < 0.5 fm), the PCAC mass properly 
improved (computed with the proper values of Csw, ca) can be used to determine 
the critical mass with residual discretization errors of O(a^) [43] without almost 
any extrapolation to the chiral limit. This is possible because the standard chirally 
breaking boundary conditions protect the spectrum of the Wilson operator from 
the appearence of very small eigenvalues, allowing simulations almost at the chiral 
point. 

What is relevant is that the cancellation of the O(aA) obtained with a properly tuned 
Csw is mass independent. This allows a tuning to full twist without a recomputation 



of the critical mass. In fact in refs. [70-72] old determination^^ [43, 73, 74] of the 
critical mass with clover improved fermions have been used. 

To summarize, all the determinations of the critical mass based on correlators which 
violate twisted parity and twisted vector symmetry, are such that the PCAC mass 
is affected at most by 0(a/iq) cutoff effects. 

If we recall the example of the polar mass (4.44) we see that now we could relax the 
constraint (4.45). It is possible to show [57] that this observation is actually true in 
general. 

In [57], it has been shown that the cutoff effects which diverge at small quark masses 
previously discussed (see eq. 4.44), so called infrared divergent (IR) cutoff effects, 
are a general property of Wtm. In general a Wtm correlator will be automatically 
0(a) improved at full twist, but could suffer from numerically large O(a^) effects 
as soon as /ir ~ aA^. The result of ref. [57] can be summarized as follows: in the 
Symanzik expansion of the lattice correlator ($) defined in eq. (4.7) at order a^'^ 
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In [70] TTT-cr has been recomputed only at one lattice spacing. 
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(A; > 1) there are terms of the kind 




2k 



k > 1 (4.49) 



proportional to the matrix element 

m\Ci\n%0))o\''. (4.50) 

It is possible to recognize the (a//iR,)^ of eq. (4.44) as the k = 1 case in eq. (4.49) 
recalling that to a first approximation oc /ir. The terms in eq. (4.49) are called 
leading infrared divergent cutoff effects, and they come from continuum correlators 
where Ci is inserted 2k times. In the following, unless specified, we use the same 
notation for the dimension 5 fields and the corresponding operators. A bit of no- 
tation is needed now: this is the continuum (see the index 0) matrix elements of 
the dimension 5 Lagrangian that appears in the Symanzik expansion of a lattice 
correlator defined by 

£i = coOo + c,Oi + (4.51) 

where Oo,i are defined in eq. (4.32) and in eq. (4.39). Actually being proportional 
to /iq the field O5 is not relevant for the following discussion. We also remind that we 
work at full twist with an unspecified estimate of the critical mass mcr.The dimension 
5 Lagrangian Ci has the quantum numbers of a neutral pion field because we recall 
we are in the twisted basis at full twist. Then it has a non-zero matrix element 
between the vacuum {Q\ and the neutral pion at rest |7r°(0)) states. 

To remove these dangerous cutoff effects, in ref. [57] it is proven that, setting the 
critical mass imposing the restoration of the twisted vector symmetry, the leading 
infrared divergent cutoff effects are removed. In particular in [57] it is suggested 
to compute for each value of /Xq the critical mass imposing that the correlator in 
eq. (4.48) vanishes for large euclidean times, and then to extrapolate to /Xq = 0. This 
leads to the result 

lim |(fi|£i|7r°(0))oP*^ = 0. (4.52) 

Incidentally in the same paper it is also suggested that an alternative possibility 
would be to use non-perturbatively improved clover fermions which would set the 
operator Oi to zero. The conclusion of ref. [57] is that with an "optimal" choice of 
the critical mass, obtained as just discussed, Wtm is automatically 0(a) improved 
for 

/iR > a^Al (4.53) 



Before ending the section we want to briefly show with an example how even without 
an extrapolation to /iq = the method of ref. [69] does not imply the existence of 
infrared divergent cutoff effects. 



45 



We set the PC AC mass to zero at a value //r = //i and as an example we consider 
again the polar mass (4.43) up to O(a^) at a value fiR = ^2 l^i 




M (/i2 + 




2 



(4.54) 



It is then clear from this example that no dangerous infrared cutoff effects appear if 
1x2 > 1^1- Even if /i2 < /ii no big enhancement are visible and if ^2 ^ /^i the cutoff 
effects take the standard form as the PCAC mass would have been set to vanish at 



All the methods we have discussed before show that all the determinations of the 
critical mass imply at most 0{afirmq) cutoff effects in the PCAC mass. We will 
see in the next section that actually from a purely theoretical point of view, all 
the determinations of the critical mass will have this property. While certainly the 
analysis of ref. [57] is correct, if we consider only 0(a) cutoff effects there is really 
no "optimal" choice of the critical mass, but they are all equivalent. We will analyze 
further this point, in particular because from a practical point of view this might 
not be the case. 



4-3.2 Wilson chiral perturbation theory 

An alternative method to analyse cutoff effects in physical observables at low energies 
is to apply the methods of chiral perturbation theory to the Wtm action, i.e. at non- 
zero lattice spacing. The discretization errors can be included systematically in a 
combined expansion in the lattice spacing and in the quark mass. 

We recall first the basic principles of chiral perturbation theory (xPT) in the contin- 
uum considering Nj = 2 flavour QCD. If the quark masses are set to zero the QCD 
action (2.1) is symmetric under the chiral group SUi^{2) x SUr{2). One then assumes 
that the theory spontaneously breaks this symmetry (see ref. [75] and refs. therein) 
down to SU\r(2) ^^\. The symmetry manifests itself in the occurrence of 2^ — 1 = 3 
pseudoscalar Goldstone bosons. In reality the QCD action contains a mass term 
which breaks the symmetry and this appears in the non-conservation of the Nother 
currents (2.18,2.19). Since the masses of the u and d quarks are small, the low energy 



Strictly speaking for Nf = 2, chiral symmetry is spontaneously broken when (m^ + 
nid) 0. In the following we will concentrate on degenerate light quarks and neglect 
physical isospin breaking effects. 



/iq = 0. 
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properties associated with chiral symmetry should show a small deviation due to the 
quark masses. A rough estimate of the size of these deviations should be given by 
the ratio between the quark masses and the intrinsic QCD scale, giving a violation 
of a few percent. 

The low energy structure of the correlation functions in QCD depends on the size of 
the quark masses. Heavy quarks play a minor role because their degrees of freedom 
are frozen at low energies. Here we consider only two flavours u and d. This is 
suitable for our purposes where we consider a lattice action for two degenerate 
flavours. In this restricted framework we are able to discuss only the dependence 
of the correlation functions on the u and d quark masses, and the remaining quark 
masses are fixed. The method we are going to briefly review [13,76] is the extension 
of the analysis carried by Weinberg [77] for the S-matrix elements, to an expansion 
of correlation functions in powers of the momenta and the quark masses. 

The method consists in adding space dependent external fields to the QCD La- 
grangian which transform accordingly in order to keep the Lagrangian invariant 
under a local chiral transformation. Then using the assumption of spontaneous 
symmetry breaking, it is possible to write a general low energy effective chiral La- 
grangian where the pion fields, i.e. the fields corresponding to the Goldstone bosons, 
are collected in a unitary matrix that automatically fulfills the requirements of chiral 
symmetry. The source terms of the QCD Lagrangian are collected in the effective La- 
grangian according to their symmetry transformations under chiral symmetry. The 
behaviour at small momenta and quark masses of the QCD correlation functions can 
be recovered matching them with the expansion of the chiral effective Lagrangian in 
powers of the derivative of the external fields and the fields themselves. We remark 
that this low-energy expansion is not a Taylor series: the pions generate poles at 
small momenta. The correlation functions admit a Taylor expansion only if the mo- 
menta are much smaller then the pion mass. So the power counting will be identified 
by the pion momentum p"^ and the pion mass (or the quark mass). In particular 
they will be treated to be both small but with the value of the ratio p'^/M^ fixed 
and unconstrained. 



The QCD Lagrangian with external fields reads 



(4.55) 



For simplicity, since not needed in the following analysis, we neglect vector and 
axial external fields and ^-terms induced by the anomaly. Co is the massless QCD 



Lagrangian that includes the gauge part)__J and the external fields s{x) and p{x 



It could in principle also include the heavier quarks. 
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are 2x2 Hermitean matrices in flavour space (we assume that p{x) is traceless) 

s{x) = s\x)l + s''(x)r^ p{x) = p"(x)r^ (4.56) 

The connected correlation functions are obtained performing functional derivatives 
with respect to the sources s{x) and p{x) on the generating functional defined as 

W[s,p] = logZ[s,p], Z[s,p] = I D[x,x]/^[?7]e-^[^'P], S[s,p] = J d'xC, 

(4.57) 

and then fixing them at their physical value: in the twisted basis the correlation 
functions for massive quarks at full twist are obtained expanding around s = and 
p = fiqT^. For example 



The local SUi^{2) x 5't/R(2) transformation of the fields is 



(4.58) 

s=0,p=/Xqr3 



X{x) ^ Vr{x)^{1 + 75)x(x) + Vl{x)^{1 - i5)x{x) (4.59) 
s{x) + ip{x) Vr{s{x) + ip{x))Vl. (4.60) 

On the other side, the pion fields are collected in a unitary matrix S which transforms 
according to the linear representation 

S(x) Vi,E{x)Vl (4.61) 

The singlet field is eliminated imposing det S = 1 , where det is applied in flavour 
space. 

The effective chiral Lagrangian will be a function of the pion fields and their deriva- 
tives together with the external fields 

/:^ = £^(S,9^S,s,p,...). (4.62) 

The order of the arguments reflects their low-energy dimensions: S counts as a field 
of order 1, (9^S as order p and s{x),p{x) as order p^. The general effective Lagrangian 
of order 1 is only a function of S and since it has to be chiral invariant it can only 
depend on det S or Tr (SS^^ where Tr is applied in flavour space, i.e. an irrelevant 
constant. We conclude that chiral symmetry implies a leading derivative coupling. 
The matrix S that collects the pion fields can be written as 



,7r"(x)r" 
7 



S(x) = Soexp(z^-^) (4.63) 
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where 7r"(x) are the pion fields, the dimensionfuU constant / is the decay constant in 
the chiral hmit normahzed to f-j^ = 93 MeV, and Eg is the vacuum expectation value 
of S, that breaks the chiral symmetry down to SUy{2). The pion fields parametrize 
the fiuctuations around Sq. 

The most general form for the effective Lagrangian consistent with chiral symmetry 
is given by 



f 



X 



(4.64) 



where the brackets (■) here indicate the trace in flavour space. The fleld col- 
lects the external fleld dependence of the effective Lagrangian according to chiral 
symmetry 

a{x) = 2Bo [s{x) + ip{x)] . (4.65) 

Analogously to what we have done for the quark Lagrangian (4.55) we can do here 
for the effective chiral Lagrangian, deflning the generating functional for connected 
correlation functions. The strategy is then to equate the correlation functions ob- 
tained in the two theories. This will express the QCD correlators as a function of the 
quark masses and the low energy constants (LEG) / and Bq. Since we are expand- 
ing around massless QCD the only scale that can appear is Aqcd, so one expects 
/ ~ -Bo ~ Aqcd- 

To include the discretization errors in effective chiral theory one proceeds in two 
steps [78]. First, one determines the continuum Symanzik action describing the in- 
teractions of quarks and gluons with momenta much smaller than vr/a. Discretization 
errors enter with explicit factors of a, and are controlled by the symmetries (or lack 
thereof) of the underlying lattice theory. Second, one uses standard techniques to 
develop a generalized chiral expansion for the Symanzik effective theory. We will 
call generically this expansion Wilson chiral perturbation theory (W^PT) having in 
mind a general form for the mass term that includes also the twisted mass case. 

The form of the Symanzik effective Lagrangian including the external sources is 
given by 

Ce{i = C + aCi + ... (4.66) 

with C given by eq. (4.55) and Ci by eq. (4.11). In order to ease the construction of 
the chiral Lagrangian in the following we keep all the terms in Ci without using the 
freedom to eliminate them through a redeflnition at 0(a) of the bare parameters. 

At this point it is useful to anticipate the power counting scheme including the 
lattice spacing a. This is 

1 > mji,fiTx,p^,a > m|,/i|,p^,a^mR/iR,mRp^amR,/iRp^a/XR,ap^ > . . . (4.67) 
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The factors of A (with A a scale of order Aqcd) necessary to make all these quantities 
dimensionless are implicit from now on unless specified. 



We recall that this approach to the description of the lattice data does not require 
a continuum extrapolation, hence the power counting scheme does not imply that 
/xr goes to zero in the continuum limit but represents only an order of magnitude 
equality. In other words the following formula are useful in describing at fixed 
lattice spacing the quark mass dependence of a given lattice correlator only in a 
region of quark masses appropriate given the power counting scheme and the value 
of the lattice spacing. 

The next step is to match the continuum effective Lagrangian (4.66) into a gen- 
eralized chiral Lagrangian. At LO the sole term that survives in Ci is the clover 
term {Oi in eq. 4.32), all the other terms being at least of NLO. The key property 
that allows this matching was noticed by Sharpe and Singleton in ref. [78] where 
they realized that the clover term transforms under chiral symmetry exactly as the 
mass term. It is the possible to add to the chiral Lagrangian (4.64) a simple source 
term in order to include the leading discretization effects. Namely the modified LO 
continuum chiral Lagrangian reads 



P 



M2) _ J 



(9^S(x)^9^S(x)) + ((T(x)S(x)t + a(x)tS(x)) + (A(x)S(x)^ + A(x)^S(x)) 

(4.68 

At the end of the analysis the sources are set to their physical value 



a{x) -> 25o(mR + ir^/in), A{x) 2Woa (4.69) 

where Wo is an unknown dimensionful constant which parametrizes the leading 
cutoff effects. We remark that if the csw coefficient would be set to its "correct" 
non-perturbative value we would have Wq = 0. This does not mean that all the 
0(a) terms will disappear, because there will be 0(a) terms at NLO, that would 
have to be cancelled by other improvement coefficients. 

Because of the chiral transformation properties of the clover term, the LO La- 
grangian (4.68) is unchanged from its continuum form if one shifts the external 
sources 

a' = a + A. (4.70) 
This at the quark level corresponds to a redefinition of the untwisted quark mass 

mn ^ + aWo / Bq = m' . (4-71) 

Recalling the definition of in eq. (4.22), this is equivalent to a shift in the critical 
mass. This observation is rather important because it means that at LO all the 0(a) 
effects in spectral quantities can be reabsorbed in the definition of the quark mass 
through the offset, or shift, aWo/Bo. This shift is not measurable, however, since 
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is not known a priori. It must be determined non-perturbatively from the simulation 
itself. The traditional definition is that m^r is the bare mass at which — > on the 
Wilson axis, i.e. /iq = 0. Since at LO the chiral Lagrangian (4.68) predicts oc m', 
we discover that this "Wilson pion" definition of m^r automatically includes the 
shift in critical mass, and chooses the untwisted quark mass to be m' . This remark 
is important because firstly it tells us that with the standard numerical definition of 
mcr, the pion and vacuum sectors are automatically 0(a) improved at LO in W^PT 
for any twist angle, and secondly it indicates that from a theoretical point of view, 
at least at LO, there is no difference in the critical line computed in this way or by 
using other methods which set m' = 0. We will come back later to this point when 
analysing numerical results. 

Analogously to what we have done in the continuum, we define a polar mass and a 
twist angle that includes this 0(a) shift 

MV"'"' = (m' + z/iRr=^). (4.72) 

In the continuum the PCAC quark mass (2.18) is the untwisted quark mass which 
appears in the Lagrangian. If we consider now the lattice PCAC quark mass defined 
in (4.46), evaluating the correlators at large distances in order to let the single pion 
state dominate one finds 

"^PCAC = m'. (4.73) 

This shows that this quantity automatically includes the 0(a) offset in the untwisted 
mass, exactly as the pion definition does. We remark that at this order it does not 
matter if the critical line is computed on the Wilson axis or with fi^y^O. 

At NLO we should add to the Symanzik expansion the O(a^) terms (see eq. 4.67). 
These can be collected in three categories. First, those that are invariant under 
Euclidean and chiral symmetries, and they simply modify, in the chiral Lagrangian, 
the leading order continuum results by a^, i.e. they lead to an O(a^) correction to the 
LEC / [79]. Second, there are four-fermion operators which violate chiral symmetry. 
In the chiral Lagrangian the corresponding operator is already present, having been 
produced by two insertions of the clover term [80]. This shows that what is relevant 
for matching are the symmetries broken by the operators (here, chiral symmetry), 
and not their detailed form. The four-fermion operators simply change the unknown 
low energy constant corresponding to a double insertion of the clover term. This 
also means that using a non-perturbatively 0{a) improved quark action, this low 
energy constant does not vanish. Finally, there are the terms violating Euclidean 
symmetry. These can be decomposed into Euclidean singlet and non-singlet parts, 
and can be shown to be of higher order [79] . 

The detailed form of the Wilson chiral Lagrangian at NLO can be found in many 
papers and reviews (see for example [81] and refs. therein). Here we want to list the 
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main results [55,56,82,83] and add some considerations. 



The charged pion mass squared at NLO reads 

+2aBoM' cos uj'{2Sw — '^vj/)+2a^w^' cos^ co' . 

(4.74) 

The first term in the r.h.s. is the continuum NLO resuh [76] while remaining terms 
show the impact of a finite lattice spacing. In particular 6w and 6^ parametrize NLO 
0(a) effects of order Aqcd- The w' term parametrizes NLO O(a^) effects of order 
Aqqj-,. We have reabsorbed the scale dependence of the chiral logs in the LEG A^r, 
and 26w is scale invariant. Chiral logs do not contain discretization corrections 
at this order because the LO discretization errors can be absorbed into a'. 



MK = 2BoM' 



'''°^'log(25oM7A- 



32n^P 



The result (4.74) shows the different possibilities for removing 0(a) errors. 

• Non-perturbatively 0{a) improve the quark action, in which case 26w — = 
and the 0(a) term vanishes. 

• Use "mass averaging" [11] in which one averages over u' and cj' + vr at fixed M'. 
This flips the sign of both m' and /xr, and thus of cosu', and cancels the 0(a) 
term. 

• Work at full twist, u' = ±7r/2. This removes the 0(a) term and, in this case though 
not in general, also the O(a^) term. This can be traced back to the exact "charged" 
twisted chiral symmetry (2.41) of the massless Wtm lattice action (2.28). 

There are two further features of the result (4.74). The 0(a) errors are determined 
for every twist angle by the combination 26w — <^^- In particular, on the Wilson 
axis, this term predicts an asymmetry in the slopes on the two sides of m^. The 
asymmetry is defined by [64] 

J,) '.,\) ry with /.R fixed. 4.75 

Neglecting chiral logs and O(a^) terms and setting /xr = in eq. (4.74) the asym- 
metry is given by 

"^'y {25w - 5^) . (4.76) 



This asymmetry has been previously observed numerically [84], in the course of 
the initial studies of the properties of Wtm. I show an example of the results in 
Fig. 3. The observed asymmetry in fig. 3 can be viewed by the difference between 
the absolute values of the slopes of the two straight lines that leads to a value for 
the asymmetry ~ 0.3 [64]. This estimate is consistent with the expected size of 
{25w — (5^) ~ Aqcd — 0.3 GeV considering that ^ 1 GeV. 
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Fig. 3. Unquenched results for (aMpg)^ as a function of k ^ 
with a 0.2 fm [84]. Straight hues are to guide the eye. 

4.3.3 Full twist 



2niQ + 8 for //q = and 



The form of the 0{a) correction in the result for in eq. (4.74), is the generic 
structure of the 0(a) terms in all the physical quantities, i.e. all the 0(a) cutoff 
effects are proportional to acosu;'. This is the way to observe automatic 0(a) im- 
provement in the framework of W^PT. To obtain automatic 0{a) improvement one 
needs cos lj' = 0(a) and thus uj' = 7r/2 + 0(a). In other words, full twist means 
at most "up to 0(a)". I want to briefly explain here in which sense the analysis 
carried out with the Symanzik expansion and W^PT gives a consistent picture for 
the determination of the full twist setup. 

In sect. 4.2 and 4.3 we have already discussed in the framework of the Symanzik 
expansion that to achieve automatic 0(a) improvement down to masses obeying the 
constraint (4.53) the untwisted quark mass has to be tuned to be zero up to errors 
at most of O(ayUq). In the power counting scheme (4.67) we assume for the W^PT 
analysis, with /xr ~ a, an 0(a) accuracy in u' requires m' = 0{a^). This is easily 
obtained considering that 

777' 

coscu' = — . (4.77) 
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If we again recall our working power counting it is easy to see that m' = O(a^) 
corresponds exactly to the requirement of having a vanishing PC AC quark mass 
up to 0(a/iq) or O(a^). So the two analyses are consistent and the result can be 
summarized as follows: to achieve automatic 0(a) improvement down to a twisted 
mass which obeys the constraint (4.53) we have to tune the theory to full twist, 
meaning that the critical mass has to be tuned such that the PCAC quark mass 
vanishes up to O(ayUq) or O(a^). 

The traditional definition of used with Wilson fermions is to extrapolate mo 
to the point where = 0. We now consider the NLO expression for the squared 
charged pion mass (4.74) and set /iq = 0. Eq. (4.74) tells us that if one is able to 
perform a perfect chiral extrapolation including the chiral log, this method is in 
principle adequate to tune m^ such that m' = O(a^) in the regime /xr ~ a. 

Alternatively to obtain m^ = in the continuum limit, it is enough to understand 
which are the symmetries that are recovered in the continuum if m^ = and to 
impose suitable identities on the lattice. For m^ ^ twisted parity (4.41) and 
twisted flavour (4.35) are broken symmetries. Enforcing this restoration in particular 
correlators for a 7^ gives a non-perturbative determination of the twist angle (or 
equivalently of mcr) to full twist. 

A way to understand how to impose the restoration of twisted parity symmetry in 
the continuum limit is to write a correlator which breaks parity in the physical basis. 
We have seen in sec. 2 which is the relation between fields in the twisted and physical 
basis. The idea is now to take either (2.12,2.13) or (2.14,2.15) as a definition of u, 
and enforce parity restoration in a particular correlator. Two examples are the lja 
method [84] 

{Vl{x)V\y)) oc (0|V>^) = 0, (4.78) 

and the up method [56] 

{S\x)Aliy))^{0\S'\n')=0. (4.79) 

The correlators are to be evaluated for x ^ y and large euclidean times in order 
to isolate the single pion contribution. Using (2.12-2.15) one can manipulate these 
criteria into results for the twist angle in terms of correlators in the twisted basis: 

{V',{x)P\y)) , ^{S\x)Al{y)) 
tantu^ = — -^7— - — ^ , , tantup = — — , , , , , . (4.80) 
{Al{x)P\y))' ^ 2{P\x)Al{y)) ^ ' 

Full twist occurs when the denominators vanish, i.e 

= 7r/2 ^ {Al{x)P\y)) = , = 7r/2 ^ {P\x)Al{y)) = . (4.81) 

The correlator in the "cup method" includes quark-disconnected contractions and is 
much more difficult to calculate in practice. The "cjyi method" is used in practice. 
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co^ method 




PCAC method 




Wilson pion 



Fig. 4. Sketch of the different methods for working at full twist. The arrows represent 
the direction one moves to approach the chiral limit (which occurs at /Uq = 0). The plot 
represents the regime in which Mr ~ a. The regime in which Mr ~ o? is represented by 
the shaded region. It is discussed in sec. 6. For the "Wilson pion" method the supposed 
behaviour given the quenched numerical data is plotted. See text for discussion. 

One fixes /Xq and varies mo until oja = vr/2. The resulting mcr(/iq) depends on the 
choice of discretization of the axial current (e.g. 0(a) improved or not), and, in 
general, upon the separation x — y. In fig. 2 a typical /iq dependence of the critical 
line is shown. At large distances, which are used in practice, the pion contribution 
dominates and the resulting becomes independent of separation. At such dis- 
tances the "cua method" is equivalent to the vanishing of the PCAC mass (4.46). We 
remind that the determination of uj at full twist does not require any computation 
of Z-factors [15]. 

Both methods can be studied in W^PT and at full twist we have [56]: 



To summarize: all the known methods to compute the critical line are theoretically 
equivalent in order to achieve automatic 0(a) improvement provided /ir > o^Kqqj^. 
Practically the situation can be very different as we are going to discuss in the next 
section. 

4-3.4 The bending phenomenon: a closed chapter 

The discussion up to this point makes clear the importance of accurate tuning to 
full twist. Initial studies of Wtm in the quenched model observed a phenomenon 
called "bending" . Although this is largely a closed chapter in the history of Wtm, it 
is worth learning the appropriate lessons. Here we assume that we work in a region 
of quark masses such that /ir ^ a^. We will enter into the region where /ir ~ a^ in 
sect. 6. 



o;^ = 7r/2 ^ uj' = Tx /2 + a5w , cup = 7r/2 => uj' = tx/2, 



(4.82) 
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I show in fig. 5, the numerical results for a/ps at a fixed value of the lattice spacing 
(a = 0.093 GeV) using several definitions for the critical mass (see the caption 
for more details). A clear "bending" is observed at smaller quark masses with the 
"Wilson pion" definition. This bending can be explained taking the LO expression 
in WxPT for the pseudoscalar decay constant 



where / is the pion decay constant in the chiral limit and /iR and mpcAc represent 
the renormalized twisted and PCAC quark masses respectively. We have already 
seen that from a theoretical point of view, independently on the definition used for 

tuning mcr, "n^pcAc vanishes or is a quantity of 0(a/iq) or O(a^). The term ( "^pcac j 



causes a deviation of the expected chiral behaviour as soon as /xr ~ '^pcac- The 
fact that for the "Wilson pion" definition a "bending" appears for /iq ~ aAQQjy, 
gives an indirect evidence that mpcAc = 0(a). In practice in a quenched simulation, 
a long extrapolation along the Wilson axis is needed because of the occurrence of 
exceptional configurations (see sec. 2.9). It is conceivable that performing this long 
extrapolation without including the chiral logs and in a region of quark masses where 
the applicability of NLO W^PT is debatable, the critical mass could be determined 
in such a way that mpcAc is of order a. If this happens, this could generate, as we 
have seen, large O(a^) corrections. A possible way to resolve the issue would be to 
perform chiral fits using NLO W^PT, including mpcAc as a fit parameter for several 
lattice spacings and then studying the lattice spacing dependence of mpcAc- 

An alternative approach is to use the clover improved Wtm action (4.18) which 
removes from the theory all the cutoff effects of the kind qAqq^) [57], and to use the 
equivalent pion and PCAC determination of rricv This is nicely shown in fig. 5 where 
both the "Wilson- clover pion" definition and the "Wilson- clover PCAC" definition 
give no evidence of bending down to masses yUq ^ o^Aq^^) [^'^] ■ 

The final answer is given by a study of the lattice spacing dependence of several 
physical quantities, with several definitions of m^r. 

The "Wilson pion" method was used in the first scaling study of Wtm in the 
quenched model [85]. This study was performed in a region of quark masses where 
most probably /xr ^ a. Therefore the "Wilson pion" definition was good enough 
to obtain a clear evidence of automatic 0(a) improvement, and no enhanced O(a^) 
effects. 

In a further publication [66] the PCAC method and the "Wilson pion" method 
have been compared, obtaining consistent results in the continuum limit but rather 
different O(a^) effects. This is clearly seen in the fig. 6 where JpsTq is plotted as a 
function of {a/r^Y pseudoscalar masses in the range 297-1032 MeV, obtained 




(4.83) 
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Fig. 5. Chiral behaviour of the pseudoscalar decay constant at fixed lattice spacing 
a = 0.093 fm for several determinations of the critical mass nicr- "Wilson-clover PCAC" (■) 
and "Wilson-clover pion" (□) [71], "Wilson pion" (o) and PCAC (•) [65,66], uja (a) [68]. 

with both definitions of m^. For small enough values of the lattice spacing, the 
values of /ps'^o show, with both definitions of rricr, a linear behaviour in {a/ro)^. 
This nicely demonstrates the 0(a) improvement for both definitions of the critical 
mass. However, for the pion definition we notice that the effects of O(a^) are rather 
large, in particular at small pseudoscalar meson masses of 297 MeV and 377 MeV. 
In contrast, the PCAC definition reveals an almost fiat behaviour as a function 
of (a/ro)^ even at these small pseudoscalar meson masses. In order to take the 
continuum limit, one should identify the scaling region where the data are well 
described by corrections linear in (a/ro)^. This turns out to start at a = 0.093 fm 
with the pion definition of mcr and at a = 0.123 fm with the PCAC definition. The 
values in the continuum limit are obtained separately by performing linear fits to 
the data in these two regions. The results of these fits are shown in fig. 6 together 
with the simulation data. It is very reassuring that these independent linear fits lead 
to completely consistent continuum values. 



4.4 The critical mass: summary 



All this discussion about the critical mass deserves a summary. In fig. 4 is sketched 
the approach to the chiral limit according to the different definitions used for the 
critical mass mcr. For the "Wilson pion" definition the supposed behaviour is shown 
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Fig. 6. ro/ps as a function of (a/ro)^ using the "Wilson pion" definition (open symbols) 
and the PCAC definition (filled symbols) of the critical mass; fits are performed with a 
linear function in (a/ro)^ separately for each set. 

given by the quenched numerical data. While they are all nominally at full twist: 
they could show different O(a^) scaling behaviour depending on the value of the 
quark mass /iR. In table (4.4) I collect all the definitions discussed and here are 
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some considerations: 

• "Wilson pion" : this is a vaUd definition to guarantee automatic 0(a) improve- 
ment down to /^R ~ a. If this method is used in the quenched case a long chiral 
extrapolation is needed because of the occurrence of exceptional configurations. 
The correct extrapolation to the chiral limit, including eventual chiral logs, is cru- 
cial in order to have a definition of the critical mass that absorbes the 0(a) offset. 
If the extrapolation is not done correctly this definition can miss the critical line 
in such a way that m' = 0(a). This is most probably the reason of the "bending 
phenomenon" . 

• "Wilson- clover pion" : this definition has the advantage with respect to the pre- 
vious one that, if one removes non-perturbatively the 0(a) effects in the lattice 
action, even if the extrapolation is needed it can miss the critical line only by 
O(a^). The nice results of ref. [71] indeed seem to indicate that this definition is 
good enough to avoid the "bending phenomenon". This somehow confirms that 
the origin of the "bending phenomenon" is the 0(a) effects present in the un- 
twisted quark mass. 

• "Wilson-clover PCAC" : to this definition the same considerations apply as for 
the "Wilson-clover pion" method. 

• ua_ '■ this definition (and the equivalent using the PCAC mass), allows a clean 
determination of the critical line mcr(/iq) for each value of fiq. It can be safely used 
in the quenched case, because the twisted mass gives a sharp infrared cutoff to the 
lattice theory. It has the drawback that the unquenched case requires a tuning for 
each value of /Xq and this could be rather expensive from the computational side. 

• LJp : this definition, theoretically very attractive, has the drawback that it requires 
the computation of quark disconnected diagrams. These diagrams are usually more 
difficult to compute so it has mainly an academic interest. 

• PCAC : this definition takes the one obtained with the ua method (or the equiv- 
alent with the PCAC mass) and then extrapolates mcr(/iq) to /Xq = 0. The advan- 
tage of this definition is that one could perform the extrapolation using a limited 
range of values of /Xq, and then the critical mass mass could be used for all the 
values of /iq one is interested in. A possible cheaper, from a computational point 
of view, alternative has been recently proposed [69] that uses the critical mass 
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determined at the lowest value of /iq used in the simulations. 

We can conclude this lengthy and technical discussion on the critical mass with the 
following statement: provided we are away from the region /iR ~ a^Aq^j-, all the 
definitions given above are theoretically equivalent in order to achieve automatic 
0(a) improvement. The "Wilson pion" definition has to be avoided in the quenched 
case because the long extrapolation to the chiral limit can induce effective 0(a) 
cutoff effects in the untwisted quark mass. 



4-5 Numerical results 



The first numerical evidence of automatic 0(a) improvement was given in ref. [85], 
where a first scaling test was performed at a fixed value of the pseudoscalar meson 
mass around Mps — 700 MeV in the quenched model. In this paper, scaling vio- 
lations of the pseudoscalar decay constant and the vector meson mass have been 
studied, confirming that, without any improvement coefficient, Wtm at full twist 
is consistent with automatic 0(a) improvement, with first indications that the re- 
maining O(a^) effects should be small. This very interesting result has triggered a 
set of quenched studies [63, 65, 66, 68, 86, 87] to further check the property of auto- 
matic 0(a) improvement [11], and to gain experience with this formulation of lattice 
QCD. An interesting quantity to compute with Wtm is the pseudoscalar decay con- 
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Fig. 7. Left plot: continuum extrapolation of /ps as a function of a^. Right plot: continuum 
limit values for fps as a function of Mpg in physical units. The empty squares are taken 
from [54]. 

stant /ps. As it was noted in [53,85,88], the computation of /ps does not require 
any renormalization constant, in contrast of ordinary Wilson fermions, and more- 
over, given automatic 0(a) improvement, does not need the computation of any 
improvement coefficient. Thus the situation for this quantity is like with Ginsparg- 
Wilson fermions. The reason for this nice property is that the continuum Ward 
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identity (2.19), remains exact on the lattice 

{d;V^{x)Om = -2^i^e"^\P\x)0m a = 1, 2, (4.84) 

(where (9* is the lattice backward derivative, and O is a local lattice operator) if one 
uses a slightly modified point-split vector current 

-If 

yil{x) = 2 - l)y (a^' ^^)x{x + afi) 

+ X{x + a/i)(7/. + ^)\u{x, ^lY\{x)Y (4.85) 
Given the fact that at full twist we have 

((^R)J(a;)pi(0))(M,,o) = ((V^R)o(^)^'(0))(n^H.m), (4-86) 

inserting a complete set of states in a standard fashion on the r.h.s. of eq. (4.86) 
and using the PCVC relation (4.84) we obtain 

/ps = |^|(0|PlP5)| a = 1,2 (4.87) 

where Mps is the charged pseudoscalar mass and \PS) denotes the corresponding 
pseudoscalar state. 

In fig. 7 (left plot) the continuum limit of ro/ps, the critical mass being computed 
with the PCAC method, is shown, from ref. [66] as a function of (a/ro)^. The scaling 
is consistent with being of O(a^), and moreover it is reassuring that the O(a^) effects 
are rather small for all the pseudoscalar masses investigated down to Mps = 272 
MeV. The right panel of fig. 7 shows the chiral behaviour of the continuum pseu- 
doscalar decay constant [66], compared with the non-perturbatively 0(a) improved 
data of [54]. We remark that this comparison is purely illustrative since it is in 
the quenched model, and the simulations with clover fermions had to stop around 
Mps — 500 MeV due to the appearance of exceptional configuration. Using a linear 
extrapolation of these data to the chiral limit (performed on the six smallest masses) 
gives the values of the pion and kaon decay constant /,r and fx (the latter in the 
SU{?)) symmetric limit). The ratio of the two gives fx/ f-w = 1-11(4), which is 10% 
smaller than that obtained experimentally. This is however consistent with what 
was observed in previous quenched calculations [89]. 

Another interesting quantity that has been computed is the vector meson mass. 
Applying the standard axial rotation (2.6) we obtain for = f 

((VR)i'^(x)(VR)i'^(0))(M„o) = ((^r)J'^(x)(Ar)J1(0))(^,,,,). (4.88) 
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Another possible interpolating field for the vector channel is given by the following 
component of the tensor current 

V = ^aoJ-^i^. (4.89) 

This current is invariant under the rotation (2.6), so the vector meson mass can be 
extracted from the correlators 

3 3 

Cl{x,) = \Y.Y.{K{^)Km « = 1>2 (4.90) 

k=l X 

CT(^o) = ^EE(^."(^raO)) « = 1,2 (4.91) 

k = l X 

One observation [66] is that the tensor correlator systematically shows smaller sta- 
tistical fluctuations. 
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Fig. 8. Left plot: continuum extrapolation of Myro as a function of (a/ro)^. Right plot: 
continuum limit values of the vector meson mass (filled circles) as a function of the pseu- 
doscalar meson mass squared. The continuum limit vector masses obtained with non-per- 
turbatively improved Wilson fermions (open squares) [54] are also plotted. 

In the left plot of fig. 8 we show the results for the vector meson mass as a function of 
(a/ro)^. Again, even for small pseudoscalar meson masses, the behaviour of the vector 
meson mass is almost flat in (a/ro)^, indicating that O(a^) lattice artefacts are also 
small for this quantity. The lines in fig. 8 (left panel) represent linear fits of these data 
as a function of (a/ro)^. The continuum extrapolated values for the vector meson 
mass are presented in the right plot of fig. 8. As a function of the pseudoscalar meson 
mass squared, they show a linear behaviour. In fig. 8 we also plot the continuum 
values obtained with non-perturbatively improved Wilson fermions [54]. The results 
of a linear extrapolation (performed with the seven smallest masses) are used to 
compute the values of Mp and of Mk* (the latter in the SU{?)) symmetric limit). As 
already observed in quenched calculations where the scale is determined through ro, 
these values turn out to be 10 — 15% larger than the experimental values. 
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4-6 Conclusions 



Despite the fact that Wtm is a relatively new lattice action, extended scaling tests 
with mesonic quantities in large volumes have shown beyond any reasonable doubt 
that automatic 0(a) improvement is at work, i.e. at the price of tuning one param- 
eter, namely the critical mass rricr, the physical correlation functions exhibits only 
O(a^) scaling violations. The issue of the choice of the critical mass can be summa- 
rized as follows: provided the definition of the critical mass is theoretically and prac- 
tically valid each definition leads to automatic 0(a) improvement. This statement 
obviously does not include regions of the bare parameters where phase transitions 
can occur. In these regions the validity of the Symanzik expansion itself can be ques- 
tioned. The importance of using both the Symanzik expansion and its application 
with xPT theory has been shown in all its power. Obviously the remaining O(a^) 
contributions of the theory are not predictable, and extensive non-perturbative sim- 
ulations have to be performed in order to understand the amount of the remaining 
scaling violations. In the quenched model the remaining O(a^) scaling violations 
are very small for a wide range of pseudoscalar masses (SOOMeV < Mps < IGeV). 
Preliminary results [90, 91] indicate that this is true also for Nf = 2 dynamical 
simulations. 
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5 The physical basis 



In the previous two sections we have seen how renormalizability and 0(a) cutoff 
effects of Wtm can be analyzed using the twisted basis. In the twisted basis the 
Wilson term takes the standard form while the mass term takes a "twisted" form. 
In this section we want to rederive some of the results already presented using the 
so called physical basis [11]. The physical basis is obtained from the twisted basis 
performing an axial rotation in the isospin directionin in the bare lattice theory. 
If one rotates the bare fields in the action and in the correlation functions, the 
physics is completely unchanged and correlation functions calculated before and 
after the rotation are the same on every single gauge configuration. This happens 
because what matters are the relative angles between the flavour-Dirac directions of 
the mass term and the Wilson term. These angles are not changed by a bare field 
rotation. 

In the physical basis the Wilson term is chirally rotated while the mass term takes 
the standard form. The advantage of this basis compared to the twisted one is that 
the dictionary of the correlation functions is unchanged since after performing the 
continuum limit, the continuum QCD action takes the standard form. On the other 
side the dictionary has to be changed from the standard one, in the renormalization 
process, because the Wilson term is now chirally rotated. 



5.1 Chirally rotated Wilson term 



The setup is identical to the one described in sect. 2.3. The only difference is the 
basis in which the fermion action is written, namely in this section the physical 
basis given by the set of fermion fields {ip,ip}. For a SU{2) flavour doublet of mass 
degenerate quarks in the physical basis the action has the form 



SF[tP,^lj,U]=a^Y.^{ 



X] 



(5.1) 



where 



Aw(^) = + V;) - are-^-^^-V;V,}, 



(5.2) 



is the twisted Wilson (tW) operator, V^, V* are the standard gauge covariant 
forward and backward covariant derivatives defined in app. A and M is the standard 
mass term. 



To ease the discussion in this section we will concentrate on the full twist case 
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uj = tt/2 and we set r = l.Then eq. (5.1) becomes 



X 



Aw + M 



(5.3) 



with 



Aw = Aw(^ = vr/2) = ^ 7^(V^ + V;) + az75r=^V;V 



(5.4) 



The choices uj = 7r/2 (and equivalently uj = — 7r/2) is particularly useful because as 
we have seen in sect. 4 it can be shown [11] that, despite the fact that the theory is 
not fully 0(a) improved, cancellation of 0(a) effects in quantities of physical interest 
(like energies and operator matrix elements) is automatic. We will prove again this 
property in this section using the physical basis. 

We shortly recall now the symmetries of the lattice action given in eq. (5.3): gauge 
invariance, lattice rotations and translations, charge conjugation C (see app. B for the 
definition), and the U{1) transformations associated with fermion number. It differs 
from standard Wilson fermions in two important ways. First, the flavour SUy{2) 
symmetry is broken explicitly by the Wilson term down to the f/v(l)3 subgroup 
with diagonal generator 

f/v(l)3: |^(-)--P(^^-^)^(-)' (5.5) 
\ip{x) ip{x) exp{—i^T'^). 

In sect. 2.4.1 we have seen that Wtm preserves a subgroup of chiral symmetry. In the 
physical basis this property is even more transparent since the lattice action (5.3) 
with zero mass (M = 0) is invariant under the standard axial transformations 



^Ij{x) ^ exp{i^-f5T^''^)ip{x), 
ip{x) ip{x) exp(z^75r"' 



f/A(l)i,2: - (5-6) 



The massless theory (eq. 5.3 with M = 0) is invariant under the group transforma- 
tions 

f/A(l)l®f/A(l)2®t/v(l)3. (5.7) 

The "charged" sector has a continuum like behaviour and the exact symmetry (5.7) 
protects the charged pion from chirally breaking cutoff effects. 

The discrete version of this symmetry 

^1,^ te) - .75r^'^^(x) ^^ g^ 

I '?/^(x) — > iil){x)'^^T ' 

is still a symmetry of the massive lattice action (5.4) if combined with a sign change 
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of the mass 

nl'^ = nl'^ X [M ~M], (5.9) 
analogously to what happens in continuum QCD. 

To understand the structure of the counterterms, we can use the symmetries of the 
lattice action (5.3) as we have done in the twisted basis. The counterterms to the 
action with dimension less or equal four are 

tr{F^.F^.}, i^75rV, M^^. (5.10) 

We notice immediately that quantum corrections imply the necessity to add to the 
lattice action a counterterm iip'j^r^ip with a linearly divergent coefficient. This is 
not surprising since this is the term that corresponds to the critical mass introduced 
in sect. 2, which is the linear divergence induced by the presence of the Wilson term. 

On the contrary because of the TZl'"^ symmetry (5.9), the operator tpip comes with 
a coefficient odd in M, and the mass is renormalized only multiplicatively. To state 
it differently, for zero quark mass the theory still preserves a remnant of chiral 
symmetry. The residual t/v(l)3 symmetry (5.5) forbids bilinears containing flavour 
matrices r^'^. The parity flavour symmetry Vp'^ defined in eq. (2.38) takes the same 
form in the two bases and thus is still a symmetry of the lattice action (5.3). This 
symmetry requires that parity and flavour are violated together so it forbids flavour 
singlet parity violating terms ip'j^ip and e^yp„F^yFp„, as well as the flavour violating, 
parity even, operator ipT^ip. It is easy to check that all the dimension four operators 
which violate parity or isospin or both are ruled out by Vp^. The lattice theory 
has to be modifled in order to include the power divergent subtraction, and to keep 
consistency between the basis we write the lattice action 

SF[^P,^, U] = a'Y.^{x) 1^ [7^(V^ + V^) + tj.r^aVy, - m,,)] + m} ^(x). 

(5.11) 

One can easily check now that this action is completely equivalent from the action 
in the twisted basis (2.28), just performing the rotations (2.6) with u = tt/2 and 
identifying mo = rricr and M = fi^- The flnal continuum theory has now the standard 
QCD form 



So = d'^x'ijj{x) 



^^Dp + Mr 



ij{x), with Mn = ZM{glaii)M (5.12) 



To understand how the choice of a twisted Wilson term influences the Ward identities 
of the theory we rewrite some of them in the physical basis. To deduce them it is 
enough to start from the Ward identities in the twisted basis, where using a mass 
independent renormalization scheme, one can make a standard analysis to construct 
the lattice flelds that are multiplicatively renormalizable and respect, up to cutoff 
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effects, the chiral multiplet structure. The final step is then to rotate all the quark 
fields back into the physical basis. 



Renormalized vector and axial currents can be taken to be 



ZyAl 



1,2, 



(5.13) 
(5.14) 

where the local currents in the physical basis are defined in eqs. (2.12,2.13) and the 
finite renormalization constants, Zy and Za, are those for the local vector and axial 
currents of standard Wilson fermions, respectively. Notice the switch between Zy 
and Za for the currents with flavour a = 1,2, due to the presence of the factor 75X3 
in front of the Wilson term in eq. (5.4). 



The expressions for the renormalized densities are 

3 ipp{aM) 



V^ = ZpV'' Vi = Zso 



5° 



Mpso{aM) 



(5.15) 



(5.16) 



where pp{aM) is a polynomial in aM, while pso{aM) is a polynomial with even 
powers of aM. The reason for this is because TZ^ is a symmetry of the lattice 
action (5.11). Formula (5.16) is rather interesting because it shows that the chiral 
order parameter is only affected by an M/a? power divergence, analogously to what 
happens with Ginsparg- Wilson fermions. 

The non-singlet Ward identities (2.20), once the local operators are properly renor- 
malized, are vahd at fixed lattice spacing up to cutoff effects 



](x)On{y)) = 2M^{V^{x)On{y)) 



with O being a generic multilocal operator and x ^ y 



(5.17) 
(5.18) 



5.2 Automatic 0(a) improvement 



In this section we are going to prove again automatic 0(a) improvement in the 
physical basis in infinite volume. 

The target continuum theory for the fermion fields will be now 



So = d'^xtp^x) 



7^D^ + A/r 



V'(x), 



(5.19) 
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The Symanzik effective action in eq. (4.3) now has correction terms given by 



^1 = / dSCM C,{y) = ^QO,(y), (5.20) 

i 

with operators of the form 

Oo = ^AVtst'V^, Oi = ?^75rV^,F^,^, 0^ = iM^^-i^T^^. (5.21) 

where A is an energy scale of the order of the QCD scale Aqcd- The operator Oi is 
the usual clover term in the physical basis, where the Wilson term is chirally rotated. 
The operator Oq parametrizes the unavoidable mass independent 0(a) uncertainties 
in the critical mass. 

We can now repeat the same steps followed in sect. 4. We consider a general multi- 
plicatively renormalizable multilocal field that in the effective theory is represented 
by the effective field 

<l>efr = $0 + a$i + . . . (5.22) 
A lattice correlation function of the field $ to order a is given by 

($) = ($0)0 -aj A($oA(y))o + a($i)o + . . . (5.23) 

where the expectation values on the r.h.s are to be taken in the continuum theory 
with action Sq. The key point is that the continuum action (5.19) is symmetric under 
isovector flavour transformations. In particular the discrete flavour rotations 

^12 . |^(a;o, x) ^ zri'2V'(xo, x) 
[?/;(xo,x) -> -i?/'(xo,x)r^'^ 

are symmetries of the continuum action. On the contrary all the operators in eq. 
(5.21), of the Symanzik expansion of the lattice action, are odd under the the discrete 
flavour symmetry (5.24). If the field $ is a lattice representation of the flavour even 
field $0 then the second term in the r.h.s. of eq. (5.23) vanishes. To show that also 
the $1 term vanishes we have to show that an operator of one dimension higher 
than the original one but the same lattice symmetries has opposite discrete flavour 
number. This goes exactly in the same way as in the twisted basis (see sect. 4). The 
gauge action is invariant under the symmetry [11] 

P = P X [M -M] (5.25) 

where 

{U{x] /i) — > U\—x — afi] n), 
X{x) e3^'^/2^(-x) (5.26) 
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while in the fermion action the terms that break flavour symmetry T^-^ symmetry 
are odd. But in particular the lattice action is invariant under JF^'^ x T>. So the 
operators in $i will necessarily have opposite flavour number to $o- Given the fact 
that the continuum action is flavour symmetric also $i vanishes. Possible contact 
terms coming from the second term amount to a redeflnition of $i as we have 
discussed in sec. 4 and so do not harm the proof. 

The same proof can be repeated using the parity symmetry V deflned in eq. (2.34) 
instead of using the discrete flavour symmetry JF^'^. 
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6 0(a2) cutoff effects 



In this section I will be mainly concerned with the O(a^) cutoff effects of the Wtm 
formulation and in particular on the cutoff effects induced by the breaking of flavour 
and parity symmetry. Some of these require a formal description of the quantum 
mechanical representation for a lattice correlator computed with the Wtm action. 
The first part of this section will be devoted to a brief introduction of the basic 
Hamiltonian formalism, while the second part will be concentrated on the typical 
consequences of these breaking effects on physical observables. A third part will deal 
with a possible strategy in order to attenuate these O(a^) effects. 

6. 1 Hamiltonian formalism 

We assume we are in a finite volume L, which is large enough to consider hadronic 
states as point-like. This is likely to be the situation for the actual lattice simulations 
with 2fm L < 4fm. The relation among euclidean correlation function, hadron 
masses and matrix elements is given in the Hamiltonian formalism by the transfer 
matrix of the lattice theory. We have shown in sect. 2 that for Wtm the transfer 
matrix exists and it is positive definite. As a consequence of physical positivity a 
lattice correlation function admits a quantum mechanical representation in term of 
states belonging to the Hilbert space of the theory. This complete set of states can 
be chosen to be the eigenstates |n, q) of the Hamiltonian HI 

m\n,q) = E^^\n,q) , 1 = Y.\n,q){n,q\ (6.1) 

n,q 

where q indicates a set of the quantum numbers corresponding to the symmetries 
of the lattice action and n the energy level given a certain quantum number. We 
choose here the following normalization 

{n\q'\n,q) = 2E^~^5n>,n5q>,q (6.2) 

The set of quantum numbers that it is possible to use to classify the states for 
Wtm are q = {j, Pp, C, Qs, p} corresponding to the symmetries of Wtm: charge 
conjugation C, parity combined with flavour exchange Vp, charge associated to the 
residual 6^/(1)3 flavour symmetry Q3, the representation of the H{3) group of spatial 
discrete rotations j and the spatial momentum p H^ . 

In the case of plain Wilson fermions (eq. 2.28 with /iq = 0) we could replace Vp and 
Qs with standard parity and isospin symmetry V, I and Is. For Wtm the set of quan- 
tum numbers q is smaller than in the continuum. This means that certain states that 

^^Vp is not used because it is not independent from Vp and Q3. 
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in the continuum have different quantum numbers, cannot be disentangled anymore 
with Wtm. A typical example is the neutral pion where we have g^o = {j, +, +, 0, 0} 
that is certainly not distinguishable from the vacuum. More generally in the quan- 
tum mechanical representation of the correlation functions computed with Wtm, 
there will be parity violating matrix elements which vanish only in the continuum 
limit. Even if technically this is slightly more complicated that analyzing correlators 
in the plain Wilson case, the state of art for the analysis of correlation functions 
requires a multi-state fit, which is what we need to extract the masses and matrix 
elements we are interested in. Moreover given the fact that these matrix elements 
vanish in the continuum limit, it could be that numerically these contributions are 
very small at sufficiently fine lattices. To be specific I will concentrate here on 2 
examples in order to show how to proceed. 



6.2 Pion channel correlator 



We take here a finite 4-d lattice with spacing a. We want to write the quantum 
mechanical decomposition of the following correlator 

C'M^o) = -a'Y.{P\^o,^)P'm ■ (6.3) 

X 

Since we sum over x, only states with vanishing spatial momentum contribute. 
Moreover we concentrate on the trivial representation of H(3) and since we neglect 
electromagnetic interactions it is more convenient to study eigenstates of charge 
conjugation C. The relevant quantum numbers are the quantum numbers of the 
charged and neutral pion 

g^i = {0,-,+,l,0} , g^o = {0,+,+,0,0} . (6.4) 

The eigenstates of charge conjugation are 

|0,vr^(p)) = -^(|0,7r+(p)) + |0,7r"(p))) \0,Ap)) = ^(|0, vr+(p))-|0, 7r"(p))) . 

(6.5) 

These two states have the same mass because of the residual t/v(l)3 isospin sym- 
metry. Using the standard representation of a correlation function in terms of the 
local operators P^, and inserting a complete set of states, we obtain (for simplicity 
we consider the lattice time extent T oo) 

C^M^o) = - Y.{0,n\P'\n,q)—-^{n,q\P'\0,Q) . (6.6) 

^ n,q 2En 

A detailed analysis of the excited states can become rather involved. Here we concen- 
trate on the fundamental and the first excited state which is a specific contribution 
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of Wtm 

n = ^ |0,7ri(0)) (6.7) 

n = 1 ^ |l,7r^(0)7r°(0)). (6.8) 

We are here implicitly assuming that in the chosen volume L the energy level of the 
two pion state with zero relative momentum is well separated from the energy levels 
of the excited two pion state. From fig. 1 of ref. [92] this seems to be certainly the 
case for the volumes which are currently simulated. 

Contrary to what happens with standard Wilson fermions, the first excited state 
is given by a two pion contribution. The parity symmetry of the the continuum 
QCD action implies that the matrix element | (0, f2|P^|l, tt^tt*^) | must vanish in the 
continuum limit. This means it is at most an 0(a) term. This implies that the 
amplitude | (0, r2|P^|l, vr^vr^) p is an 0(o?) effect independently of whether we set the 
untwisted quark mass to zero or not l£J. The first two terms of eq. (6.6) read 



|(0,fi|pi|l,7ri(0))| 



2M^i 

+ ... , (6.9) 



\ (0,l]|Pi|l,7ri(0))P E2^{L) 

where £'271 (-^) represents the lowest energy level for the two pion state. At finite 
lattice spacing a the excited state correction increases when the quark mass goes to 
zero. In fact, up to finite size corrections, the gap is proportional to the neutral pion 
mass which, as we will see in sects. 6.4.1 and 6.5.1, vanishes at very small values of 
the twisted mass. This could be in principle a problem for a reliable extraction of 
the charged pion mass. Practically the problem is absent because |(0, fi|P^|l, tt^tt*^)!^ 
is of O(a^) and, like all the other parity violating contributions to the correlator, it 
is highly suppressed at sufficiently small lattice spacing. An additional suppression 
factor is given by '"^^^ . We remark that, independently on the lattice QCD formu- 
lation used, the gap between the fundamental and the first excited state in the pion 
channel will always go to zero in large volumes, because the first excited state will 
be at most a 3 pions state with vanishing momentum. 

Since Wtm does not preserve parity at finite lattice spacing, it would be possible 
in principle to have as excited state a charged scalar. But the scalar that can be 
created must have opposite isospin component because parity "times" isospin in 
direction 1 and 2 {Vp"^, defined in eq. 2.38) is a symmetry of the lattice action. This 
isospin component of the scalar current has opposite charge conjugation. The resid- 
ual t/v(l)3 symmetry (2.40) excludes mixing with the remaining isospin components 
of the scalar particle. 



16 



I acknowledge very useful discussions with G. Herdoiza and R. Frezzotti on this point. 
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The conclusion of this discussion is that the first excited state contributions typical 
of Wtm are either absent, like the charged scalar, or highly suppressed like the two 
pion state. This is nicely confirmed by the recent unquenched simulations performed 
by the ETM Collaboration [93]. In fig. 9, I show a typical effective mass for Nf = 2 
dynamical quarks. A long plateau in euclidean time is observed that allows a very 
precise determination of the pseudoscalar mass. We conclude that Wtm does not 
seem to make more difficult extracting the pseudoscalar mass, compared with other 
lattice actions which preserve parity symmetry. 
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Pseudoscalar channel a = 0.087 fm /j. = 0.004 
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Fig. 9. Effective mass for the pseudoscalar channel for a lattice spacing of a = 0.087 fm and 
a corresponding mass of 292 MeV. The masses using 3 different interpolating operators 
are shown. 



6.3 Charged scalar channel correlator 



The same analysis can be repeated for the charged scalar correlator. In particular we 
want to study the correlator for the charged ao scalar meson. The relevant quantum 
numbers are 

gaj = {0, +,+,1,0}, = {0, -, +, 0, 0} . (6.10) 

The 7] we consider here is the pseudoscalar isosinglet in the SU (2) framework, some- 
times labelled as 772 . The leading contributions of the charged scalar correlator 

1 ^ „-Ei^''xQ 

CssM = - Y.{0,n\S'\n,q)—-y-^{n,q\S'\0,n) , (6.11) 
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are 



n = 0^ \0,al{0)) 
n = 1 ^ |l,a^(0)7r°(0)), |1 



(6.12) 
(6.13) 



The |l,?77r^) state is the natural decay channel for the charged scalar particle. If the 
values of the quark masses in the simulations are such that the decay threshold is 
not open, this term contributes as an excited state with a gap from the fundamental 
state which could be very small. 

The parity violating term in the spectral decomposition is, analogously to the pseu- 
doscalar channel, proportional to 



with an amplitude squared of O(a^). 

In fig. 10 we see an example of a charged scalar correlator, from a recent analysis of 
Nf = 2 dynamical configurations [93]. 

Using the same argument we have use in the previous section to exclude a charged 
scalar as a possible excited state of the charged pseudoscalar, we can exclude the 
presence of a single pion intermediate state. As a confirmation of this fact there is no 
sign in fig. 10 that the correlator dips down to the value of the pseudoscalar mass. 

6.4 Isospin violation 

The fact that Wtm at finite lattice spacing breaks isospin symmetry, while retaining 
the isovector subgroup t/v(l)3, allows a splitting between neutral and charged pions. 
This is not a physical splitting and in the continuum limit of a A'f = 2 theory with 
degenerate quarks we will obtain again a degenerate triplet of pions (if we neglect, as 
we do here, the electromagnetic interactions). But at finite lattice spacing the masses 
of the charged and neutral pion are not the same. This is a very important issue to 
understand and rather surprisingly only recently numerical simulations started to 
investigate this issue [71,94,95]. Other examples of isospin violating cutoff effects are 
the splitting in the kaon sector or in the A-baryon multiplet investigated numerically 
in [68]. In general the isospin splitting is an O(a^) effect independently of the value 
of the twist angle. This can be shown using the symmetry Vp'^ (2.38) which is the 
product of parity and flavour exchange. This is a symmetry of Wtm (2.28) indepen- 
dently on the value of the twist angle. This symmetry of the lattice action implies 
that all the dimension 5 counterterms of the Symanzik effective action will have to 

1 2 

be even under Vp . For parity conserving correlation functions this automatically 



|(0,l]|g^|l,a;(O)vr°(O))p ^_^^_(^)_^^j 



(6.14) 




74 



Charged scalar channel a = 0.087 fm yu.^ = 0.004 
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Fig. 10. The same as fig. 9 but for the charged scalar channel and with the two best 
correlators. 



implies that their 0(a) insertion will contribute in the same way to the correlation 
functions independently on their flavour structure. Thus the isospin splitting has to 
be at least of O(a^) independently on the value of the twist angle. In other words the 
0(a) affecting the neutral and charged pion correlation functions out of full twist are 
the same. This is confirmed by the results obtained in W^PT [83,96] and discussed 
in the following section. 



6.4-1 Pion splitting 

To fix the notation we recall some basic definitions. The charged pseudoscalar den- 
sities are given by 



ri ± it' 



(6.15) 



and a possible interpolating field for the neutral pion, with u = 7r/2, is the scalar 
current 



(6.16) 
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The charged and neutral pseudoscalar masses can be extracted by the following 
correlators 

a+(xo) = -a'Y.{P^i-^)P^m CAxo) = f E(5°(^)^°(0))- (6-17) 

X ^ X 

If we separate the flavour components and we perform the fermionic contractions 
we obtain 



tr 



ao(a;o) = -E{( 

^ X 

+ tr x) 



G'„(0,x)G,(x,0) 



tr 



G,(0,x)Grf(x,0) 



+ tr 



Gd(x,x] 



Gu{0,0) 



+ tr 
+ tr 



G^{x^ X ^ 
Gd(x,x) 



Gd{0,0) 



(6.18) 



where Gu,d{x,y) are the fermionic propagators for the the up and down quarks. 

The neutral pion correlator (6.18) contains fermionic disconnected diagrams (the last 
four terms in eq. 6.18), which are notoriously more difficult to compute on the lattice. 
One way to circumvent this problem is to consider only the fermionic connected part 
(the first two terms of eq. 6.18). The question would be now if it is still possible 
to interpret the fermionic connected piece as coming from a correlation function 
between local operators which represent a pion field. The answer is positive if we 
consider the valence quarks of the theory discretized with the so called Osterwalder- 
Seiler (OS) action [97,98]. To be specific we consider now the following model: the 
sea quarks are described by the standard Nf = 2 degenerate Wtm action (2.28) 
while the valence quarks, that we call u and d' have the following actions 



X 

SS = a' E {d'i^) [Dw + m„ + ^/il^^Ts] d'{x)}. 



(6.19) 



(6.20) 



These two quarks have the same mass /i^ = /i^/, thus they have identical actions, dif- 
fering only because they have different fiavours. We can then define as interpolating 
field for the neutral pion 



(vr°)'fx) 
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u{x)u{x) — d'{x)d'{x) 



(6.21) 



It is easy to see that, performing the axial rotations (2.6) corresponding to the 
actions (6.19,6.20) 



u[xj — > exp(-z|75)u(x) 
u{x) — > u{x) exp(— 2|75), 




exp(-i|75)d'(x) 
c!'(x)exp(-if75) 



(6.22) 
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with uj = 7r/2, this is a vahd interpolating field for the neutral pion. Performing 
standard Wick contractions we obtain 



X 

where the disconnected terms cancel because the u and d' quarks have identical 
propagators. Since parity "times" fiavour exchange is a symmetry of the lattice 
action we have 

7oGd(0,xp)7o = G«(0,x) = (7^/(0, x), with xp = (xo,-x). (6.24) 

The connected part of the neutral pion correlator with Wtm (6.18) is equal to the 
correlation function (6.23) for a model with OS valence quarks. This argument tells 
us that the connected part of eq. (6.18) is not the tt" propagator with Wtm, but it 
is still a 7r° propagator with valence quarks having the same OS actions (6.19,6.21). 
In the continuum limit the two neutral pion correlators will coincide. 

The OS action can be used only as a discretization for valence quarks because a par- 
ity violating and isospin singlet counterterm like FF could be otherwise generated 
in the process of renormalizing the theory (see sect. 2). This freedom of choosing a 
different lattice action between valence and sea quarks will be further used to ease 
the renormalization of local operator (see sect. 7). Here we just use this freedom as 
a tool to interpret the fermionic connected part as a theoretically valid expression of 
a correlation function between local operators. We remark that the fermionic con- 
nected part is not the neutral pseudoscalar meson of Wtm, but it is an interesting 
quantity to study with precise data on its own, in view of a possible use of mixed 
actions (the OS action for the valence quarks and Wtm for the sea quarks). More- 
over in the quenched model the fermionic determinant is neglected anyhow, and this 
puts the two models (quenched OS and quenched Wtm) on the same footing. 

We also note that in the same model we could define a charged pion made up of 
a u and d' quark. In this case the correct interpolating field would be ud' that has 
identical correlation function with the neutral pion field (vr°)'. This leads to a model 
with NO isospin splitting among the pions. The drawback of this choice would be, 
apart from O(a^) unitarity violations, that this charged pion field is not the field 
that is protected by the exact lattice Ward identity (7.2). 

The choice of valence actions different from the sea action allows freedom to solve 
some problems which appear in the unitary formulation. We will see in sect. 7 how 
this freedom can simplify the renormalization of local operators. It remains to be 
seen how the O(a^) unitarity violations affect the theory and the approach to the 
continuum limit. This is an open question which could be answered only performing 
detailed scaling tests. 



Gu{0,x)Gu{x,0) -ti Gd'{0,x)Gd'{x,0))}, (6.23) 
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In [94] , to which I refer for all the technical details, a pilot quenched study has been 
performed to study flavour breaking effects with Wtm. In this paper the scaling 
behaviour is studied in the quenched model of the pion splitting using either only 
the connected correlator (YMpg)Qg — (Mpg)^), or including also the disconnected 
terms ((Mpg)^ — (Mpg)^j^. The results show O(a^) scaling violations for both the 
pion splittings with indications that the neutral pseudoscalar meson for Wtm (with 
the inclusion of the disconnected correlator) has reduced splitting with the charged 
meson, within the rather large statistical errors. It is possible to give a very rough 
estimate of the pion splitting ro((Mpg)^ — (Mpg)^) ^ c(a/ro)^ with c ^ 10 (with 
large errors). Comparing to a quenched simulation for naive staggered fermions with 
Wilson gauge action [99] , one finds a similar size of the flavour splitting encountered 
for the pion mass at a similar lattice spacing with a value c ^ 40. For dynamical 
improved staggered fermions a value of c ~ 10 has been found [100]. 

The pion splitting has been investigated in the quenched model in [71], where a 
comparison between the impact on the methods used to determine the critical mass 
has been made. The outcome of this study is that if one considers only the fermionic 
connected diagrams (OS neutral pion) the inclusion of the clover term into the Wtm 
lattice action (4.18) reduces the pion splitting in comparison with the Wtm action 
without clover term. This result is rather interesting but not conclusive, in fact the 
addition of the disconnected terms can have a tremendous impact as we have just 
seen in the quenched model. Moreover recent dynamical simulations [69] indicate 
that the inclusion of the disconnected terms changes the sign of the pion-splitting. 
To be precise the outcome is that a first analysis at a value the pseudoscalar mass 
Mps — 300 MeV, taking the disconnected contribution in the neutral channel fully 
into account, shows that the neutral pseudoscalar meson is about 20% lighter than 
the charged one. Expressed differently one obtains rQ((Mpg)^ — (Mpg)^) = c(a/ro)^ 
with c = —4.5(1.8). This coefficient is, in absolute value, a factor of 2 smaller than 
the value found in quenched investigations [101]. If one considers only the connected 
term, i.e. considering OS valence quarks, the charged pion is lighter or in other words 
the sign of c is positive. 

From the theoretical side a tool to address and understand this problem is again 
WxPT. Adopting the same power counting scheme used in sect. 4 given in eq. (4.67) 
the result [83,96] at NLO order is 

Mi-Mj.^-2.V-4-,. (6.25) 

The splitting is as expected an O(a^) effect and vanishes on the Wilson axis (/ir = 0) 
as it should. It is moreover necessarily even in /iq, from the combined violation of 



^ In both computations the "PCAC method" is used for the determination of the critical 
mass. 
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parity and isospin argument given in sect. 6.4. The splitting is maximized when 
m' = as it should be since this is the condition that minimizes the parity violation. 
The sign of the LEG w' tells us which pion is heavier. We will see in sect. 6.5 that w' 
also determines the nature of the vacuum structure of the theory around the chiral 
point. 

Given the non-degeneracy among the pion multiplet, several issues would have to 
be investigated in the future. The usage of a mixed action approach with OS va- 
lence quarks (or another isospin conserving action like the overlap discretization) 
eliminates the problem of isospin splitting at the valence level. Of course the isospin 
breaking cutoff effects of the sea Wtm action will still generate a splitting, but only 
among the virtual pions. One possible consequence of this phenomenon is the change 
of the finite size effects in a physical quantity induced by the pion cloud: the domi- 
nant contribution comes from the lightest pion of the theory at finite lattice spacing, 
that, as we have seen, is not necessarily the charged one I . 

6.4-2 Kaon splitting 

With an interest in the phenomenology of hadrons built of n, d and s quarks, it 
is useful to explore Wtm including heavier quarks. We have already seen in sect. 2 
that there is no unique way to introduce the s quark into the calculation. The 
method investigated in ref. [95] considers a pair of quark doublets (m, d) and ("c", s), 
following the proposal of ref. [39] discussed in sect. 3.2. In particular the partner 
of the s quark does not play an active role and should not be thought of as the 
physical charm quark; moreover with this method no mass splitting is introduced 
within either doublet. 

To be specific one introduces Wtm with two degenerate quark doublets. 



referred to as the light and heavy doublets respectively. 

As we have seen in sect. 2 this construction can be extended to include non- 
degenerate quarks in a single doublet [39]. Even if with this approach the fermion 
determinant is not real, and so it would not be practical to perform dynamical simu- 
lations, this action could still be useful in the spirit of using different lattice actions 
for the valence and the sea quarks. 

The two-doublet lattice action is simply a block- diagonal version of two copies of 
I thank C. Michael for bringing also my attention to this issue. 




(6.26) 
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the one-doublet theory (2.28). 




where and V* are the usual covariant forward and backward lattice derivatives 
respectively, and 

° ]. ' ]■ (6.28) 

mh,ol2j \ IJ^h,qr:iJ 

The method used in [95] to fix the critical mass is the "cu^ method", with the 
untwisted bare quark mass tuned so as to achieve full twist at each particular twisted 
mass value. 

NLO WxPT [95] predicts a formula for M|-o — M^± very similar to the formula for 
the pion splitting (6.25). In particular this formula, analogously to the pion case, 
predicts for the kaon splitting at NLO no dependence on the twisted quark mass if 
the twist angles for the light and the heavy sector are u = 7r/2 for all the values of the 
twisted mass ("cua method"). The same formula, again analogously to the pion case, 
predicts that the kaon splitting should depend linearly in and it should vanish in 
the continuum limit. Modulo the unknown higher order effects, figs. 11 and 12 are 
in reasonable agreement with these expectations. In fig. 11 are plotted the masses 
squared of the charged and neutral kaons, i.e. the ground state pseudoscalar mesons 
containing one s quark from the heavy doublet and one u or d quark from the light 
doublet, as a function of the sum of the valence twisted quark masses. Fig. 12 shows 
the lattice spacing dependence of the squared mass differences for four mass values. 
To summarize: the approximate mass independence of the kaon splitting is evident 
from fig. 11, and the dependence on in fig. 12 is approximately linear, though 
a linear fit misses the massless prediction at a = by a few (statistical) standard 
deviations. 

It should be noted from fig. 12 that even at the smallest lattice spacing, the mass 
splitting of rrixo — m/^± ~ 50 MeV is significant relative to the kaon mass itself. 
However, in terms of the difference of mass squared, these results are consistent 
with the pseudoscalar meson mass sphttings discussed in ref. [94]. 

The splitting in the Kaon system that we have just discussed involves only twisted 
quarks. Different O(a^) cutoff effects appear if one studies the splitting between pseu- 
doscalar correlators computed with twisted and untwisted quarks. One can consider 
a twisted (m, d) doublet, and flavour- singlet Wilson strange (and charm) quarks. 
This approach could be a viable one for doing full dynamical simulations. On the 
other side automatic 0(a) improvement would be lost and there would be the need 
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Fig. 11. Pseudoscalar meson mass squared as a function of the sum of quark and antiquark 
twisted mass parameters. Subscripts I and h indicate the light and heavy doublets. Results 
are taken from ref. [95] which uses the lva method and and from ref. [94] which uses the 
PCAC method. The results from [94] have equal masses for the quark and anti-quark. 



to introduce suitable improvement coefficients. In ref. [72] a precise investigation has 
been performed, in the quenched model, of this particular splitting between pseu- 
doscalar correlators induced by the fact that twisted and untwisted actions have 
been used for the quarks. The twisted and untwisted actions of [72] are both clover 
improved, and all the relevant improvement coefficients have been used to improve 
the correlators. The twisted and untwisted quark masses are always matched in or- 
der to have the same continuum QCD theory with meson masses around the Kaon 
mass. The critical mass rricr used in [72] has been taken from refs. [43,73,74] and it 
has been recomputed at one lattice spacing. The quantities studied in this paper are 
the pseudoscalar masses Mtt, Mtw, ^ww and decay constants Fu, Ffw, Fww for 
mesons which are respectively made of a doublet of twisted quarks tt, a twisted and 
untwisted quark tW and two untwisted quarks WW. The meson made of twisted 
quarks is the charged pseudoscalar correlator of eq. (6.17). The outcome of this pre- 
cise quenched study is that the splitting between pseudoscalar meson masses and 
decay constants ranges, at the coarsest lattice spacing a = 0.093 fm, between 5-13% 
and it goes to zero in the continuum limit. The splitting seems to depend also on 
the value of the improvement coefficient ca which is needed to improve the axial 
current with Wilson fermions. In this paper the issue of the pseudoscalar splitting 
induced by the twisted mass term is not addressed. 
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Fig. 12. The difference between charged and neutral squared pseudoscalar meson masses 
as a function of squared lattice spacing, for selected values of the charged meson mass. 
Plot taken from ref. [95] 



An exploratory dynamical study with two light degenerate quark and two heavier 
non-degenerate quarks {Nf = 2 + 1 + 1) has been carried out in ref. [41], following 
the approach explained in sect. 2 with off-diagonal flavour splitting [38]. As we have 
already explained this approach has the big advantage of being well suited for dy- 
namical simulation. Moreover it has been noted in [41] that with this formulation 
the masses in the kaon doublet (and D-meson doublet) are exactly degenerate. This 
follows from an exact symmetry of the lattice action defined in eq. (3.20). To un- 
derstand how this works, we first write two interpolating fields for the kaon doublet 
in the physical basis (see sect. 5): 

K+{x) = m(x)75s(x), K°{x) = d{x)j5s{x). (6.29) 

Performing the rotations (3.2,3.5) we obtain in the twisted basis 

i 1 1 - i - 

K^{x) = --j=u{x)c{x)+-j=u{x)-f5s{x), K^{x) = --j=d{x)-f5c{x)--j=d{x)s{x) 

(6.30) 

Both the action for degenerate u and d quarks in eq. (2.28) and for non degenerate 
s and c quarks in eq. (3.20) are symmetric under parity "times" flavour exchange 
Vp defined in eq. (2.38). Under this symmetry transformation the single flavoured 
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quarks transform as 



M(Xo,Xj 

rf(xo,x) 
c(xo,x) 
[s(xo,x) 



z7os(xo, - 
z7oc(xo, - 



-x) 
-x) 
-x) 



(J(xo,x) 
c(xo,x) 
s(xo,x) 



The effect of tliis exact symmetry transformation is 
i^+(xo,x) ^ J^O(xo,-x), ir°(xo,x) 



-id{xo, -x)7o 
-iM(xo, -x)7o 
-zs(xo, -x)7o 
-zc(xo, -x)7o 



(6.31) 



(6.32) 



Hence the equahty of the masses within kaon doublets follows. It is straightforward 
to repeat the same argument to show that also the D-doublet is mass degenerate. 



6.4-3 Baryon splitting 

In principle O(a^) isospin breaking cutoff effects can appear also in baryon cor- 
relators. Here we concentrate on the spin-parity (1/2)^ and (3/2)^ baryons and 
we always consider Wtm with two degenerate quarks as specified by the action in 
eq. (2.28). 

We write now the interpolating fields for the baryons in the physical basis 



^ = eABc[u\C"^l^dB\uc, (6.33) 

^ = eABc[dlC-'^,UB]dc, (6.34) 

A'l^ = eABc[uAC~'^lkUB]uc, (6.35) 

2 1 

= -^eABc[uAC~'^'ykdB]uc + -^(-ABcWAC~^lkUB\dc, (6.36) 

2 1 

A° = —^eABcWA'^~^lkUB\dc + -j^eABcWA'^~^lkdB\uc, (6.37) 

Afc = eABc\dl<0~^lkdB\dc, (6.38) 



where C is the charge conjugation matrix defined in app. B. To get the corresponding 
interpolating fields in the twisted basis it is enough to apply the rotation in eq. (2.6) 
to the quark fields. In app. F, I show how this works with the proton field. Here I just 
want to prove that the proton and the neutron with Wtm are degenerate even at 
finite lattice spacing. This was shown in ref. [68] and here we give a slightly different 
argument based on the symmetries of the lattice action. In particular we know that 
parity "times" flavour exchange Vp, defined in eq. (2.38) is a symmetry of Wtm. 
This symmetry is independent on the basis of choice, in fact also the action in the 
physical basis (5.3) is invariant under the Vp symmetry transformation. Applying 
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this symmetry to the baryon fields we obtain 

^(xo, x) ^ -27o^(xo, -x) (6.39) 
A++(xo,x) ^ -z7oA^(xo,-x) (6.40) 

A+(xo, x) ^ -z7oA°(xo, -x). (6.41) 
If we consider a baryon field SS and the corresponding correlation function 

a3^(i^(a;o,x)if(0,0)), (6.42) 

X 

we conclude that proton and neutron are degenerate as are A^^ with A^ and A^ 
with A°. There is on the other hand an O(a^) splitting between A^^ and A^. This 
splitting it is also numerically simpler to study than the pion splitting since it does 
not require computation of disconnected diagrams. This was in fact noticed and 
studied in ref. [68] in the quenched model. Figs. 13 and 14 show the results for 
the splitting in physical units in the A channel, as obtained respectively from the 
"Wilson pion" definition of the critical mass and the "tu^ method". No significant 
difference is found between the two definitions of full twist for this observable. Based 
on the larger statistics of fig. 13 (1000 configurations rather than only 300), there is 
some evidence that m/^\+- — m^+,o decreases as a ^ 0, as expected. 
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Fig. 13. Flavour splitting within the A(1232) multiplet as a function of the pseudoscalar 
meson mass squared, for two lattice spacings: a = 0.123 fm (o) and a = 0.093 fm (•). The 
critical mass was fixed using the "Wilson pion" method. 
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Fig. 14. Same as fig. (13) but with the critical mass fixed using the "a;^ method". 

The extension of W^PT, with a twisted mass term, to include the baryons has been 
worked out in ref. [102]. In this work, beside showing how to include the baryons 
in the effective theory, several interesting results have been obtained. It has been 
confirmed automatic 0(a) improvement at full twist for baryon masses, it has been 
discussed the splittings among the delta-resonances and a mass splitting formula 
has been given for a genric twist angle uo. 



0.2 



6.4-4 Conclusions 

Isospin breaking cutoff effects are most probably the most delicate issue concerning 
Wtm. In the pseudoscalar sector they are similar in magnitude with the staggered 
taste violating cutoff effects, and they are even smaller if we include the discon- 
nected diagrams in the pion sector. The possibility to use the OS action for the 
valence quarks partially mitigates this problem. The mixed action approach, while 
removing completely the isospin breaking at the valence level, cannot cure com- 
pletely the problem. For physical quantities which involve pion scattering, isospin 
breaking cutoff effects are an issue which has to be further investigated. In the baryon 
sector isospin breaking cutoff effects seem to be under control. Moreover choosing 
as a lattice action for the heavy non-degenerate quarks the action proposed in [38] 
the splitting in the K and D sector is absent. 
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6.5 Phase structure of Wilson fermions 



In the first part of tliis sect. 6, we liave discussed tlie O(a^) isospin breaking cutoff 
effects. Now we want to analyze the impact of the O(a^) discretization errors the 
shape of the chiral phase diagram. The O(a^) chirally breaking cutoff effects change 
the form of the chiral phase transition. Thus the knowledge of the shape of the chiral 
phase diagram is an important prerequisite to perform large scale simulations with 
a given lattice action. In this section I show that both numerically and analytically 
there is strong evidence that Wilson type fermions show a peculiar phase diagram 
around the chiral point. The structure of such a phase diagram is due to lattice 
artifacts. We will see in sect. 8 that the approach to the chiral and continuum limit in 
numerical simulations shows a substantial slowing down in algorithmic performance, 
so in general a careful continuum limit and chiral extrapolation is required in order 
to compare lattice results with experimental measurements. The knowledge of the 
position of the phase transition point as a function of the quark mass and the 
lattice spacing, becomes then very important, given also the fact that it is not clear 
whether HMC-like algorithms can correctly sample the configuration space in such 
extreme situations. In particular the shape of the phase structure is such that there 
is a minimal pion mass that can be simulated at fixed lattice spacing. This bound 
depends on the details of the Wilson-like action used in the simulations (e.g. on 
the gauge action). Lattice simulations are performed in a finite volume and strictly 
speaking there are no phase transitions in a finite box, but if the lattice size is large 
enough, the shadow of the would-be phase transition could still be visible in the 
numerical simulations. 



6.5.1 WxPT analysis 

Before going into the study of the phase diagram of Wilson-like fermions let us 
recall which is the situation in continuum QCD. In the continuum chiral symmetry 
is spontaneously broken in the massless limit, and the order parameter of such a 
transition is the chiral condensate. The structure of the phase diagram is thus simply: 
a singular point at the origin of the "mass plane" where we plot the untwisted quark 
mass along the horizontal axis and the twisted mass along the vertical axis (see 
fig. 15). The aim of this section is to understand how this picture is modified by 
the inclusion of lattice artifacts in the chiral Lagrangian. We have already discussed 
in sect. 4.3.2 the construction of the chiral Lagrangian including lattice spacing 
effects. In particular our starting point was the power counting scheme specified 
in eq. (4.67). We have seen in sect. 4.3.2 that at LO the inclusion of 0(a) lattice 
artifacts can be fully reabsorbed in a shifted quark mass m', defined in eq. (4.71). 
We can then immediately conclude that if we are in the power counting scheme 
given by eq. (4.67) the phase structure is continuum like. To start to be sensitive to 



86 



the modifications of tlie pliase structure in presence of lattice artifacts we have to 
further lower the quark masses at fixed lattice spacing. In particular the appropriate 
power counting is given 

1 > m',/iR,p^a2 > . . . (6.43) 

This region of quark masses is usually called the Aoki region because it is the region 
where a possible Aoki phase appears [12]. The Aoki phase, as we will see below, is 
characterized by a region in the bare parameters space where parity and flavour are 
spontaneously broken. 

Since we are changing the power counting scheme, the careful reader could wonder 
if this does not imply a breakdown of the expansion. The reason why this is not the 
case is that there is a rearrangement in the ordering of the coefficients of the chiral 
Lagrangian. It is possible to show [78] that the neglected terms of the expansion are 
at most of order ~ a'^, and thus suppressed by one power of a. Loop corrections are 
also suppressed, since they are quadratic in m! and /xr (up to logarithms) and thus 
~ a^. The reordering of the expansion is possible only because the leading order 
discretization error has exactly the form of a mass term and so can be completely 
absorbed into m', to all orders in the chiral expansion [78]. 

Upon neglecting the derivative interaction terms, as we are interested in the vacuum 
state, the potential of the effective chiral Lagrangian at LO [55,83,96, 103] is given 
by 

= -^(S + St) + ^(S + St)2 + |(z(S - St)r3) (6.44) 

where the explicit forms of the coefficients, and their sizes in our power-counting 
scheme, are 

Cl = 2Bofm' ~ m', 

p2 / 2 2 

C2 = — / w a ~ a , 

C3 = 25o/Vr ~ /^R, • (6.45) 

All the Cj coefficients are independent at non-zero lattice spacing and the only extra 
term introduced by the twisted mass is the one proportional to C3. 

In the Aoki region the coefficient C2 is of the same size as the coefficients ci^a; all 
are of O(a^). This implies a competition of the mass terms and the O(a^) term in 
the shape of the potential that causes a non-trivial vacuum structure. Minimizing 
the potential (6.44) gives rise to two possible scenarios [78,83,96,103] for the phase 
diagram of Wilson-like fermions: 

• the Aoki scenario [12]; 

• the Sharpe-Singleton scenario [78]. 
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Note that both m' and /xr must be of O(a^) in order for such competition to occur; 
if either m' or /xr is of 0{a) then one is in the continuum-hke region. I would 
hke to remark that the presence of two possible phase structures was foreseen by 
Creutz [104]. 

The order parameter of the continuum chiral phase transition is the chiral con- 
densate. The condensate in W^PT is determined by the minima of the potential 
energy (6.44). We parametrize the chiral field in the standard way: S = So + ■ r" 
with real Sq and Sq satisfying So ■ Sq + ■ = 1, so that So, S^ e [—1, 1]. Sim- 
ilarly, the condensate (the value of S at the minimum of the potential) is written 
Sq™''' = Sq"^^ + iS^*"^ ■ r". In the Aoki region the potential is then 

= -ciSo - C3S3 + C2^l . (6.46) 

We have now to compute the values of Sq"*'' and S3™'' . We will not enter in the details 
of the computations that are nevertheless rather simple. One finds that there two 
possible solutions depending on the sign of C2 (or w'). Thus it is the sign of C2 (or 
w') that determines the possible scenario for the chiral phase structure of Wilson- 
like fermions. We summarize here the results illustrating them with plots [96]. For 
convenience we define the following rescaled mass variables 

m = A = - ^7a' (6.47) 

\C2\ a''\w\ 

n = — = ~ hr a . (6.48) 

\C2\ a^\w'\ 

They are of 0(1) in the region of interest. 

In fig. 15 I plot the lines corresponding to a first order phase transition which have 
second order phase transitions end points, where the x-axis is m and the y-axis is n. 

I first discuss the structure of the phase diagram for C2 positive or negative and then 
I will describe the physics one finds in the two situations. 

In the left panel of fig. 15 the Aoki scenario that appears if C2 > {w' < 0) is 
depicted. In this scenario along the Wilson axis m there is a region (indicated by 
the thick line) where parity and flavour are spontaneously broken: this is the Aoki 
phase [12]. At the endpoints of this line the three pions are massless. Inside the 
Aoki phase the charged pions stay massless while the neutral one becomes massive. 
A non-zero value of the twisted mass washes out the Aoki phase introducing an 
explicit breaking of flavour and parity symmetry. 

In the right panel of fig. 15 the Sharpe-Singleton scenario that appears if C2 < 
{w' > 0) is depicted. In this scenario the first order phase transition line extends 
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into the twisted direction. The transition ends with a second order phase transition 
point, where the neutral pion mass vanishes. 
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Fig. 15. Left panel: the phase diagram of Wilson diagram according to Wx PT for C2 > 0. 
Right panel: as in the left panel but for C2 < 0. The x-axis is proportional to m' /a^ and 
the y-axis to /xr/o^ (6.47,6.48). 

We go now to study the mass dependence of the condensate and of the pion masses 
inside the Aoki region. 

We begin with results for C2 > 0. Fig. 16 shows the behaviour of the identity com- 
ponent of the condensate, Sg'"^ = (S^''"-*)/2, for n = 0, 1, 2 and 3 as a function of m. 
The corresponding pion masses are shown in fig. 17. In the untwisted theory (n = 0) 
there are second order transitions at m' = ifc ° ^ , as shown by the discontinuity in 

the derivatives respect to m in Sq™"* and the vanishing of the pion masses. In the Aoki 
phase, with £3™'' 7^ 0, the SUy{2) vector flavour symmetry is spontaneously bro- 
ken into [/v(l)3, and correspondingly there are two Goldstone bosons between these 
second-order points. Once n is non-vanishing, however, the transition is smoothed 
out into a crossover, and the pion masses are always non-zero. Flavour is broken for 




Fig. 16. The global minimum of the potential, Sq , as a function of m, for C2 > and 
n = 0, 1, 2, 3. 
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Fig. 17. Mass of the pions as a function of m, for C2 > and n = 0, 1, 2, 3. 

all m, with the charged pions lighter than the neutral pion by O(a^), as given by 
eq. (6.25). For the special case of m = (full twist), Sg"^-* vanishes and (Sg*"-*)^ = 1, 
so the mass-squared splitting is 2a^|w'| for all n (which becomes a difference of 2 in 
the units in the plots). 

If we keep —2 < m < 2 fixed, and we lower n, it is possible to pass through the Aoki 
phase. When we change the sign of twisted mass, i.e. passing through n = 0, there is 

a first-order phase transition, since Sg"^'* jumps from +\Jl — (Sq™'')^ = y^l — 
to - (4™^)^ = - ^1 - mV4. 

We now consider the case where C2 is negative. Fig. 18 shows Sq™''' and the pion 
masses as a function of m for fixed values of n. Fig. 18(a) and fig. 18(b) show the 
results for /ir = 0. The condensate jumps from Eq = 1 (and thus Eg'"'* = 1, S^™) = 0) 
for m > to So = -1 (and thus = -1, S^"") = 0) for m < 0. This is a first order 
transition without flavour breaking, so all pions remain massive and degenerate. 

The rest of fig. 18 shows what happens at non-zero twisted mass. The effect of /iR 
is to twist the condensate, so that there is a non-zero rg component t!"^\ There is. 
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Fig. 18. Global minimum Sq and pion masses as a function of m, for C2 < and 
n = 0, 1, 2, 3. The dashed lines are for 7ri^2 and the solid lines are for vra. 

however, still a first order transition at which Sg™^ flips sign between ±(1 — n/2) 
(assuming n > 0). The neutral pion is now lighter than the charged pions due to 
explicit flavour breaking. The neutral pion has a mass M^^ = 2a^|u;'|(l — at 
the transition, while, as noted above, the charged pions have a n independent mass 



given by m: 



2a \w'\. The transition weakens as Inl increases, and ends with a 
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Fig. 19. Mass of the pions as a function of n, for C2 < and m = 0, 1, 2, 3. 

second order transition point at n = ±2, at which the neutral pion is massless [see 
Fig. 18(f))]. For larger n the transition is smoothed out. Note that, once away from 
the transition, for m = the mass-squared splitting between charged and neutral 
pions is still 2a^|w'| and independent from n. 

These plots illustrate the general fact that the C2 < case can be obtained from 
that with C2 > by a 7r/2 rotation and appropriate redefinitions. Indeed, the results 
for C2 > can alternatively be viewed as the plots for C2 < at fixed m with n 
varying, and vice versa, with the exception of the charged pion masses, which differ 
by a constant offset of 2|c2|//^. To illustrate this latter point I plot, in fig. 19, the 
pion masses for C2 < as a function of n for fixed values of m. Comparing to fig. 17, 
we see the equality of the neutral pion masses and the constant offset in the charged 
pion masses. 

6.5.2 Numerical results 

The analysis carried out with W^PT is very predictive but cannot choose between 
the two possible shapes of the potential. The occurrence of one of the two scenarios 
depends on the sign of the coefficient C2 proportional to the O(a^) terms of the chiral 
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Lagrangian, and this coefficient C2, similarly to the other LEG, it is not determined 
by chiral symmetry, but it depends on the details of the lattice action one is using, 
like the gauge action, the presence in the lattice action of the clover term and of 
course it depends on the value of the bare gauge coupling, i.e. of the lattice spacing. 

The lattice spacing dependence of the phase structure can be studied in detail us- 
ing lower dimensional models, allowing a better analytical control over the com- 
putations. An analysis with Wilson fermions of the two dimensional Gross-Neveau 
model [105] indicates that indeed both the scenarios describe the phase structure of 
Wilson fermions depending on the value of the couplings of the model. The analysis 
shows that at strong coupling there is an Aoki phase while at weak coupling the 
Sharpe-Singleton scenario sets in. This analysis has been recently extended for the 
twisted mass case [106], indicating even more complicated structures, like a coexis- 
tence of the two scenarios for certain values of the bare couplings. 

In QGD the situation is more involved. The prediction of the existence of an Aoki- 
phase was made years before the analysis of W^PT [12, 107, 108]. The existence of a 
region of bare parameters where parity and flavour were spontaneously broken was 
a possible mechanism to generate a massless pion at finite lattice spacing. Quenched 
studies [109-111] found evidence of such a phase structure. Recently a quenched 
computation [94] of the charged and neutral pseudoscalar masses using Wtm has 
shown that the charged pion is lighter than the neutral, indicating that in the 
quenched case we are in an Aoki scenario. 

On the other hand it is interesting to revisit some old controversial results about 
the chiral phase structure of Wilson fermions, using the standard Wilson gauge ac- 
tion (2.26). In [112, 113] from a finite temperature study there was an indication 
of difficulties in observing a phase with spontaneous breaking of flavour and parity 
symmetry (Aoki phase) at /3 > 4.8. In [114] the MILG collaboration found a sur- 
prising bulk ffist order phase transition at /5 ~ 4.8. In a recent investigation [115] of 
the Aoki phase for values of the coupling j3 < 5 the authors con&m evidence for an 
Aoki phase but only for values (3 < 4.6. For values of > 5 with Wilson plaquette 
gauge action a systematic investigation was missing. 

This gap has been filled in a set of publications that have thoroughly investigated 
the chiral phase structure of Wilson-like fermions. In [116] the ffist study of Wtm 
with A'^f = 2 dynamical fermions was performed and rather surprising results were 
found: the existence of the Sharpe-Singleton scenario. 

The action used in this study is Wilson gauge action combined with Wilson fermions 
with and without twisted mass. In particular at a lattice spacing of a ^ 0.16 fm, 
strong evidence of a ffist order phase transition was found for a rather large range of 
values of twisted masses going from zero twisted mass to /iq ~ 100 MeV. This study 
reveals also that the phase transition tends to disappear on increasing the value of 
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/iq, it persists for /iq = and it is volume independent. 



This study was extended in [117]. In fig. 20 I show scans of the average plaquette 
and mpcAC fixed twisted mass (with /iq roughly fixed in physical units) for three 
decreasing lattice spacings (at increasing P values). The plaquette and mpcAC both 
have a discontinuity, and show hysteresis. Both effects decrease as one approaches 
the continuum limit, qualitatively consistent with expectations. The fact that mpcAc 
has a minimum away from zero is a manifestation of the non-zero minimum in the 
pion masses. 
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Fig. 20. Average plaquette and mpcAC plotted vs. (2k) ^ = mo + 4 for fixed fiq [117]. 



Our present understanding of the lattice QCD phase diagram with Wilson gauge 
action and Wilson-like fermions can be summarized as follows. For values of the 
lattice spacing much coarser than a = 0.15 fm, there is an Aoki phase [12, 115, 
118]. For smaller values of the lattice spacing, a first order phase transition appears 
(Sharpe-Singleton scenario) [84,116,117,119] that separates the positive quark mass 
from the negative quark mass phase, and extends into the twisted direction. This 
first order phase transition is reminiscent of the continuum phase transition when 
the quark mass is changed from positive to negative values with the corresponding 
jump of the scalar condensate which is the order parameter of spontaneous chiral 
symmetry breaking. The generic phase structure of lattice QCD is illustrated in 
fig. 21 and discussed in refs. [84,116,119]. 



6.5.3 The gauge action 

In the Sharpe-Singleton scenario, the pseudoscalar mass mpg cannot be made ar- 
bitrarily small. When lowering the quark mass from the twisted direction there is 
a minimal pion mass accessible with numerical simulations given directly by the 
extension of the first order phase transition line, even if the twisted mass provides 
a sharp infrared cutoff in the sampling performed by the simulation algorithm. 

An important result of W^PT is that both the size of the phase boundaries and the 
isospin splitting for pions are determined by the same parameter, w'. It thus makes 
sense to try and tune the gauge and fermion actions to reduce \w'\. Note that this 
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Fig. 21. Current knowledge of the Wilson lattice QCD phase diagram as function of the 
inverse gauge coupling (3 oc the hopping parameter k = (8 + 2amo)~^ and the twisted 
mass parameter /Xq. 

tuning is not the same as a systematic improvement program, but it is nevertheless 
very important. 



In [117] the lattice spacing dependence of the first order phase transition with Wilson 
gauge action has been studied, taking, as a measure of its strength, the gap between 
the two phases in the plaquette expectation value and in the PCAC quark mass. The 
qualitative estimate for the lattice spacing, where a minimal pion mass m,r — 300 
MeV could be reached without being affected by the first order first transition is 
0.07-0.1 fm. 

It is suggestive to interpret, at the microscopic level, the occurrence of the Sharpe- 
Singleton scenario with a massive rearrangement of the small eigenvalues of the 
Wilson-Dirac around the first order phase transition line. This rearrangement could 
be suppressed by the use of a renormalization group improved or O(a^) improved 
gauge action, and indeed results from [120] indicate that metastabilities in the av- 
erage plaquette observed for iVf = 3 dynamical Wilson fermions with a clover term 
can be suppressed replacing the Wilson gauge action with the Iwasaki action [121]. 

The dependence of the phase diagram on the gauge action used and on the lat- 
tice spacing has been studied in a set of papers [84, 116, 117, 119] (see also [122] 
for a detailed summary of these results). The gauge actions studied so far can be 
parametrized by 



^ 6 



bo tr(l - P^^^(x; z/)) + 6itr(l - P^^^(x; (6.49) 
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with the normahzation condition bo = 1 — 861. In particular three gauge actions have 
been investigated: 

• Wilson action [1] =^ bi = 

• tree-level Symanzik action [123] ^ ^1 = 

• DBW2 action [124] ^ b^ = -1.4088 

The result of these studies is summarized in fig. 22 which shows that the discontinu- 
ities in the plaquette are significantly reduced going from the Wilson to the DBW2 
gauge action. 

By varying the coupling bi which multiplies the rectangular plaquette term, one can 
interpolate between various actions and this allows to understand in more detail 
the properties of the phase structure, in particular how the strength of the transi- 
tion depends on the additional term and how this infiuences the approach to the 
continuum limit. 

That even a small value of bi can already have a large impact on the phase structure 
is illustrated in fig. 22 which shows the average plaquette value as a function of the 
hopping parameter k for three different actions, i.e. values of bi, namely bi = 
(Wilson), bi = -1/12 (tlSym) and fei = -1.4088 (DBW2). As one moves k from 
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Iv 

Fig. 22. Hysteresis of the average plaquette value as k is moved across the critical point, 
for Wilson, tlSym and DBW2 gauge action at a ~ 0.17 fm. 

the negative or positive side across the critical point, where the PCAC quark mass 
vanishes, a hysteresis in the average plaquette value develops whose size and width 
are indicators of the strength of the phase transition. Both the width and the size 
of the gap in the plaquette value decreases considerably as we switch on bi to bi = 
— 1/12 (tlSym action). Decreasing bi further down to bi = —1.4088 (DBW2 action) 
still seems to reduce the size of the gap, but the effect is surprisingly small despite 
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the large change in bi\^ The results in fig. 22 are for a lattice spacing a ~ 0.17 fm 
that is roughly the same for all three actions. I remark also that the strength of the 
first order phase transition weakens rapidly when the lattice spacing is made finer. 
This is illustrated in the left plot of fig. 20 for the case of the Wilson plaquette gauge 
action. 

A satisfactory setup for dynamical simulations with, say, Nf = 2 fiavours of quarks 
would be to reach pion masses of about 250 — 300 MeV and a box size of L > 2 fm. At 
the same time, one should stay at full twist to realize automatic 0(a) improvement. 
From all the results presented here it is conceivable that for the tlSym action this 
can be achieved with a reasonable computer time at a ~ 0.1 fm on L/a = 24 
lattices. Although for the DBW2 action the situation might be somewhat better, 
the advantages of the tlSym action, such as good convergence of perturbation theory 
and not big scaling violations as found in [125] for the DBW2 gauge action, suggests 
the tlSym gauge action as the action of choice. 

Of course other options could be advocated like sticking to the Wilson gauge action 
but simulating maybe at an even smaller lattice spacing, or modifying some parts 
of the fermionic action like the form of the covariant derivatives or by adding the 
clover term. In [71] it has been found that using the clover improved fermion action 
reduces the pion isospin splitting and thus w'. This result cannot be considered 
however conclusive for two reasons: the disconnected diagrams were neglected in 
the computation of the neutral pion; it was a quenched study. We have seen that 
these are critical issues because going from quenched to unquenched simulations 
and/or including the disconnected diagrams can have dramatic impacts in the phase 
diagram like even changing the sign of w', i.e. greatly affecting the neutral pion mass. 



^ In fact, the reductions are greater than appears, as the results with the Wilson action 
are for fj,q ^ 0, where the transition is weaker than for /_fq = 0, which is the value used for 
the other actions. 
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7 Renormalization and weak matrix elements 

In the previous sections we have extensively discussed the issue of the continuum 
limit and we have analyzed cutoff effects of order a and with Wtm. 

Renormalization is necessary in order to perform the continuum limit and correctly 
evaluate hadronic matrix elements. Here we will not discuss the way how the renor- 
malization is performed but only the mixing patterns of relevant physical quantities 
according to the lattice action used, i.e. according to the way the theory is regu- 
larized. Throughout the section we will assume that a mass independent renormal- 
ization scheme has been used, and that the scheme allows for a non-perturbative 
renormalization. Examples of such schemes are the so called RI-MOM scheme [126], 
and the Schrodinger functional (SF) [50,51]. 

Whatever is the chosen strategy, the renormalization is a difficult problem, that 
becomes even more difficult if one considers four-quark operators which appear in 
the effective weak Hamiltonian. The details of the renormalization patterns strongly 
depend on the symmetries preserved by the regularization adopted, i.e. on the lat- 
tice action. Regularizations which preserve chiral symmetry are clearly advanta- 
geous [127-129], and are increasingly used despite their high computational costs [63]. 
If computationally cheaper quarks of the Wilson type are used instead, the oper- 
ator renormalization is complicated, mainly due to the presence of power diver- 
gences [35,130] and/or mixings induced by the explicit chiral symmetry violation. It 
is then more than welcome to make use of the possibility to ease the renormalization 
pattern of relevant physical quantities using Wilson-like quarks. 

In this section I will give two examples of how Wtm can be used to reach this 
goal: the pseudoscalar decay constant and Bk- This result can be obtained using 
the relative freedom of choosing the twist angle depending on the flavour of the 
quark field appearing in the four-fermion operators. This freedom has been used 
in different ways in [39,98] leading to milder renormalization patterns for the four- 
fermion operators relevant for the A/ = 1/2 rule. The basic ideas on how to achieve 
these simplifications can be understood in an easier way by analyzing in detail the 
renormalization of Bk- This is presented in sect. 7.2. A detailed treatment of the 
A/ = 1/2 rule goes beyond the scope of this report. We thus refer to the original 
papers [39,98] and references therein for a detailed analysis on this particular process. 

7.1 Decay constants 

The determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements is 
of phenomenological and theoretical interest and to test the unitarity of the CKM 
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matrix, it is crucial to have a precise determination of the matrix element \Vus\ that 
parametrizes the coupling of the m — > s transition in the Standard Model (SM). 

Lattice QCD computations can address this issue determining in a very precise way 
the decay constants of the pion and the kaon. 

In fact the knowledge from first principles of fx/fn allows for an accurate determi- 
nation of |Kis|/|Kid| from the ratio of leptonic kaon decay rates [131]: 



2 N 2 



where the prefactor K on the r.h.s. takes into account radiative corrections. By 
combining the experimental result r(i^ — >■ /xi/^(7))/r(7r — > /ii/^(7)) with the Lattice 
QCD determination of fx/ fn and the experimental value of \Vud\ one obtains \Vus\- 

At the moment, the lattice errors dominate the uncertainty in the determination 
of \Vus\, thus it is important to have available a lattice action that can minimize 
possible sources of systematic uncertainty. 

With Wtm it is possible to determine the pseudoscalar decay constant without any 
computation of renormalization factors and without the knowledge of any improve- 
ment coefficient, still with only O(a^) cutoff effects. It is only necessary to correctly 
tune the untwisted bare quark mass to full twist as extensively detailed in sec. 4. 

Due to the exact flavour symmetry of massless Wilson quarks, the situation with the 
vector Ward identity is the same as with Ginsparg- Wilson fermions, i.e. the classical 
PCVC relation (2.19) holds exactly on the lattice 

d;{Kix)Om = -2f^^e'^'>{P'>ix)0m (7.2) 

(where 9* is the lattice backward derivative) with the point-split vector current 

~ 1 ( t"" 

(7.3) 

and the local pseudoscalar density. This implies that the multiplicative renormal- 
ization constants of composite flelds which belong to the same isospin multiplet 
must be identical in order to respect the vector Ward identities. An example is the 
renormalized pseudoscalar density which has the structure 

P^ = Zp(P'^ + 5"=^^). (7.4) 

The mixing with the identity operator appears because standard parity is not a 
symmetry of the massive lattice theory, but while P^'^ are still protected by the 



99 



symmetries Vp'^ defined in eqs. (2.38), is not. The vector Ward identity (7.2) 
here imphes that Zp is the same for all flavour components, and 

Zp = ^. (7.5) 

Moreover since standard parity combined with a change of the sign of the twisted 
mass V defined in eq. (2.39) is a symmetry of the lattice action, the mixing pattern 
of is actually only with a quadratically divergent term ^ similarly to what 
happens with overlap fermions. 

The existence of the point-split vector current can be used as a tool to determine 
the finite renormalization constant Zy that normalizes the local vector current. One 
possibility would be to determine Zy through 

with the local vector current defined in eq. (2.10). This method requires the com- 
putation of a correlation function involving the point-split current. This can be 
bypassed substituting with ZyV^ in eq. (7.2) and the relation is now valid up to 
order lattice artifacts. This relation allows for a determination of Zy through 

E.{a^vy(i)p'(o)> 

where 9^ is the symmetric lattice derivative. 

In fig. 23 there is an example of Zy as a function of the quark mass for a quenched 
simulation [66] determined using eq. (7.7). It turns out that Zy is to a good approx- 
imation linear in (afiq)'^. 

In contrast, the axial Ward identity does not hold in the bare theory. Axial Ward 
identities therefore provide normalization conditions which determine finite renor- 
malization constants such as Z^, or finite ratios of scale dependent renormalization 
constants, such as Zs/Zp [35]. 



As Wtm and standard Wilson quarks are not related by a lattice symmetry, the 
counterterm structure for composite fields with the same physical interpretation 
depends upon the twist angle u, and a particular choice of ui can lead to substantial 
simplifications. A typical example is the computation of the pseudoscalar decay 
constant from the axial coupling to the pion. We have argued in sec. 2 that for 
renormalized correlation functions, the relation between physical and twisted basis 
can be inferred from the corresponding classical relations if the renormalization 
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scheme adopted is mass independent. With Wilson fermions usually the pseudoscalar 
decay constant is extracted from the correlation function 

«'E((^r)S(^)^r(0)). (7.8) 

X 

This requires first the determination of the renormalized axial current and eventually 
all the relevant improvement coefficients. Performing the axial rotation (2.6) we 
obtain that in the continuum limit we should have (we set the isospin component 
a = 1) 

{{AK)l{x)Pl{y))iM^,0)=COs{u){{A^)l{x)P^{y))(^:rnn^^^^ 

(7.9) 

If we set = 7r/2 the r.h.s of eq. (7.9) only contains the vector current which is 
protected against renormalization. Moreover using the exact PCVC relation on the 
lattice (7.2) it is easy to show [53,85,88] that 

/ps- ^^1^ a- 1,2, (7.10) 

where Mps is the charged pseudoscalar mass and |PS) denotes the corresponding 
pseudoscalar state. As we see, the neat result is that Wtm at full twist allows the de- 
termination of an automatically 0(a) improved pseudoscalar decay constant without 
the knowledge of any renormalization constant and any improvement coefficient. 
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In fig. 24, I summarize Wtm results coming from quenclied [66,95] {Nf = 0) and un- 
quenclied [69] {Nf = 2) simulations, for the pseudoscalar decay constant over a range 
of pseudoscalar masses around [300 — 700] MeV. We recall that for Wtm, to produce 
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Fig. 24. Summary of the Wtm numerical results at io = ■k/2 for the pseudoscalar decay 
constant /ps as a function of Mpg. The bending phenomenon is visible with the "Wilson 
pion" definition of the critical mass for Nf = 0. The corrected data for Nf = 2 (slightly 
displayed for clarity) include the effect of the non vanishing mpcAC (^-H). The perfect 
consistency of the data dispels all the doubts about the possible origins of the mass de- 
pendence for /ps. 

such a plot, no renormahzation constant are needed. Even if the lattice spacing for 
the quenched (a = 0.093 fm) and the unquenched (a = 0.087 fm) simulations are 
not exactly the same, it is interesting to notice the effect of the dynamical quarks 
in the mass dependence of the pseudoscalar decay constant. Moreover 1 also plot for 
the dynamical case the pseudoscalar decay constant corrected using eq. (4.83) and 
using the PCAC quark mass determined in the same set of simulations 



^corr 
/PS 



/] 



PS 



(7.11) 



In the statistical errors the two results are perfectly consistent showing that the 
mass dependence of fps for iV/ = 2 is a genuine /Xq quark mass dependence and 
and the curvature is not artificially induced by a non exactly zero PCAC mass. 
As a comparison 1 show also the Nf = results obtained with the "Wilson pion" 
definition of the critical mass, where, on the contrary, the quark mass dependence 
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Fig. 25. Comparison of the continuum limit (filled symbols) at roMps = 1.2544 obtained 
in Wtm at full twist without [66] (red squares) and with [72] (black circles) improvement 
coefficients. 

is largely given by the cutoff effects, as discussed in sect. 4. 

In ref. [66,72] the continuum limit of the quenched data is performed at = 7r/2 
producing consistent results. In fig. 25 we see an example of the two continuum limits 
for a pseudoscalar mass fixed at the kaon mass (roMps = 1.2544). While in ref. [66] 
no improvement coefficients were used, in ref. [72] the whole set of improvement 
coefficients were used (csw, ca,--.) in the computation of the decay constant from 
various lattice definitions. Some of the improvement coefficients are known only at 
one loop, but obviously at full twist they are actually irrelevant. The data shown in 
fig. 25 corresponds to a definition involving only csw and by (this is the coefficient 
multiplying the 0(a) counterterm in eq. (D.15)). From the plot it is evident how 
automatic 0(a) improvement makes the usage of the improvement coefficients not 
relevant for the cancellation of the 0(a) discretization errors. In ref. [66] the contin- 
uum limit is performed over a wide range of quark masses as discussed in sect. 4.5 of 
this review (e.g. see fig. 7). In ref. [72] several definitions of the decay constant allow 
to perform a constrained continuum limit with a final result which has a relative 
error smaller than 2%. 

The very precise unquenched data down to pseudoscalar masses around 300 MeV [69] 
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allow a discussion of whether xPT formulae can reproduce the mass dependence for 
amps and a/ps- One possible source of systematic uncertainty is the finite size ef- 
fects, and this can be taken into account using xPT. In particular the lowest two 
quark masses turn out to be significantly affected. Preliminary results at a smaller 
lattice spacing [90,91] show small discretization errors, indicating that discretization 
errors are under control. Therefore one could try to use continuum xPT to describe 
consistently the dependence of the data both on the finite spatial size (L) and on 
the quark mass /iq. 

The fit to the raw data for Mps and /ps has been performed simultaneously with 
the appropriate {Nf = 2) xPT formulae [132, 133] 



l+eiog(25o^q/A^) , (7.12) 



/ps(L) = F[l-2e^i(A)] l-2eiog(25o/iq/Ai) , (7.13) 



where 



e = 2So/iq/(47rF)2 , A = ^2Bofic,L^ . (7.14) 
The finite size correction function ^i(A) was first computed by Gasser and Leutwyler 
in ref. [132] and is also discussed in ref. [133]. In eqs. (7.12) and (7.13) NNLO 
xPT corrections are assumed to be negligible. The formulae above depend on four 
unknown parameters, Bq, F, A3 and A4, which will be determined by the fit. For 
details about the data analysis I refer to the original paper [69]. 

For the lightest four values of a/Xq, an excellent fit to the data on /ps and Mps is 
found (see figures 26 and 27). The fitted values of the four parameters are 

2aBo = -4.99(6) , 

aF = -0.0534(6) , 
log(a2A2) = -1.93(10), 

\og{a^Al) = -1.06(4) . (7.15) 
The data are clearly sensitive to A3 as visualized in figure 26(a). 

The value of afiq, afi^^, at which the pion assumes its physical mass, is deter- 



mined [69] requiring that the ratio [y[Mpg{L = 00)]/ fps{L = 00)] takes the value 
(139.6/130.7) = 1.068. From the knowledge of a/i^ one can evaluate 

[3,4 = \0g{AlJM^s)\Mps=m^ 

afi^ = 0.00078(2), [3 = 3.65(12), k = 4.52(06). (7.16) 

These results compare nicely with other determinations (for a review see ref. [134]). 
The inclusion of the results from OyUq = 0.0150 in the fit gives an acceptable descrip- 
tion of Mpg but misses the data for /ps, as shown in figures 26(b) and 27. Note, 
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Fig. 26. In (a) I plot (aMps)V(a/Uq) 

as a function of a//q. In addition the fit with 
eq. (7.12) to the data from the lowest four values of /iq is compared with linear and 
quadratic fits. In (b) I plot (aMps)^ as a function of a//q with the corresponding chiral fit. 
In both figures (a) and (b) the raw and the finite size corrected (L — > oo) data are plotted. 



however, that in Eqs. (7.12, 7.13), and thus in the fit results (7.15, 7.16), a number 
of systematic errors, as discussed below, are not included. 

The values presented here should hence be taken as a first estimate, the validity of 
which has to be checked in the future. Nevertheless, the statistical accuracy achieved 
implies that there is a very good prospect of obtaining accurate and reliable values 
for the low-energy constants from Wtm fermion simulations. 

Based on the physical value of fj,, one gets 

a = 0.087(1) fm. (7.17) 

Using the value of ro/a reported in ref. [69], this lattice calibration method yields 
ro = 0.454(7) fm. 

We now discuss the possible sources of systematic error. This analysis is based on 
lattice determinations of properties of pseudoscalar mesons with masses in the range 
300 to 500 MeV on lattices with a spatial size slightly above 2 fm. Systematic errors 
can arise from several sources: 

(i) Finite lattice spacing effects. Preliminary results at a smaller value of the lattice 
spacing that were presented in refs. [90,91] suggest that 0{a) improvement is nicely 
at work and that residual O(a^) effects are small. 

(ii) Finite size effects (FSE). Even if the results presented here are obtained from a 
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Fig. 27. Plot of a/ps (before and after finite size corrections) as a function of a^Uq together 
with a fit to formula (7.13). This fit is compared with a quadratic fit including all the 
data points and excluding the heaviest mass. A comparison with the impact of different 
physical values for tq is also shown. 



box of physical size L >2 fm, to check that next to leading order (continuum) xPT 
adequately describes the FSE, a run on a larger lattice is required. 

(iii) Mass difference of charged and neutral pseudoscalar meson. In the appropriate 
WxPT power counting for our values of the lattice spacing and quark masses, i.e. 
/^R ~ ~ (4.67), the pion mass splitting is a NLO effect (M^ri)^ — (M^o)"^ = 
0{a^AQQY))j from which it follows that to the order we have been working the effects 
of the pion mass splitting do not affect, in particular, the finite size correction 
factors for Mps and fps- In spite of these formal remarks, it is possible, however, 
that the fact that the neutral pion is lighter than the charged one (by about 20% 
at a/iq = 0.0040, see ref. [69]) makes inadequate the continuum xPT description of 
finite size effects adopted in the present analysis. This caveat represents a further 
motivation for simulations on larger lattices, which will eventually resolve the issue. 

(iv) Extrapolation to physical quark masses. We are assuming that xPT at next 
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to leading order for the Nf = 2 case is appropriate to describe the quark mass 
dependence of Mpg and fps up to ~ 450-500 MeV. The lattice data are consistent 
with this, but it would be useful to include higher order terms in the xPT fits as 
well as more values of a/iq to check this assumption. The effect of heavier quarks in 
the sea should also be explored, and preparatory studies with Nf = 2 + 1 + 1 [41] 
dynamical quarks show that the inclusion of heavier quarks in dynamical simulations 
is accessible using current algorithms and machines. 



7.2 Bk 



Indirect CP violation in K ^ titt decays is measured by the parameter Ek-, defined 
in terms of kaon decay amplitudes as 

where / is the total isospin of the two-pion state. The long distance non-perturbative 
QCD contribution to \sk\ is provided by the following matrix element 

{K^\Or-'\K^) ..... 

where O^'^^^ is the effective four-quark interaction renormalized operator, with bare 
operator 

0^'=' = {-siidTs^id) , (7.20) 
where 7^ = 7^(1—75). Bk largely dominates the uncertainty on the standard model 
(SM) value for {ekI and improving the accuracy of Bk is essential in order to derive 
stringent bounds on the amount of non-SM CP violation in kaon decay. 

One of the important sources of uncertainty in lattice QCD computations of Bk 
with Wilson fermions arises from operator renormalization. The operator O^^^^ is 
usually split into parity-even and parity-odd parts as 

qAS-2 ^ Oyy+AA — OvA+AV, (7-21) 

where with obvious notation we have 

Ovv+AA = {s'y^d){s'j^d) + ( 37^,75 rf) (37^750?) , (7.22) 

OvA+Av = is-ff,d){s'jf,j5d) + {s'jf,-f5d){s'y^d) . (7.23) 
Since parity is a QCD symmetry, the only contribution to the K^-K^ matrix ele- 
ment comes from Ovv+aa- In regularizations which respect chiral symmetry, the lat- 
ter operator is multiplicatively renormalizable. If chiral symmetry is not preserved, 
Ovv+AA mixes with four other dimension-6 operators [135-137] with positive parity: 
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(Or)vv+aa(/u) = ^vv+aa(5'o, O/W) 



Ovv+aa(5'o) + J2^ii9o)Oi{go) 



i=l 



(7.24) 



The operators Oi{go) belong to different chiral representations than Ow+aa- The 
mixing coefficients Ai{gQ) are finite functions of the bare coupling, while the renor- 
malization constant Zw+aa diverges logarithmically in a/i, where fi is the renor- 
malization scale. 

There have been several proposals to eliminate this operator mixing with Wilson 
fermions, all based on the observation [137, 138] that, even in the absence of chiral 
symmetry, the operator Ova+av is protected from finite operator mixing by discrete 
symmetries, and thus it renormalizes multiplicatively 

(C'R)vA+Av(/i) = ^VA+Av(fi'0, a/^)C'vA+Av(fi'o) • (7.25) 



In order to show this, it is convenient to work with massless Wilson fermions, also 
having in mind a mass independent renormalization scheme. 

First of all we note that O'^^^^ cannot mix with operators of lower dimensionality, 
because such operators do not have the four-flavour content of the original one. 
Thus O^^^^ can mix with any other dimension-six operator, provided it has the 
same quantum numbers, i.e. with any operator which has the symmetries of O'^^^^ 
and of the lattice action. The generic QCD Wilson lattice action with 4 massless 
quarks is symmetric under parity V defined in eq. (2.34), and charge conjugation C 
(see app B for the definition). We define the generic four fermion operators 

Ora)r(2) = (xir(^)x2)(X3r('^X4) (7.26) 

for all r*^^) and T^"^^ combinations of interest. Moreover, there are other useful 
(flavour) symmetries of the action, namely the switching symmetries S' and S" 
defined by [138] 

S': ^^2, ^^27) 
(X3 ^ Xi 

S": ^^4' (7.28) 

(X2 ^ X3- 

In Table 1 I classify the operators Or(i)r{2), for all T^^^ and T^'^^ with negative parity, 
according to the discrete symmetries C, S' and S" , with notation 

Or(i)r(2)_|_r(2)r(i) = Or(i)r(2) ±0r(2)r(i)- (7.29) 

Parity violating operators, for which CS" is not a symmetry, have been symmetrized 
or antisymmetrized in order to obtain eigenstates of CS" . 
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Or(i)r(2) 


V 


CS' 


cs" 


CVS' 


CVS" 


A+AV 


— 1 


— 1 


— 1 


+1 


+ 1 


( Jt / A AT / 

A— Av 


-1 


-1 


+ 1 


+1 


—1 


Osp~ps 


— I 




— I 


— I 




Osp+ps 


-1 


+1 


+1 


-1 


-1 




-1 


+1 


+1 


-1 


-1 



Table 1 

Classification of parity violating four-fermion operators Op{i)p(2) according to useful prod- 
ucts of discrete symmetries V, C, S' and S" . Note that Off = Ott a.nd Orpf = Ofj.. 

The parity violating four- fermion operators listed in table 1 do not all have identical 
CVS' and CVS" values. It is straightforward to see that the operator Ova+av cannot 
mix with the other parity violating operators . 

We are going now to discuss two possible approaches [98, 139] which use Wtm to 
extract the renormalized (-S'°|Ovv+aa|-^'^°) by relating it to the matrix element of 
a parity violating operator. As the Wtm action differs from the standard Wilson 
fermion action by the mass term, the renormalization properties of composite oper- 
ators in mass independent renormalization schemes are not modified. In particular, 
OvA+Av remains multiplicatively renormalizable, with the same renormalization con- 
stant and running as with Wilson fermions. Thus finite subtractions are avoided in 
the Wtm determination of Bk- In order to disentangle the two approaches we will 
call the proposal made in [70, 139] the "mixed twist" method while the proposal 
in [98] the "mixed action" method. Hopefully the reasons for these names will be- 
come clear in the following. 



7.2.1 "Mixed twist" method 

The first variant of the "mixed twist" method consists in choosing the following 
fermionic action 

S'^''^^ = [x(a^)Pw,sw + mi+ ilbT'^ l^i)x{x) + s(x)(L'w,sw + mjs(x)]. (7.30) 

X 

where D^^^w is the clover improved Wilson operator (4.18,4.19) and x collects the 
light doublet 



X 



u 
d 



"^^ It could still mix with the Fierz transformed in Dirac space operator Oyyi_|_^y It turns 
out that Oy^_^_^y = OyA+AV (see ref. [137] for details). 
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The label tt/2 refers to the choice of the twist angle that can be obtained setting 
"^R,/ = 0. We see that while the twist angle of the first doublet is 7r/2, the twist 
angle of the strange quark is zero, hence the name "mixed twist". 

The second variant is based on the action 

S'^'^^ = a^Y^ [m(x)(L'w,sw + mu)u{x) + x(a;)(Dw,sw + mi + i'^^r^ ^ii)x{x)] . (7.31) 

X 

where now x collects a doublet made of a down and a strange quark 




To set now the twist angle to ir/A requires /xr^; = mpj,^;. The variant with action (7.31) 
assumes a priori that the s and d quarks have degenerate physical masses. This is 
necessary if one would like to extend this variant to unquenched simulations: a non 
degenerate diagonal twisted doublet would lead to a complex determinant [39]. As 
long as this action is used in quenched QCD this restriction is not needed but all the 
computations carried out with it are performed with degenerate s and d quarks. The 
action in eq. (7.30), on the other hand, is perfectly well suited for an unquenched 
computation. 

It is important to stress that in both the variants not all the quark flavours are fully 
twisted, i.e. automatic 0(a) improvement [11] does not apply. To have full 0(a) 
improvement of the matrix element it would be necessary to subtract a number of 
dimension-seven counterterms from the four-fermion operator. Such a procedure is 
highly impractical, and has not been pursued. Hence leading cutoff effects in Bk 
with the "mixed twist" methods are expected to be linear in a. 

We have already discussed in sect. 2, that to relate renormalized correlation functions 
computed with tmQCD to QCD, it is enough to perform the needed change of 
variables (axial rotations). This explains the choice of the twist angles for the two 
variants. In fact performing the rotation in eq. (2.6) for the two actions, we obtain 
in both cases 

{K^\ (0^)vV+AA |i?°)QCD = (0«)vA+AV |i?°)tmQCD , (7.32) 

which holds in the continuum limit for the two versions of Wtm under considera- 
tion. From this identity, Bk can be extracted from a K^-K^ matrix element of the 
multiplicatively renormalizable operator Ova+av- 

Both these variants has been used to compute Bk with quenched fermions. 

In particular the non-perturbative renormalization has been performed in the SF 
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scheme [140], and the matrix element has been computed using both the variants 
in quenched simulations [70,72]. The twist angle has been tuned to u = 7i/2 setting 
'^R,/ = in the clover improved theory. This corresponds to the full twist tuning 
discussed in sect. 4.3.1. To tune the twist angle to 7r/4, the untwisted quark mass 
has been tuned such that /iR,; = m^j, which via eqs. (4.22,4.23) and eqs. (4.25,4.26) 
translates into 



anig^i = ^ani{l + ^(6^ - &r«) - a/^/} (7.33) 

For a given choice of a/i;, k = (2am/ + 8)~^ is tuned so that am^ / satisfies one 
of the two above relations. This requires the knowledge of all the renormalization 
constants and improvement coefficients. The final result for the renormalization 
group invariant matrix element Bk being 

Bk = 0.735(71). (7.34) 

The only systematic uncertainty that affects this computation is the usage of quenched 
fermions. 

In fig. 28 we summarize results for B^ obtained by several collaborations using 
different lattice actions. In particular we observe that Wtm allows a very precise 
determinations of Bk which is competitive with the determination obtained with 
other discretizations. The main drawback of the "mixed twist" approach is the fact 
that the leading cutoff effects in B^ are expected to be hnear in a. 



7.2.2 "Mixed action" method 

The second method we are going to discuss has all the quark flavours at full twist, 
thus retaining automatic 0(a) improvement, but it makes use of different lattice 
actions for valence and sea quarks. In particular the valence actions are chosen in 
order to remove the unwanted mixings for Ova+av- To obtain this goal the number 
of valence quarks has to be extended, i.e. some of the valence quarks involved in the 
correlation function under investigation have to be doubled with a slightly different 
lattice action. The resulting lattice theory is not unitary thus the approach to the 
continuum limit will be affected by O(a^) effects coming from unitarity violations. 
With the "mixed action" a detailed understanding of these particular discretization 
errors have to be achieved in order to correctly perform the continuum limit. This 
does not pose any problem of principle because once the renormalized quark masses, 
or equivalently the corresponding pseudoscalar meson masses, are matched in the 
continuum limit, one obtains the desired and correct renormalized matrix elemental 

This statement could be put at a more formal level using a Ginsparg- Wilson regular- 
ization. 
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Fig. 28. Summary plot of the results for Bk obtained by several collaborations using 
different lattice actions. Prom the top to the bottom the refs. are [72,141-153]. The Nj = 
Wtm result [72] is shown with a filled square. 



For the light (/) and the heavy {h) sea doublet quarks the actions are 



^(0 



.4 Y^: 



(7.35) 



Swim = Xhi^) Dw + mo + i^ihlbT^ + e/,r^ 



Xh{x) 



(7.36) 



For the particular case of we introduce the following valence quarks: s and 
s' . The d and the s valence quarks have the Osterwalder-Seiler (OS) action 



qif) _ .4 V- 



■'OS 



D^f^ + mo + 2/i}75 



Xf{x) 



(7.37) 



In principle it is possible to choose also a non-degenerate quark action for the light 
sector. 
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with f = d,s while the s' quark has the same OS action 



s (x) 



(7.38) 



but with an opposite sign for the mass term. While this certainly does not change 
the sign of the physical quark mass, it changes the leading discretization errors of the 
actions, in such a way that to extract Bk it is possible to use a specific correlation 
function that renormalizes multiplicatively and it is automatically 0(a) improved. 

The OS action differs from Wtm by the fact that it does not violate isospin. This 
action is not used for dynamical fermions since, as we have discussed in sect. 2, it 
could generate an FF in the renormalization process through vacuum polarization 
diagrams. Because of the peculiar flavour structure, this does not happen for Wtm, 
but anyhow the OS action can be still used for valence quarks. To correctly define 
the theory we would need to introduce for each valence quark, the corresponding 
ghost action which cancels exactly the contribution to the determinant [154]. We 
will assume that this it has been done. All the considerations which follow will not 
involve correlation functions of ghost fields, thus the standard Symanzik expansion 
will apply [155], and all the symmetries we are going to use will be naturally extended 
to the ghost fields. 

We reiterate that the valence and sea quark masses have to be correctly matched 
to obtain the correct continuum limit. We also recall that to obtain automatic 0(a) 
improvement the untwisted quark mass has to be tuned to the critical value, and so 
the physical quark mass is totally given by the twisted quark mass. This is obtained 
with the following matching prescription 

(/^d)R = (/^d)R, 

(/i.)R = (/i:)R = (/i^OR- (7.39) 



With the choice of the valence quark actions specified in eqs. (7.37,7.38), it is easy 
to show, performing the usual axial rotations, that the original matrix element 

{lf\0^'='\K%CB = \MlflBK (7.40) 
is equivalent to the matrix element 

{k"'\2{0^)va+av\K'') = ^MKMK'fxfK'BK (7.41) 

where this second matrix element on the lattice is computed with the model specified 
before. In particular the bare operator reads 

OvA+Av = i [(s7mC?)(s'7m75C?) + {87^,15(1) (s'-f^d)] (7.42) 
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and the interpolating fields for the kaons are 



K\x) = td{x)s{x), k'\x) = s (x)75rf(x). (7.43) 

The matrix element (7.42) can then be extracted by the asymptotic behaviour of 
the correlation function 

CK'OK{xo,yo) = 2a''Y.(idi5s'){x)OvA+Avmds)iy)), (7.44) 

with Ho < and xq > 0. The values of Mk.Mk'J'k and fx', can be extracted in a 
standard fashion from the 2-point correlation functions involving the interpolating 
fields in (7.43). We repeat here briefly the reason why Ova+av is multiplicatively 
renormalized. Mixing with lower dimensional operators is excluded by flavour sym- 
metry of the massless lattice actions, in particular Ova+av can mix only with op- 
erators such that As = As' = — A(i/2 = 1. Parity symmetry of the action in the 
chiral limit forbids mixing with operators of the opposite parity. Possible dimension 
6 operators of negative parity that could mix with Ova+av are 

OvA-Av = {s-ff,d){s'-f^-f5d) - {s-f^-f5d){s'-ff,d) (7.45) 

OsP±Ps = (sd) (s'^^d) ± (s^^d) (s'd) (7.46) 

Oj^f = e^uXp{scr^i,d){scrxpd) (7.47) 

To rule out these operators is enough to use CVS' and CVS" symmetries defined in 
eqs. (7.27,7.28), as we have done previously. 

The main advantages of this approach is that all the quarks are at full twist so to 
prove automatic 0(a) improvement it is enough to generalize the form of the twisted 
parity used in sect. 4 

f/(xo, x; 0) — > t/(xo, — x; 0), t/(xo, x; /c) — f/^^(xo, — x — a/c; /c), = 1, 2, 3 
X/,h(xo,x) loiil5T^)xi,hi^o, -x) 
: { Xi,hixo, x) ^ Xi,hixo, -x) {il5T^ho 
g/(xo,x) ^ 7o(275)g/(xo, -x) 
^g/(xo,x) qf{xo, -x)(z75)7o, 

(7.48) 

where qf generically indicates all the valence fields. Moreover non-degenerate quarks 
can be introduced retaining the nice property of a real and positive definite deter- 
minant [38]. This approach requires care in the matching procedure between valence 
and sea quarks. 

An obvious alternative to avoid renormalization problems consists in using regu- 
larizations with exact chiral symmetry. However, the computational costs involved 
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make it difficult to perform continuum limit extrapolations and study finite vol- 
ume effects. A ffist attempt using Wtm sea and overlap valence quarks has been 
presented at the last lattice conference [156]. 

We reiterate that one drawback of the mixed action approach is that the resulting 
theory is not unitary at non-zero lattice spacing. This introduces peculiar O(a^) 
cutoff effects, that nevertheless can be described by mixed action chiral perturbation 
theory [157,158]. This cutoff effects have been investigated thoroughly in refs. [159, 
160], and formulae, valid up to O(a^), have been given for correlators computed with 
valence overlap fermions on a set of different sea quark actions. 

The preliminary results of ref. [156] indicate that the whole procedure is limited by 
the statistical accuracy achievable on a large volume with overlap fermions, in order 
to correctly perform the matching procedure. A possible way to solve this problem 
is to match the current quark masses in a finite volume. 
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8 Algorithms for dynamical fermions 



At present the only practical way to simulate numerically a 4-D euclidean quantum 
field theory with fermions is to perform the Grassmann integral on the fermion 
fields analytically and then to apply Monte Carlo methods in the resulting effective 
bosonic theory. 

After integrating out the fermion fields x^X the partition function (2.22) of Wtm 
for Nf = 2 degenerate flavors is given by 

Zoc J VUdet{Q^Q)e-^'^^^^ = J VUe-^^«^^^ , (8.1) 

with 

S,s[U] = SG[U]-TT\ogQ^Q . (8.2) 

The reality of the effective action is only guaranteed by the positivity of the quark 
determinant. 

8.1 Hybrid Monte Carlo 

The determinant det((5^(5) can be expressed in terms of the so called pseudofermion 
complex flelds 

det{Q^Q) oc J V(j)'^V(j)exp(^-SpF[U,(f)\(f)]^ , (8.3) 

where Spf[U, (p] = {Q~^(f), Q^^cfy is the pseudofermion action with standard scalar 
product (u, v). The pseudofermion flelds are formally identical to the fermion flelds 
X, but follow the statistics of bosonic flelds. 

A way to obtain a "good" global update with the non-local effective bosonic action 
S'pF is to simulate a microcanonical ensemble with a suitably deflned Hamiltonian. 
The Hamilton equations have to be integrated with a suitable integration scheme, 
and it is then possible to correct for the the flnite step errors with a flnal stochastic 
accept-reject step. This is schematically what is called the Hybrid Monte Carlo 
(HMC) algorithm [161]. To be more speciflc, the (f) version of the HMC algorithm is 
based on the Hamiltonian 

H{Ii, U, 0, 0^) = ^ E Tr [n(x, fxf] + Sg[U] + Spp[U, 0, 0t] , (8.4) 

where we introduced traceless hermitian momenta Il{x,fi). The HMC algorithm is 
then composed by the following steps 
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• Global heat-bath for momenta and pseudofermion fields: 

the initial momenta are randomly chosen according to a Gaussian distribution 
exp(— n^/2), and the random fields R are produced from a distribution like 
exp(— i?'''i?)with the initial pseudofermions computed as (p = Q ■ R. 

• Molecular dynamics evolution: 

propose a gauge configuration U' and a momentum 11' integrating the Hamilton 
equations of motion for the gauge field U and the momentum 11 at fixed pseud- 
ofermion field (f). 

• Metropolis accept /reject step: 

The proposals U' and 11' are accepted with probability min{l, exp(— Ai^^)}, where 
AH = if(n', U', 0, (j)^) —H(Il, U, 0, 0'''). This step is needed because of the numer- 
ical inexact integration of the equations of motion. 

It possible to prove that the HMC algorithm satisfies the detailed balance condi- 
tion [161] and hence the configurations generated with this algorithm correctly rep- 
resent the intended ensemble. Since the Hamiltonian is conserved up to finite step 
errors, the integration can be set up such that the gauge configurations are globally 
updated keeping the acceptance rate high. 

8.2 Molecular dynamics evolution 

In the molecular dynamics part of the HMC algorithm the gauge fields U and the 
momenta H need to be evolved in a fictitious computer time r. Differentiation on a 
compact Lie group is defined, given a generic function G, by 

G'[[/e-«T«] ^ ^ iUaS:G[U] + 0{iu') , (8.5) 

where U is group element, u is an infinitesimal vector with dimension equal to that 
of the adjoint representation and are the group generators. Since the variation 
6^G takes values in the Lie algebra it is natural to let the momenta H(x, /i) be an 
element of the Lie algebra. With respect to r, Hamilton's equations of motion read 

^U{x, fi) = H(x, i2)Uix, /i) , -^H(x, /i) = -F{x, /i) , (8.6) 

where the force F{x, fi) is obtained by differentiation with respect to the gauge links 

{uj,F) = Su,S, S^U{x, jj.) = uj{x, jj,)U{x, jj.) (8.7) 

and S = Sq + Spp. Since analytical integration of the former equations of mo- 
tion is normally not possible, these equations must in general be integrated with a 
discretized integration scheme that is area preserving and reversible. The discrete 
update with integration step size Ar of the gauge field and the momenta can be 
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defined as 

Tu(Ar): U ^ U' = E {zAtU{x, fi)) U , 

Ts(Ar) : n ^ n' = n - iAtF , ^ ' ' 

where E{-) stands for tlie SU{3) exponential function. To integrate tlie equation of 
motion (8.8) it is possible to use many different integrators. One of the simplest is 
the so called leap-frog integrator. Given eq. (8.8) one basic time evolution step of 
the leap-frog reads 

T = Ts(Ar/2) Tu(Ar) rs(Ar/2) , (8.9) 
and a whole trajectory of length r is achieved by a number of molecular dynamics 
Nmd = t/At successive applications of the transformation T. 

8. 3 Preconditioning 

For each step along a trajectory of length At the force F has to be computed. 
The computation of the force F is the most expensive part in the HMC algorithm 
since the inversion of the Wilson-Dirac operator is needed. In order to improve the 
integration (i.e. to increase the step size At or equivalently to reduce the number of 
steps A^md) it would be better to have small forces along the trajectory. In general 
this is achieved by preconditioning the HMC algorithm. A good strategy would then 
be to split the forces in such a way that those which require more computer time 
to be computed are those with the smallest magnitude. When the chiral limit is 
approached, the quark forces tend to increase, which requires the step sizes to be 
adjusted accordingly [162], and the choice of the preconditioner can have an influence 
on this behaviour. A parameter that cannot be predicted in a dynamical simulation 
is the autocorrelation time of the quantity one is interested in. Extended simulations 
are needed in order to understand if the autocorrelation time is under control or not. 

In general, preconditioning is always associated with a factorization of the quark de- 
terminant into the determinants of certain different operators. The number and the 
type of operators depend on the chosen preconditioner. A possible efficient precon- 
ditioning is obtained using a domain decomposition [163] (DD). The preconditioner 
we are going to discuss is the so called Hasenbusch acceleration or mass precondi- 
tioning [164]. From now on I will discuss only the case of the plain hermitian Wilson 
operator Qw, because all the conclusions can be easily extended to the Wtm oper- 
ator. It was realized in ref. [164] that using the identity 

det Ql = det (qI^ + p") det i J ' (^-^O) 

with an adjustable mass shift p, can speed up the HMC algorithm. Each of the two 
determinants on the r.h.s. of eq. (8.10) is treated by a separate pseudofermion field 
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0j and a corresponding pseudofermion action SpPi ■ The Hamiltonian can be written 
as 

H=Iy. Tr[n(x, + Sg + SpF, + ^PF, . (8.11) 

where 5'pfi and Spp^ are the pseudofermion actions respectively related to det (Qvi/ + P 
and det 



2^ 



In ref. [165,166] it was argued that the optimal choice for p is given by = y/XuKn- 
Here Am (Am) is the maximal (minimal) eigenvalue of Q^. The reason for that choice 



is obtained from minimizing the sum of the condition numbers K = Ki + K j ^1 for 



the operators appearing in Spp-^ and Spp^ ■ With the optimal pl^^ = y/XuKn the two 
condition numbers Ki and K2 are equal to ^Xu/^m, both of them being smaller 
than the condition number of Q^r which is Xu/Xm- 

Since the force contribution in the molecular dynamics evolution is supposed to be 
proportional to some power of the condition number, the force contribution from 
the pseudofermion part in the action is reduced and therefore the step size At can 
be increased, in practice by about a factor of 2 [164, 165]. 

This preconditioning can be very easily adapted to the Wtm operator observing that 
Qw + = Q^Q with p being the twisted mass and Q being 75 "times" the Wtm 
operators (2.66). 



8.4 HMC with multiple time scale integration and mass preconditioning 



In the paper by Liischer [163] impressive acceleration factors were obtained with a 
DD preconditioned HMC compared with a plain HMC. One of the reasons for this 
impressive result is the observation that the resulting forces after the precondition- 
ing have a magnitude which decreases as the computer time needed to invert the 
corresponding operator increases. It is then beneficial to integrate each single force 
with different time steps. This hierarchy of forces is obtained with a preconditioner 
that provides a strong infrared cutoff and separates low- and high-frequency modes 
of the system. In [167] the idea was explored of using as an infrared cutoff the mass 
in the Hasenbusch acceleration method, which could also give a similar hierarchy of 
forces that could then be combined with a suitable integrator. Therefore it might be 
advantageous to change the point of view: instead of tuning p a la refs. [164,165], i.e. 
minimizing the condition number of the operators appearing in the pseudofermion 
action, rather to exploit the possibility of arranging the forces by the help of mass 

The condition number is the ratio of the maximal over the minimal eigenvalue of a 
given operator. 
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preconditioning with the aim to reach a situation in which a multiple time scale 
integration scheme is favorable, i.e. tuning p to achieve a hierarchy of forces as with 
the DD-HMC. 

The idea to combine a multiple time step integrator with a separation of infrared 
and ultraviolet modes was already proposed in ref. [168]. This idea was applied to 
mass preconditioning by using only two time scales in refs. [169,170] in the context 
of clover improved Wilson fermions. However, a comparison of results presented 
in the next section to the ones of refs. [169, 170] is not possible, because volume, 
lattice spacing and masses are different. The gain reported in [169] compared to the 
performance of the method described in [164, 165] was at most 20%. 

In order to generalize the leap frog integration scheme (8.9) we assume, in the 
following, that we can bring the Hamiltonian to the form 

1 " 
H = -Y.TT[Ii{x,pf] + Y.Su{U], (8.12) 

with n > 1. For instance with n = 1 5*0 might be identified with the gauge action 
and 5*1 with the pseudofermion action of eq. (8.4). 

Clearly, in order to keep the discretization errors small in an algorithm like leap 
frog, the time steps have to be small if the driving forces are large. Hence multiple 
time scale integration is a valuable tool, if the forces originating from the single 
parts in the Hamiltonian (8.12) differ significantly in their absolute values. Then 
the different parts in the Hamiltonian might be integrated on time scales inversely 
proportional to the corresponding forces. 

The leap frog integration scheme can be generalized to multiple time scales as has 
been proposed in ref. [171] without loss of reversibility and the area preserving prop- 
erty. The scheme with only one time scale can be recursively extended by starting 
with the definition 

To = Ts„(Aro/2) ru(Aro) Ts,{Aro/2) , (8.13) 

with Tu defined as in eq. (8.8) and where Ts^(Ar) is given by 

Tg,(Arfc) : P ^ P-iAnFk. (8.14) 

As Atq will be the smallest time scale, we can recursively define the basic update 
steps Tfc, with time scales Ar^ as 

T, = Ts,(Ar,/2) [T,_i]^-^ Ts,(Ar,/2) , (8.15) 

with integers A^^ and < k < n. One full trajectory r is then composed by [T„]^". 
The different time scales Arfc in eq. (8.15) must be chosen such that the total number 
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of steps on the k-th time scale A^md^, times Ar^ is equal to the trajectory length r 
for all < k < n: A^MD^Arfc = r. This is achieved by setting 

where iVMDfc = Nn ■ iV„-i ■ ... ■ Nk. 

In ref. [171] also a partially improved integration scheme with multiple time scales 
was introduced, which reduces the size of the discretization errors. 

To summarize we decompose the pseudofermion action as 

S = j:S,[U] (8.17) 

k=0 

with 

So = Sg[u], s, = (grVfc, QrV^), (s.is) 

Sk = {Qk'Qk^Kl^k^Qk'Qk-M l<k<n (8.19) 

where 

Qk = Qw + iPk Pk < Pk-i (8.20) 
and Qn is the operator we are interested in: Qw if we want to perform dynamical 
simulations with Wilson fermions or Qn = Qw+ipq if we want to perform simulation 
with Wtm (/iq < pkik). The strategy is then to tune pk and n in eq. (8.17) such that 
the more expensive the computation of a certain Fk is, the less it contributes to the 
total force. The different parts of the action can then be integrated on different time 
scales Arfc chosen according to their force magnitude F^, guided by Ar^Fk = const 
for all k. 



8.5 Results 



In ref. [167,172] it was shown that this idea proves to be useful in practice. Good per- 
formances were found for the mtM-HMC (multiple time scales mass preconditioned- 
HMC) compared to the DD-HMC of ref. [163] and to a plain HMC as used in 
ref. [173]. The algorithm was tested using the standard Wilson gauge action (2.26) 
and Nf = 2 degenerate Wilson fermions (eq. 2.28 with /Xq = 0). The simulations 
were done on 24^ x 32 lattices with /3 = 5.6 and estimated pseudoscalar masses of 
mps = 665 MeV, 485 MeV, 380 MeV and 300 MeV (runs A, B, C and D). De- 
tails of the algorithm parameters as well as results for several quantities such as the 
plaquette expectation value or the vector mass my can be found in ref. [167]. 

The first important observation from this investigation is that for all four aforemen- 
tioned simulation points the preconditioning masses and time scales can be tuned 



121 



V [172] (mtM-HMC) u [163, 175] (DD-HMC) v [173] (HMC) 



A 


0.15750 


0.09(3) 


0.69(29) 


1.8(8) 


B 


0.15800 


0.11(3) 


0.50(17) 


5.1(5) 


C 


0.15825 


0.23(9) 


0.62(23) 




D 


0.15835 


~ 0.35 


0.74(18) 




Table 2 



Comparison of i' values from ref. [172], ref. [163] (with updates from [175] and ref. [173]). 
The z^- value for simulation point D is only based on an extrapolation of Ti,^^{P) in 1/mpg. 

such that simulations are stable. Examples for Monte Carlo histories of the plaquette 
expectation value or AH can be found in ref. [167]. 

In order to compare the performance of the mtM-HMC to other HMC variants 
one could choose two different measures. The first is the performance figure v = 
10~'^(2A^„ + 3)rint(-P) as introduced in ref. [163]. rint(-P) is the integrated autocorre- 
lation time of the plaquette and n is the number of integration steps for the physical 
operator Q^r necessary for one trajectory, i/ represents the number of inversions of 
the operator Qw in thousands needed in order to obtain one independent configu- 
ration. It is clearly algorithm and machine independent, but it does not account for 
the preconditioning overhead, which is, at least for the mtM-HMC, not completely 
negligible. The estimate was that this overhead was roughly a factor of two, and this 
has been recently confirmed in [174]. 

The results for the i^- values are summarized in table 2 and, while the //-values for the 
mtM-HMC and the DD-HMC are comparable, they are significantly smaller than 
the values extracted for the plain HMC algorithm used in ref. [173]. 
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Concluding remarks 



I have presented an overview of tlie tlieoretical properties and numerical results of 
Wilson twisted mass QCD. To warm up, 1 started by considering a classical theory 
with a twisted mass term (tmQCD). In the continuum a twisted mass term can 
always be rotated away by a non-anomalous change of fermion variables in the 
functional integral. 

On the lattice the twisted mass cannot be rotated away: the standard Wilson lat- 
tice action and the Wilson twisted mass (Wtm) lattice action are not related by a 
change of variables, and as a consequence have different discretization errors. The 
difference between the two discretizations is governed by the amount of disagline- 
ment, also called twist, in the chiral space between the Wilson term and the mass 
term. After describing the theoretical properties of Wtm QCD I have shown how 
the choice of working at full twist, i.e. maximally disaligning the Wilson term and 
the mass term, is rewarding in many respects. At full twist, the dimension five oper- 
ators describing the leading discretization errors in physical quantities break chiral 
symmetry in an "orthogonal" direction respect to the mass term. As a consequence 
their insertions in correlation functions of multiplicatively renormalizable operators 
vanish. Physical quantities are then automatically 0(a) improved if we stay at full 
twist at finite lattice spacing, without the knowledge of any of the improvement 
coefficients required by the implementation of the Symanzik improvement program. 
After discussing many technical issues concerning the choice of the critical mass, 
I have rederived some of the results using the so called physical basis, where the 
twist is applied to the Wilson term, keeping the mass term in the standard chiral 
direction. 

The disaglinement in chiral space between the mass term and Wilson term causes, at 
finite lattice spacing, the breaking of isospin and parity symmetry. Both symmetry 
breakings appear as O(a^) cutoff effect. The consequences of the breaking of parity 
symmetry at finite lattice spacing can be analyzed in the correlation functions with 
standard techniques. Numerical results seem to indicate that parity breaking cutoff 
effects are well under control. Isospin breaking cutoff effects induce splittings among 
flavour multiplets. In particular the neutral and the charged pions are not degenerate 
at finite lattice spacing, even if they contain degenerate quarks. The splitting is an 
O(a^) effects but numerically not negligible. I have shown that one way to partially 
mitigate this problem is to use a discretization for the valence quarks which do 
not break flavour symmetry like the Osterwalder-Seiler action or a Ginsparg- Wilson 
action. Carefully matching the renormalized quark masses for the valence and the sea 
quark action, it is then possible to eliminate the isospin splitting at the valence level. 
At the level of virtual quarks and virtual pions the problem cannot be eliminated, 
and this is certainly an issue that has to be investigated both theoretically and 
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numerically in the future. Isospin splitting in other sectors, like in baryon multiplets, 
seems to be numerically under control. 

The O(a^) cutoff effects of a Wilson-like theory deform the diagram of the chiral 
phase transition. The chiral phase diagram can be analyzed using a generalization 
of chiral perturbation theory (xPT) at finite lattice spacing to include the 0(a) and 
O(a^) cutoff effects. This generalization called W(ilson)xPT predicts two possible 
scenarios for the chiral phase diagram: the Aoki scenario, the Sharpe-Singleton sce- 
nario. There are numerical evidences that close enough to the continuum limit, with 
a Wilson lattice gauge action and Wtm lattice fermion action the Aoki scenario 
applies for the quenched model while the Sharpe-Singleton scenario applies for dy- 
namical quarks. A consequence of the Sharpe-Singleton scenario is that even if the 
twisted mass term gives a sharp infrared cutoff to the lattice theory on each gauge 
background, there is a minimal value of the twisted mass which can be simulated 
before encountering a second order phase transition point where the neutral pion 
is massless. The position of this second order point goes to zero in the continuum 
as O(a^) and depends on all the details of the lattice action used. In particular I 
have shown how using a slightly modified gauge action, like the tree-level Symanzik 
(tlSym) improved gauge action can mitigate this problem allowing simulations at 
smaller quark masses given a certain lattice spacing. In particular the tlSym gauge 
action and the Wtm fermion action allow simulations at a lattice spacing of a ~ 0.1 
fm with ^ 300 MeV. 

A further advantage of Wtm is the possibility to ease the renormalization patterns of 
phenomenologically relevant physical quantities. The pseudoscalar decay constant, 
contrary to what happens with Wilson and clover improved fermions, can be com- 
puted without the knowledge of any renormalization constant, because at full twist, 
the decay constant is related to the vector current which is protected from renor- 
malization, because of the exact flavour symmetry of Wtm in the massless limit. 
Recent results with Nj = 2 dynamical simulations show that the chiral behaviour 
in the pseudoscalar sector can be analyzed very precisely and, within the current 
understanding of all the systematic effects, it is consistent with NLO xPT. A second 
example is the provided by Bx which parametrizes the non-perturbative contribu- 
tion to the indirect CP violation in the kaon sector. Two different strategies have 
been presented which remove the mixing of the relevant four-fermion operator, both 
based on the observation that Bk can be evaluated also from a parity violating four 
fermion operator which renormalizes multiplicatively and can be extracted using 
twisted mass regularizations. 

I have concluded, in the last section, with a brief digression on recent algorithmic 
developments to simulate dynamical Wilson-like fermions. 

Wtm can be used to describe also non-degenerate quarks. I have discussed two 
different discretizations which describe non-degenerate quarks in the context of tm 
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QCD. First results with up, down, strange and charm dynamical quarks show that 
simulations of realistic QCD with realistic physical parameters are within reach of 
the current machines and algorithms. 

Although presently not all aspects of Wtm are fully investigated, Wtm is a powerful 
discretization of lattice QCD, and it certainly belongs to the pool of well founded 
fermion actions that ought to be used to control the continuum limit of physical 
quantities of interest. 

I hope that the ongoing theoretical and numerical investigations related to Wtm 
might help in the quest for solving QCD. 
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A Definitions and conventions 



A.l Index conventions 



Lorentz indices are taken from the middle Greek alphabet fi,^, . . . and run from to 
3. Latin indices k, I, . . . are used to label the components of the spatial vectors and 
run from 1 to 3. For Dirac indices we use letters a, (3, ... from the beginning of the 
Greek alphabet and they run from 1 to 4. Colour vectors in the fundamental repre- 
sentation of SU(A^c) carry indices A,B, . . . ranging from 1 to Nc, while for vectors in 
the adjoint representation capital letters a, b, . . . running from 1 to A*"^ — 1 are em- 
ployed. By abuse of notation for flavour vectors in the fundamental representation 
of SU(A'f) we use latin indices . . ranging from 1 to A'^f while for vectors in the 
adjoint representation indices a,b,. . . running from 1 to — 1 are used. It will be 
clear from the context to which case we are refering to. Repeated indices are always 
summed over unless stated and scalar products are taken with euclidean metric. 



A. 2 Dirac matrices 



We choose a chiral representation for the Dirac matrices, where 



7m 



The 2x2 matrices are taken to be 







et 



(A.l) 



with cTfc the Pauli matrices 
It is then easy to check that 




0-2 




7/. = 7i, {7^.7^ = 25 



(A.2) 



Furthermore we define 



75 = 7o7i7273 



75 
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In particular 

75 = ll ll = 1, (A.3) 

and the hermitian matrices 

f^M^ = ipf.,lu], cr^u = crl^, (A.4) 

are explicitly given by 

/ 



^cxk 



O'Ok 



\ -cTfc J y 0-fc 



(Tk 



where eijk is the totally antisymmetric tensor with €123 = 1. 
A.3 Gauge group 

A representation of the Lie algebra of SU(A'c) is given by complex x matrices 
Xab which satisfy 

= -X, tr{X} = Xaa = 0. (A.5) 

It is possible to choose a basis in this matrix space T", a = 1,2,. . . A'^^ — 1 that 
satisfies 

ti{T''T^} = (A.6) 
For Nc = 3 the standard basis is 

T« = ^, a = l,...8, (A.7) 

where A*^ denote the Gell-Mann matrices. With this convention the structure con- 
stants f"''"^, defined by 

[T^T^] = r^^T^ (A.8) 
are real and totally antisymmetric. 

A.4 Lattice derivatives 



Ordinary lattice forward and backward derivatives are diagonal in colour space and 
are defined by 

d,fix) = ^[fix + afi)-f{x)], (A.9) 
dlf{x) = -[f{x)-f{x-afi)], (A.IO) 
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where /t denotes the unit vector in direction /z. The gauge covariant derivatives 

acting on a quark fields are not trivial anymore in colour space and they are defined 
by 

V^x(a;) = -[U{x,fi)x{x + afi) - x{x)], (A.ll) 

V*x(a;) = -[x{x) -U{x- afi,^)~^x{x - afi)]. (A.12) 
The left action of the lattice derivative operators is defined by 

x{x)V^ = -[x{x + afi)U{x, fi)-' -x{x)], (A.13) 

CI 

V* = -[xix) - x{x - afi)U{x - afi, /i)]. (A.14) 



A. 5 Continuum gauge fields 



The gauge field in a continuum gauge theory belongs to the algebra of the gauge 
group and may be written as 

G,{x) = G;(x)T« (A.15) 

with real components The field strength tensor 

F^u = d^G,{x) - d,G^{x) + [G^{x), (A.16) 

can be also decomposed in the same way. The right and left covariant derivatives 
are defined by 

D^x{x) = {d, + GM^) (A.17) 

x{x)D^.= x{x)(d^-G^). (A.18) 

We stress that 9^ in formulffi(A.17,A.18) are continuum partial derivatives, and do 
not have to be confused with the lattice derivative (A. 9). 



B Symmetries in the twisted basis 



We list here the form of the relevant symmetries for a generic angle u. The SU(2) 
axial and vector twisted transformations take the form 

QTT fo\ J^^-^ - exp(-zf75r3)exp(2^r'^)exp(zf75r3)x(a;), 

exp(2^75r'^) exp(-z^r") exp(-z|75r'^). 
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err fo\ ) — ' exp(-2f 75r3) exp(i^75r«) exp(z|75r3)x(x), 
X[x) — > x(x)exp(«275r'')exp(«^75r'')exp(-?275r''). 

The twisted discrete symmetries that involve axis reflections (parity and time rever- 
sal), using the gamma matrix representation given above, are 

{f/(xo,x;0) — > f/(xo, — x; 0), f/(xo,x;A;) — > U^^{xq, —li. — ak;k), /c = l,2,3 
x(xo, x) — > 7o exp(zu;75r3)x(a;o, -x), 
x(xo, x) — > x{.Xq, -x) exp(zu;75r^)7o, 

(B.3) 

{f/(xo,x;0) — > f/"^(— Xo — a, x; 0), f/(xo,x;/c) — f/(— Xq, x; fc), /c = l,2,3 
x(xo,x) — ^ z7o75exp(zcL;75r3)x(-Xo,x), 
x(xo,x) — > -2x(-a;o,x)exp(it^75r^)757o. 

(B.4) 

Charge conjugation takes a form that is invariant under the change of variables (2.6) 

r?7(x;/i) — > U{x;i2)*, 
C: L(x)-^C-'xixV, (B.5) 
[xix) — > -xixfC, 

where C satisfies 

-lj: = C^,C~\ ^, = C^,C-\ (B.6) 
In the chosen gamma matrix representation a possible choice is C = z7o72 so that 

c = c^ = c-\ 



C Transfer matrix 



We use here the original notation of ref. [28]. The transfer matrix as an operator 
in Fock space and as an integral kernel with respect to the gauge fields has the 
structure 

To[U, U'] = fl{U)K,[U, U']fp{U'), (C.l) 
with pure gauge kernel Kq and the fermionic part 

fpiU) = det{2KBy/^exp{f]^P^Cfi) exp^-rf^^Mfi). (C.2) 

Here, the operators (x) are canonical {i is a shorthand for colour, spin and flavour 
indices) 

{r},(x),r}](y)} = 5,,a-35,y, (C.3) 
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and B and C are matrix representations of the difference operators 



k=l 

(C.4) 



C(x,y; 



Aa;BP 



E(7, 

k=l 



1 



U (x; /c)Ai?5(x + ak, y) - U (y; k)/j^S{y + ak, x) 



(C.5) 



The matrix M is defined as 



M 



2 V2k 



(C.6) 



In order to prove the positivity of the transfer matrix it is enough to show that Tp 
is bounded and invertible and that Ko[U, U'] is positive. Ko[U, U'] depends only on 
the gauge action and the proof of its positivity for the Wilson gauge action can be 
found in [28]. What remains to be proven is the positivity of B. Let's indicate with 
Ui the following unitary operator 



Ui{'il>){x) = Ui{x)ip{x + ai) 



(C.7) 



then the operator whose kernel is B (C.12) can be written as 



b = i-ky: 



Ui+U, 



We can then write the eigenvalues of this matrix in the following way 



\{B) = 1-kY. 



i=l 



2MeA(Wi) 



(C.9) 



and given the unitarity of Ui we can conclude 

1-6k< \{B) < 1 + 6k 



(C.IO) 



From this we can conclude that if |k| < | then the operator B is positive definite. 

We remark here that in principle physical positivity can be violated at finite lattice 
spacing, but one should have available a clean definition of the transfer matrix and 
the corresponding eigenvalues. Moreover positivity should be lost only at the cutoff 
scale, in order to become unimportant in the continuum limit. In ref. [176] it has 
been shown indeed that this is the case for improved gauge actions. 
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In ref. [29] it was shown that Wtm for degenerate quarks has also a well defined and 
positive tranfer matrix. The twisted mass term can be added to the antihermitian 
C matrix without changing the given proof 



3 1 . _ . 

k=i ^ 



AB 



(C.ll) 



except that /Xq must be real. 



But if we want to add a twisted term for non-degenerate quarks as given in eq. (3.20) 
one needs to add this to the hermitian matrix B. To keep the matrix B positive the 
constraint on the values of k has to be changed In fact the matrix B is given now 
by 



k=l 



SijhiYl ^(^i A;)AB5(x+aA;, y)+f/(y; A;)^]j(5(y+aA;, x) 



(C.12) 

Repeating the same argument given in eqs. (C.8- CIO) we obtain the new constraint 

, , 1 



K < 



6 + 2aen 



eq > 0. 



(C.13) 



D 0(a) improvement 



In this appendix we give more technical details of the dimension 5 operators de- 
scribing the 0(a) effects of the Wtm lattice action. In sect. 2 we have analyzed the 
form of the continuum action 5*0 using the symmetries of the lattice theory. 

We now construct Si in eq. (4.3), which contains dimension five operators. To classify 
all the possible operators, we have to use again the symmetries of the lattice action, 
and make use of partial integration in (4.5) given the space integration over y. As 
for 5*0 the residual f/v(l)3 fiavour symmetry forbids bilinears with r^'^. A number of 
potential terms can be excluded by the P = P x (/iq — /^q) symmetry, because it 
implies that parity violating dimension 5 fields have to be multiplied with a twisted 
mass term to an odd power. We can exclude then terms like mq/iqXX, "^qX^TsTsX? 
/^qXi75r3X, X-D^z75r3X and x^cr^i/-F)jjy«75r3X- The last of these, the "twisted Pauli 
term", requires a factor of fiq and thus appears only in 5*2. Charge conjugation 
symmetry excludes a term like ifiqixj^r^-f^D^x + xD^n^ilbT^x]- 
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The list of the possible fields contributing to Ci{y) then is [9,29] 



Oi = ix(y^uF^,yX. (D.l) 

O2 = mqtr{F^,F^,}, (D.2) 

O3 = mlxx, (D.3) 

O4 = m^fi^xibT^X: (D.4) 

a = /iqXX, ^ (D.5) 

Og = rriqix-f^D^x - xDf,if,x}, (D.6) 

O7 = {xD.D^x + xDfi.x}- (D.7) 



To apply the Symanzik improvement program one needs to improve also the local 
fields in the correlation functions. To be concrete we give here the example of the 
currents A^, V^^ and P'^ which appear in the Ward identities (2.18,2.19). The local 
operators contributing to the leading correction term in the effective axial current 
Al are [9,10,29]: 

{Ofi)l = Xl^—a^^D^x - xDu<yt,ulr>^X, (D-8) 
mi = X75 + x5^75 yX, (D.9) 

{Oro)l = ^qX7M75 yX, (D.IO) 

{OnT, = l^^e^^'xiJ^X , (D.ll) 
for the pseudoscalar density P"^ we need 

(O12)'' = mqX75yX , (D.12) 

and for the vector current 

(Oi3)^ = X^^u^-^D^x + xDuCr^J-^x. (D.13) 

{0,,)l = m^XlJ^X. (D.14) 
(Oi5)^ = /iqe''^' X7M75yX • (D.15) 



In general we consider a generic connected lattice correlation function ($) (4.7) at 
Xi ^ X2 ^ ■ ■ ■ ^ Xn- In the effective theory up to order a, ($) is given by eq. (4.8). 
Without contact terms the continuum equation of motions can be used in order 
to constrain the set of operators Oi. Since we consider correlation functions with 
Xi ^ X2 ^ ■ ■ ■ ^ Xni the only contact terms arise when the argument y of Ci{y) 
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is equal to the argument of one of the fields in the correlation function. But as 
we have discussed in sect. 4.1 this simply amounts to a redefinition of 0i, so the 
equations of motion of the continuum theory can be used to simplify the field basis 
for the effective action and operators. But it is important to remember that the 
coefficients appearing in the linear combination (pi will depend on the choice of the 
basis and on the value of the coefficients of Ci. 

The list of operators given above for the 0(a) lattice correction can then be reduced 
using the field equations. Formal application of the field equations gives 

Ci - C7 + 2O5 + 2C3 = (D.16) 

Ce + 2C3 + 2C4 = 0, (D.17) 
which allows to eleminate Oq and C7. 

A technical remark has to be made here. The naive relations (D.16, D.17) are ob- 
tained at tree-level of perturbation theory. At higher order they should be replaced 
by linear combinations of all basis elements with coefficients that depend on the 
coupling. Nevertheless the simple existence of such relations allows us to eliminate 
Oe and Oj. 

The same procedure can be adopted to eliminate redundant terms for the operators 
describing the leading 0(a) corrections to the pseudoscalar density and the axial 
and vector current, similarly to what we have done for the operators contributing 
to the effective action. 

For example for the axial current the operator (Os)^ can be eliminated using the 
equations of motion. 



E Automatic 0(a) improvement 



In this appendix we will give the basic idea of the first proof of automatic 0(a) 
improvement [11]. We do not give the whole proof and we invite the reader to look 
at the original paper [11] for more details. 

The form of the unimproved lattice action is S = Sq + Sp and the Wtm quark 
action Sp is defined in eq. (2.28). The form of the gauge action is not relevant 
because it has leading O(a^) discretization errors. A general Symanzik expansion 
of a connected lattice correlation function (4.8) of multiplicatively renormalized 
fields contains coefficients (e.g. in the linear combinations which define Ci) that will 
depend on the bare parameters of the lattice action. 
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We consider first the plain Wilson action setting /iq = in eq. (2.28), and we define 



here the discrete axial-vector transformatioE 



jx(a;o,x) ^ 75x(a:o,x) 

The lattice action is invariant if after a TZ^ transformation on the fields we change also 
the sign of the Wilson term and the mass term. In particular the lattice action (2.28), 
with /iq = 0, is invariant under the following transformation 

nf = 7^5 X [r ^ -r] X [mo -mo]. (E.2) 

This happens because both Wilson and mass terms are odd under the TZ^ symme- 
try. It is also interesting to note that both the Wilson and the mass term have odd 
dimensions, if we exclude the bare quark mass mo and the lattice spacing a in the 
counting. As a consequence of this, all the terms of 0(a) in the Symanzik expan- 
sion (4.8) will have opposite definite properties under r — > — r, such that averages of 
correlation functions computed on the lattice with opposite values of r have a faster 



approach to the continuum limit, i.e. the 0(a) effects are cancelled 



If the theory is tuned to full twist, u = 7r/2, this averaging procedure is automatic. 
In fact if we consider the action (5.2) the transformation r —r is equivalent to 

U UJ + 71. 



F Nucleon operators in the twisted basis 



In this appendix we show how the nucleon operators change in the twisted basis. 
This is a rather simple exercise and it is shown here just as an example how to 
proceed in a general case. The starting point is the lattice action in eq. (2.28). The 
connection between correlators in standard QCD and tmQCD can be inferred from 
classical consideration. We consider the doublet 




The axial rotation in eq. (2.6) for single flavour takes the form 

u{x) ^ exp(i'^-f5^u{x) (F.l) 

This transformation is not anomalous since is the product of 2 non-anomalous vector 
and axial transformations in the same isospin direction. 

Ideas in this direction were already put forward in [12, 177]. 
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dix) ^ exp[- t'^^,)d{x) (F.2) 
and the corresponding for the antiquark fields. The proton field 

^ = UA[ulC~^-f5dc]eABC, (F.3) 

transforms in 

— > pxn ( t ^ 



exp ( i^75 ) ^. (F.4) 



In QCD the proton correlation function 

G^(xo) = j d3x(^(xo,x)^(0,0))QCD, (F.5) 

has the general decomposition (for simplicity we consider only the leading exponen- 
tial contributions) 

P+^(-xo)A^-e*-^^-^° - P_^(xo)A^-e-^^-^°, (F.6) 

where 

^± = ^(l±7o). (F.7) 

The backward propagating contributions correspond to the antiparticles of the for- 
ward propagating states with opposite parity. The desired state is obtained using 
the appropriate projection operator P± and the appropriate direction of propaga- 
tion. Given the classical correspondence between QCD and tmQCD (2.23) the mass 
of the desired state in the proton channel can be extracted from the exponential 
behaviour of the following correlator in tm QCD 



Tr 

where 



P±G%\xo) (F.8) 



G*^(xo) = /d3x/e('^^^)^(xo,x)^(0,0)e('^^0\ . (F.9) 

J \ I tmQCD 

The correlation function in eq. (F.9) can now be computed on the lattice with the 
Wtm lattice action (2.28). 
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